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Statistical physics has proven to be a very fruitful framework to describe phenomena outside the 
realm of traditional physics. The last years have witnessed the attempt by physicists to study 
collective phenomena emerging from the interactions of individuals as elementary units in social 
structures. Here we review the state of the art by focusing on a wide list of topics ranging 
from opinion, cultural and language dynamics to crowd behavior, hierarchy formation, human 
dynamics, social spreading. We highlight the connections between these problems and other, 
more traditional, topics of statistical physics. We also emphasize the comparison of model results 
with empirical data from social systems. 
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I. INTRODUCTION 

The concept that many laws of nature are of statisti- 
cal origin is so firmly grounded in virtually all fields of 
modern physics, that statistical physics has acquired the 
status of a discipline on its own. Given its success and 
its very general conceptual framework, in recent years 
there has been a trend toward applications of statistical 
physics to interdisciplinary fields as diverse as biology, 
medicine, information technology, computer science, etc.. 
In this context, physicists have shown a rapidly grow- 
ing interest for a statistical physical modeling of fields 
patently very far from their traditional d omain of invest i- 
gations fjChakrabarti et a/-.2006 : Stauf fer et aLl . [20"06b( ) . 
In social phenomena the basic constituents are not par- 
ticles but humans and every individual interacts with 
a limited number of peers, usually negligible compared 
to the total number of people in the system. In spite 
of that, human so cieties are c harac terized by stunning 
global regularities ( Buchananl |2007( ). There are transi- 
tions from disorder to order, like the spontaneous forma- 
tion of a common language/culture or the emergence of 
consensus about a specific issue. There are examples of 
scaling and universality. These macroscopic phenomena 
naturally call for a statistical physics approach to social 
behavior, i.e., the attempt to understand regularities at 
large scale as collective effects of the interaction among 
single individuals, considered as relatively simple entities. 

It may be surprising, but the idea of a physical model- 
ing of social phenomena is in some sense older than the 
idea of statistical modeling of physical phenomena. The 
discovery of quantitative laws in the collective properties 
of a large number of people, as revealed for example by 
birth and death rates or crime statistics, was one of the 
factors pushing for the development of statistics and led 
many scientists and philosophers to call for some quan- 
titative understanding (in the sense of physics) on how 
such precise regularities arise out of the apparently er- 
ratic behavior of single individuals. Hobbes, Laplace, 
Conite, Stuart Mill and many o thers shared, to a differ- 
ent extent, this hue of thought (jBalll . 1200 j). This point 
of view was well known to Maxwell and Boltzmann and 
probably played a role when they abandoned the idea 
of describing the trajectory of single particles and in- 
troduced a statistical description for gases, laying the 
foundations of modern statistical physics. The value of 
statistical laws for social sciences has b een foreseen also 
by M ajorana in his famous tenth article ( Maioranal . ll942l 
[loos'). But it is only in the past few years that the idea 
of approaching society within the framework of statisti- 
cal physics has transformed from a philosophical declara- 
tion of principles to a concrete research effort involving a 
critical mass of physicists. The availability of new large 
databases as well as the appearance of brand new social 
phenomena (mostly related to the Internet) and the ten- 
dency of social scientists to move toward the formulation 
of simplified models and their quantitative analysis, have 
been instrumental for this change. 



In this review we mostly discuss several different as- 
pects of a single basic question of social dynamics: how 
do the interactions between social agents create order out 
of an initial disordered situation? Order is a translation 
in the language of physics of what is denoted in social 
sciences as consensus, agreement, uniformity, while dis- 
order stands for fragmentation or disagreement. It is 
reasonable to assume that without interactions hetero- 
geneity dominates: left alone, each agent would choose 
a personal response to a political question, a unique set 
of cultural features, his own special correspondence be- 
tween objects and words. Still it is common experience 
that shared opinions, cultures, languages do exist. The 
focus of the statistical physics approach to social dynam- 
ics is to understand how this comes about. The key factor 
is that agents interact and this generally tends to make 
people more similar (although many counterexamples ex- 
ist). Repeated interactions in time lead to higher degrees 
of homogeneity, that can be partial or complete depend- 
ing on the temporal or spatial scales. The investigation 
of this phenomenon is intrinsically dynamic in nature. 

A conceptual difhculty immediately arises when try- 
ing to approach social dynamics from the point of view 
of statistical physics. In usual applications, the elemen- 
tary components of the systems investigated, atoms and 
molecules, are relatively simple objects, whose behavior 
is very well known: the macroscopic phenomena are not 
due to a complex behavior of single entities, rather to 
nontrivial collective effects resulting from the interaction 
of a large number of 'simple' elements. 

Humans are exactly the opposite of such simple enti- 
ties: the detailed behavior of each of them is already the 
complex outcome of many physiological and psychologi- 
cal processes, still largely unknown. No one knows pre- 
cisely the dynamics of a single individual, nor the way 
he interacts with others. Moreover, even if one knew the 
very nature of such dynamics and such interactions, they 
would be much more complicated than, say, the forces 
that atoms exert on each other. It would be impossi- 
ble to describe them precisely with simple laws and few 
parameters. Therefore any modeling of social agents in- 
evitably involves a huge and unwarranted simplification 
of the real problem. It is then clear that any investiga- 
tion of models of social dynamics involves two levels of 
difficulty. The first is in the very definition of sensible 
and realistic microscopic models; the second is the usual 
problem of inferring the macroscopic phenomenology out 
of the microscopic dynamics of such models. Obtaining 
useful results out of these models may seem a hopeless 
task. 

The critique that models used by physicists to describe 
social systems are too simplified to describe any real sit- 
uation is most of the times very well grounded. This 
applies also to highly acclaimed models introduced by 
social s cientists, as Sch elling's model for urb an segre- 
gatio n ( Schellinj . Il97lh and Axelrod's model ( Axelrodl . 
I1997D for cultural dissemination. But in this respect, 
statistical physics brings an important added value. In 
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most situations qualitative (and even some quantitative) 
properties of large scale phenomena do not depend on 
the microscopic details of the process. Only higher level 
features, as symmetries, dimensionality or conservation 
laws, are relevant for the global behavior. With this con- 
cept of universality in mind one can then approach the 
modelization of social systems, trying to include only the 
simplest and most important properties of single indi- 
viduals and looking for qualitative features exhibited by 
models. A crucial step in this perspective is the com- 
parison with empirical data which should be primarily 
intended as an investigation on whether the trends seen 
in real data are compatible with plausible microscopic 
modeling of the individuals, are self-consistent or require 
additional ingredients. 

The statistical physics approach to social dynamics is 
currently attracting a great deal of interest, as indicated 
by the large and rapidly increasing number of papers de- 
voted to it. The newcomer can easily feel overwhelmed 
and get lost in the steadily growing flow of new publica- 
tions. Even for scholars working in this area it is difficult 
to keep up on the new results that appear at an impres- 
sive pace. In this survey we try to present, in a coherent 
and structured way, the state of the art in a wide subset 
of the vast field of social dynamics, pointing out moti- 
vations, connections and open problems. Specific review 
articles already exist for some of the topics we consider 
and we will mention them where appropriate. We aim 
at providing an up-to-date and - as much as possible - 
unified description of the published material. Our hope 
is that it will be useful both as an introduction to the 
field and as a reference. 

When writing a review on a broad, interdisciplinary 
and active field, completeness is, needless to say, a goal 
out of reach. For this reason we spell out explicitly what 
is in the review and what is not. We focus on some con- 
ceptually homogeneous topics, where the common thread 
is that individuals are viewed as adaptive instead of ratio- 
nal agents, the emphasis being on communication rather 
than strategy. A large part of the review is devoted to the 
dynamics of opinions (Sec. IIIip and to the related field of 
cultural dissemination (Sec. lIVp . Another large section 
describes language dynamics (Sec. |V]), intended both as 
the formation and evolution of a language and as the 
competition between different languages. In addition we 
discuss some other interesting issues as crowd dynam- 
ics (Sec. IVl)) . the emergence of hierarchies (Sec. IVIip . 
social spreading phenomena (Sec. lIXp . coevolution of 
states and topology (Sec. |X| and what is becoming es- 
tablished as 'human dynamics' (Sec. IVIIip . Although 
it is often very difficult to draw clear borders between 
disciplines, we have in general neglected works belong- 
ing to the field of econophysics as we ll as to evolution- 
ary theory ( Blvthe and McKanell2007l ) and evolutionary 
game theory, except for what concerns the problem of 
language formation. On such topics there a re excellent 
books and reviews (lBoucha.ud and Pottersl 120001 : iLuxl . 
I2OO7I : [Mantegna and Stanlev . ll999[) to which we refer the 



interested reader. We leave out also the physical investi- 
gation of vehicu lar traffic, a rathe r well established and 
successful field (IChowdhurv et all I2OOOI : iHelbind . I2OOII: 
'Naga tanl |2002| ). though akin to pedestrian behavior in 
crowd dynamics. The hot topic of complex networks 
has a big relevance from the social point of view, since 
many nontrivial topological structures emerge from the 
self-organization of human agents. Nevertheless, for lack 
of space, we do not to discuss such theme, for which 
we refer to (lAlbert an d Barabasi'. '2002'; B occaletti et all 
I2OO6I: iDorogovtsevan d Mcndes, 2002; Nc wmanl . l2003a " 
Networks will be considered but only as substrates where 
the social dynamics may take place. Similarly, we do not 
revi ew the intense recent activity on epidemics spre ad- 
ing ( Anderson a nd May. 'l9 9ll: iLlovd and May. '2G 0lD on 
netw orks ( Mavi . Pastor-Satorras and Vespignanil . 

I2OOII) . though we devote a section to social spreading 
phenomena. It is worth remarking that, even if we have 
done our best to mention relevant social science literature 
and highlight connections to it, the main focus of this 
work remains the description of the statistical physics 
approach to social dynamics. 

The reader will realize that, with the exception of some 
specific sections, there is a striking unbalance between 
empirical evidence and theoretical modelization, in favor 
of the latter. This does not correspond to our personal 
choice: it is a rather objective refiection of a great dis- 
proportion in the literature on social dynamics. In this 
respect things are very different from other related fields 
like complex networks or econophysics. In those cases the 
great bursts of activity occurred recently have been essen- 
tially data-driven: the availability of unprecedented huge 
new datasets has spurred first their thorough empirical 
characterization and, based on that, an intense theoret- 
ical activity followed for their modeling, well rooted in 
the comparison between theory and data. For some of 
the fields that will be reviewed below, things have gone 
the other way around. The introduction of a profusion 
of theoretical models has been mainly justified by vague 
plausibility arguments, with no direct connection to mea- 
surable facts. Very little attention has been paid to a 
stringent quantitative validation of models and theoreti- 
cal results. The contribution of physicists in establishing 
social dynamics as a sound discipline grounded on empir- 
ical evidence has been so far insufficient. We hope that 
the present review will be instrumental in stimulating the 
statistical physics community towards this goal. The lat- 
est developments of the field, especially those related to 
so-called Social Web, hint at a change of perspective in 
this direction. 



II. GENERAL FRAMEWORK: CONCEPTS AND TOOLS 

Despite their apparent diversity, most research lines 
we shall review are actually closely connected from the 
point of view of both the methodologies employed and, 
more importantly, of the general phenomenology ob- 
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served. Opinions, cultural and linguistic traits, social 
status, kinematic features are always modeled in terms of 
a small set of variables whose dynamics is determined by 
social interactions. The interpretation of such variables 
will be different in the various cases: a binary variable 
will indicate yes/no to a political question in opinion dy- 
namics, two synonyms for a certain object in language 
evolution, two languages in competition, whether some- 
body has been reached by a rumor or not, etc.. Other 
details may differ, but often results obtained in one case 
can immediately be translated in the context of other 
sub- fields. In all cases the dynamics tends to reduce the 
variability of the initial state and this may lead to an 
ordered state, with all the agents sharing the same fea- 
tures (opinion, cultural or linguistic traits, velocity) or 
to a fragmented (disordered) state. The way in which 
those systems evolve can thus be addressed in a unitary 
way using well known tools and concepts from statisti- 
cal physics. In this spirit some of the relevant general 
questions we will consider in the review include: What 
are the fundamental interaction mechanisms that allow 
for the emergence of consensus on an issue, a shared cul- 
ture, a common language, collective motion, a hierarchy? 
What favors the homogenization process? What hinders 
it? 

Generally speaking the drive toward order is provided 
by the tendency of interacting agents to become more 
alike. This effect is often t ermed 'social influenc e' in 
the social science literature (jFestinger et 'dl. Il95rf> and 
can be seen as a counterpart of ferromagnetic interaction 
in magnets. Couplings of anti- ferromagnetic type, i.e., 
pushing people to adopt a state different from the state 
of their neighbors, are also in some cases important and 
will be considered. 

Any modelization of social agents inevitably neglects a 
huge number of details. One can often take into account 
in an effective form such unknown additional ingredi- 
ents assuming the presence of noise. A time-independent 
noise in the model parameters often represents the vari- 
ability in the nature of single individuals. On the other 
hand a time-dependent noise may generate spontaneous 
transitions of agents from one state to another. A crucial 
question has then to do with the stability of the model 
behavior with respect to such perturbations. Do sponta- 
neous fluctuations slow down or even stop the ordering 
process? Does diversity of agents' properties strongly af- 
fect the model behavior? 

An additional relevant feature is the topology of the in- 
teraction network. Traditional statistical physics usually 
deals with structures whose elements are located regu- 
larly in space (lattices) or considers the simplifying hy- 
pothesis that the interaction pattern is all-to-all, thus 
guaranteeing that the mean field approximation is cor- 
rect. This assumption, often also termed homogeneous 
mixing, generally permits analytical treatment, but it is 
hardly realistic in a social context. Much more plausible 
interaction patterns are those denoted as complex net- 
works (see Sec. III. Bp . The study of the effect of their 



nontrivial topological properties on models for social dy- 
namics is a very hot topic. 

One concept playing a special role in many social dy- 
namic models and having no equally common counter- 
part in traditional statistical physics is 'bounded confi- 
dence', i.e., the idea that in order to interact two individ- 
uals must be not too different. This parallels somewhat 
the range of interaction in physics: if two particles are 
too far apart they do not exert any influence on each 
other. However let us stress that the distance involved in 
bounded confidence is not spatial, rather being defined 
in a sort of opinion space. We will discuss in the review 
several instances of this general principle. 

Let us finally clarify some problems with nomencla- 
ture. Being a strongly interdisciplinary field, in social dy- 
namics there is a natural tendency towards a rather free 
(or sloppy) use of terms. This heterogeneity is in some 
cases very confusing as it happens for some words (like 
polarization) that have been used with opposite meaning. 
For the sake of clarity we specify that in this review con- 
sensus indicates a configuration where all agents share 
the same state. When many possible states are possi- 
ble but only two of them survive in the population we 
speak of polarization. Fragmentation indicates instead a 
configuration with more than two surviving states. 



A. Order and disorder: the Ising paradigm 

In the previous section we have seen that the com- 
mon theme of social dynamics is the understanding 
of the transition from an initial disordered state to a 
configuration that displays order (at least partially). 
Such type of trans i tions abound in traditi onal statisti- 
cal physics (|Huand . Il987l : iKubo et all Il985t) . It is worth 
summarizing some important concepts and tools used in 
that context, as they are relevant also for the investiga- 
tion of social dynamics. We will illustrate them using 
a paradigmatic example of order-disorder transitions in 
physics, t he one exhibi t ed by the Ising model for ferro- 
magnets ( Binnev et all Il992| ). Beyond its relevance as 
a physics model, the Ising ferromagnet can be seen as 
a very simple model for opinion dynamics, with agents 
being influenced by the state of the majority of their in- 
teracting partners. 

Consider a collection of N spins (agents) that can 
assume two values ±1. Each spin is energetically pushed 
to be aligned with its nearest neighbors. The total energy 
is 



H 



(1) 



where the sum runs on the pairs of nearest-neighbors 
spins. Among the poss ible types of dynam i cs, th e most 
common (Metropolis) ( Landau and Bindeil . l2005l ) takes 
as elementary move a single spin flip that is accepted with 
probability exp(— AE'/fc^T'), where /S.E is the change in 
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FIG. 1 Snapshots of equilibrium configurations of the Ising 
model (from left to right) below, at and above Tc. 



energy and T is the temperature. Ferromagnetic inter- 
actions in Eq. ([1]) drive the system towards one the two 
possible ordered states, with all positive or all negative 
spins. At the same time thermal noise injects fluctua- 
tions that tend to destroy order. For low temperature T 
the ordering tendency wins and long-range order is estab- 
lished in the system, while above a critical temperature 
Tc the system remains macroscopically disordered. The 
transition point is characterized by the average magne- 
tization m = i/NJ2i{^i) passing from for T > Tc to 
a value m(T) > for T < Tc- The brackets denote the 
average over different realizations of the dynamics. This 
kind of transitions is exhibited by a wealth of systems. 
Let us simply mention, for its similarity with many of 
the social dy namic mo dels discussed in the review, the 
Potts model (|wiil . ll982D . where each spin can assume one 
out of q values and equal nearest neighbor values are en- 
ergetically favored. The Ising model corresponds to the 
special case q = 2. 

It is important to stress that above Tc no infinite-range 
order is established, but on short spatial scales spins are 
correlated: there are domains of -t-1 spins (and others 
of —1 spins) extended over regions of finite size. Below 
Tc instead these ordered regions extend to infinity (they 
span the whole system), although at finite temperature 
some disordered fluctuations are present on short scales 
(Fig.©. 

Not only the equilibrium properties just described, 
that are attained in the long run, are interesti ng. A much 
investigated and nontrivial issue (|Bravl . ll994D is the way 
the final ordered state aX T < Tc is reached, when the 
system is initially prepared in a fully disordered state. 
This ordering dynamics is a prototype for the analogous 
processes occurring in many models of social dynamics. 
On short time scales, coexisting ordered domains of small 
size (both positive and negative) are formed. The sub- 
sequent evolution occurs through a coarsening process of 
such domains, which grow larger and larger while their 
global statistical features remain unchanged over time. 
This is the dynamic scaling phenomenon: the morphol- 
ogy remains statistically the same if rescaled by the typ- 
ical domain size, which is the only relevant length in the 
system and grows over time as a power-law. 

Macroscopically, the dynamic driving force towards 
order is surface tension. Interfaces between domains 
of opposite magnetization cost in terms of energy and 
their contribution can be minimized by making them as 



straight as possible. This type of ordering is often re- 
ferred to as curvature-driven and occurs in many of the 
social systems described in this review. The presence of 
surface tension is a consequence of the tendency of each 
spin to become aligned with the majority of its neighbors. 
When the majority does not play a role, the qualitative 
features of the ordering process change. 

The dynamic aspect of the study of social models re- 
quires the monitoring of suitable quantities, able to prop- 
erly identify the buildup of order. The magnetization of 
the system is not one of such suitable quantities. It is 
not sensitive to the size of single ordered domains, while 
it measures their cumulative extension, which is more or 
less the same during most of the evolution. The appro- 
priate quantity to monitor the ordering process is the 
correlation function between pairs of spins at distance 
r from each other, C{r,t) — {si{t)si+r{t)) — {si{t))'^, 
where brackets denote averaging over dynamic realiza- 
tions and an additional average over i is implicit. The 
temporal variable t is measured as the average number of 
attempted updates per spin. The dynamic scaling prop- 
erty implies that C (r, t) is a function only of the ratio 
between the distance and the typical domain size L{t): 
C{r, t) = L{tYF[r/ L{t)]. L{t) grows in time as a power- 
law t^/^ . The dynamic exponent z is universal, indepen- 
dent of microscopic details, possibly depending only on 
qualitative features as conservation of the magnetization 
or space dimensionality. In the Glauber/Metropolis case 
z = 2 in any dimension. Another quantity often used is 
the density of interfaces na{t) = Na{t)/Np, where Np is 
the total number of nearest neighbor pairs and Na the 
number of such pairs where the two neighbors are in dif- 
ferent states: na = 1/2 means that disorder is complete, 
while Ha — indicates full consensus. 

Finally, a word about finite size effects. The very con- 
cept of order-disorder phase-transitions is rigorously de- 
fined only in the limit of a system with an infinite number 
of particles (thermodynamic limit), because only in that 
limit truly singular behavior can arise. Social systems are 
generally composed by a large number N of agents, but 
by far smaller than the number of atoms or molecules in a 
physical system. The finiteness of N must play therefore 
a cr ucial ro le in the analysis of models of social dynam- 
ics (jToral a nd Tessonc, 2007). Studying what happens 
when N changes and even considering the large- A'^ limit 
is generally very useful, because it helps characterizing 
well qualitative behaviors, understanding which features 
are robust, and filtering out non-universal microscopical 
details. 



B. Role of topology 

An important aspect always present in social dynamics 
is topology, i.e., the structure of the interaction network 
describing who is interacting with whom, how frequently 
and with which intensity. Agents are thus supposed to 
sit on vertices (nodes) of a network, and the edges (links) 
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define the possible interaction patterns. 

The prototype of homogeneous networks is the un- 
correlated random graph model proposed by Er dos and 
Renyi (ER model) jErdos and Rfevill959lll960f) . whose 



construction consists in drawing an (undirected) edge 
with a fixed probability p between each possible pair out 
of N given vertices. The resulting graph shows a bino- 
mial degree distribution, the degree of a node being the 
number of its connections, with average (fc) ~ Np. The 
degree distribution converges to a Poissonian for large 
N. If p is sufficiently small (order 1/N), the graph is 
sparse and presents locally tree-like structures. In order 
to account for degree heterogeneity, other constructions 
have been propose d for random g;raphs with arbitrary de- 
gree distributions (I Aiello and Lu ^ 20011 Catanzaro et all 
\200A iGoh et al. I. l200ll:lMollov and Reedl . ll995lll998D . 

A well-known paradigm, especially for social sciences, 
is that of "small- world" networks, in which, on the 
one h and, the averag e distance between two agents is 
small ( Milgranilll967l ). growing only logarithmically with 
the network size, and, on the other hand, many trian- 
gles are present, unlike ER graphs. In order to recon- 
cile both properties, Watts and Stro ^atz have introduced 
the small-world network model (Wa tts and Strogat j . 

[1998) , which allows to interpolate between regular low- 
dimensional lattices and random networks, by introduc- 
ing a certain fraction p of random long-r ange connections 
into an initially regular lattice ( Ne wman and Wattd . 

[1999) . In (jWatts and Strogati [19981 ) two main quanti- 
ties have been considered: the characteristic path length 
L{p), defined as the number of edges in the shortest path 
between two vertices, averaged over all pairs of vertices, 
and the clustering coefficient C{p), defined as follows. If 
a node i has k connections, then at most fc(fc— l)/2 edges 
can exist between its neighbors (this occurs when every 
neighbor of i is connected to every other neighbor). The 
clustering coefficient C{p) denotes the fraction of these al- 
lowable edges that actually exist, averaged over all nodes. 
Small-world networks feature high values of C'{p) and low 
values of L{p). 

Since many real networks are not static but evolv- 
ing, with new nodes entering and establishing connec- 
tions to already existing nodes, many models of growing 
networks have also been introduced. The Bara basi and 
Albert model (BA) (jBarabasi and Alberlll999l ). has be- 
come one of the most famous models for complex hetero- 
geneous networks, and is constructed as follows: starting 
from a small set of m fully interconnected nodes, new 
nodes are introduced one by one. Each new node selects 
TO older nodes according to the preferential attachment 
rule, i.e., with probability proportional to their degree, 
and creates links with them. The procedure stops when 
the required network size N is reached. The obtained 
network has average degree (fc) = 2m, small clustering 
coefficient (of order 1/A^) and a power law degree dis- 
tribution P{k) ~ k~'^, with 7 = 3. Graphs with power 
law degree distributions with 7 < 3 are referred to as 
scale-free networks. 



An extensive analysis of the existing network models 
is out of the scope of the present review and we refer the 
reader to the huge literature on the s o-called complex 
netw o rks ( Albert and Barabasi. 20021: Boccaletti et al. . 
2006t ICaldarellL 2007: Dorogovtsev and Mended. 12005 : 



Newman. 2003 



!in. zuu/: uorogovtsev ana ivienaeg . uuudt 
3 iPastor-Satorras and Vespignanil . |2004 ). 



It is nevertheless important to mention that real net- 
works often differ in many respects from artificial net- 
works. People have used the social network metaphor 
for over a century to represent complex sets of relation- 
ships between members of social systems at all scales, 
from interpersonal to international. A huge amount of 
work has been carried out about the so-called social net- 
work analysis (SNA), especially in the soc i al science liter- 
ature dFreeman', '2004': Granovetter, '1973', 1983; 'Moreno], 
I1934 IScott. 2000; Wasserman and Faust . ,1994). Re- 
cently the interest of physicists triggered the investiga- 
tion of many different networks : from the network of 
scienti fic colla borations ((Barabasi et all [2002; Newman , 
2001 al ibi. I2OO4I ) to that of sexual contacts dLilieros et al. , 
2OOII) and the ongoing social relat ionships (iHolme 



2003 ), from email exchanges networks (lEbel et all 2002 ; 



Eckmann et all |2004 [Newman etai]. 120(321) to the dat- 
ing community network ( Holme et all |2004|) a nd to 



mobile communi cation networks ( Onnela et all 120071 : 
iPalla et aZI . l2007l ). just to quote a few examples. From 
this experimental work a set of features characteriz- 
ing soci al networks have been id entified. It has been 
shown ( Newman and ParkL I2003D how social networks 
differ substantially from other types of networks, namely 
technological or biological. The origin of the difference is 
twofold. On the one hand they exhibit a positive corre- 
lation between adjacent vertices (also called assqrtativ- 
ity), while most non-social networks ([Newman! . l2003bl : 
[Pastor-Satorras et Ml, l200lD are disassortative. A net- 
work is said to show assortative mixing if nodes with 
many connections tend to be linked to other nodes with 
high degree. On the other hand social networks show 
clustering coefficients well above those of the correspond- 
ing random models. These results opened the way to a 
modeling activity aimed at reproducing in an artificial 
and controlled way the same fe atures observed in real so- 
cial networks ( Jin et all [200 1[ ) . We cannot review here 
all these attempts but we have quoted some relevant ref- 
erences all along the review when discussing specific mod- 
eling schemes. It is important to keep in mind that future 
investigations on social dynamics will be forced to take 
into account in a more stringent way st ructural and d; 



mto account m a more strmgent way st ructural ana gy - 
namic properties of real social networks ( Roehneil . [2003) • 



When applying models of social dynamics on spe- 
cific topologies several non-trivial effects may arise, po- 
tentially leading to important biases for the dynamics. 
For instance, on networks with strongly heterogeneous 
degree distributions the dynamics of models with bi- 
nary asymmetric interaction rules, i.e., where the two 
selected agents have different roles, may be affected 
by the order in which the interaction partners are se- 
lected. This is a consequence of the fact that, in a 
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network with degree distribution -P(fc), if one picks a 
random node and then one of its neighbors, it is hkely 
that the second node has a higher degree than the first. 
More precis ely, the degree distribution of t he second is 
kP{k)/{k) (jPastor-Satorras and Vespignanil . 12004) . For 
scale-free degree distributions this may have strong ef- 
fects (as is the case of the voter model, as seen in 
Sec. nil. Bl and of the Naming Game, as seen in Sec. IV. Bp . 



C. Dynamical systems approach 

One of the early contribution of physicists to the study 
of social systems has been the introduction of methods 
and tools coming from the theory of dynamical systems 
and non-linear dynamics. Thi s development goes un- 
der the name of socio dynamics { Helbind . Il99ll : IWeidlichl . 
119911 : IWeidlic"hl . l2002D . The term sociodynamics has been 
introduced to refer to a systematic approach to mathe- 
matical modeling in the framework of soci al sciences. 

Sociodynamics is a branch of synergetics (lHakenl. [l978l ) 
devoted to social systems, featuring a few important dif- 
ferences. In synergetics one typically starts with a large 
set of microscopic equations for the elementary compo- 
nents and performs a reduction of the degrees of free- 
dom. This is not the case for social systems, for which 
equations at the microscopic level are not available. In 
this case one has to identify relevant macro- variables and 
construct directly equations for them, based on reason- 
able and realistic social hypotheses, i.e., informed by so- 
cial driving forces. The typical procedure consists in 
defining probabilistic transition rates per unit of time for 
the jumps between different configurations of the system 
corresponding to different values of the macro- variables. 
The transition rates are used as building blocks for set- 
ting up the equations of motion for the probabilistic evo- 
lution of the set of macro- variables. The central evolution 
equation in sociodynamics is the master equation, a phe- 
nomenological first-order differential equation describing 
the time evolution of the probability P(m,i) for a sys- 
tem to occupy each one of a discrete set of states, defined 
through the set of macro- variables m: 



dP(m,t) _ ^ 



[W^'.mPK, t) - I^™.™'P(m, t)l (2) 



where Wm,m' represents the transition rate from the state 
m to the state m'. The master equation is a gain- loss 
equation for the probability of each state m. The first 
term is the gain due to transitions from other states m', 
and the second term is the loss due to transitions into 
other states m'. 

While it is relatively easy to write down a master equa- 
tion, it is quite another matter to solve it. It is usually 
highly non-linear and some clever simplifications are of- 
ten needed to extract a solution. In general only nu- 
merical solutions are available. Moreover, typically the 



master equation contains too much information in com- 
parison to available empirical data. For all these reasons 
it is highly desirable to derive from the master equation 
simpler equations of motion for simpler variables. One 
straightforward possibility is to consider the equations of 
motion for the average values of the macro-variables m, 
defined as: 



nikit) = ^mfcP(m, i). 



(3) 



The exact expression for the equations of motion for 
mk{t) does not lead to simplifications because one should 
already know the full probability distribution P{m,t). 
On the other hand, under the assumption that the dis- 
tribution remains unimodal and sharply peaked for the 
period of time under consideration, one has: 



P(m,O^P(m(i)), 



(4) 



yielding the approximate equations of motions for Wik (t) , 
which are now a closed system o f cou pled differential 
equations. We refer to ( Weidlichl . |2002| ) for a complete 
derivation of these equations as well as for the discus- 
sion of several applications. The appr oach has also be en 
applied to mo del behavioral changes ( Helbind . Il993al lbl: 
lHelbing| . [T99l . 



D. Agent-based modeling 

Computer simulations play an important role in the 
study of social dynamics since they parallel more tra- 
ditional approaches of theoretical physics aiming at de- 
scribing a system in terms of a set of equations, to be 
later solved numerically and/or, whenever possible, an- 
alytically. One of the most successful methodologies 
used in social dynamics is agent-based modeling. The 
idea is to construct the computational devices (known 
as agents with some properties) and then simulate them 
in parallel to model the real phenomena. In physics 
this technique can be traced back to m olecular dynam- 
ics ( Alder and Wainwrightl. 1957 . 19591) an d Metropolis 
and Monte Carlo (jMetropolis et all Il953^ simulations. 
The goal is to address the problem of the emergence 
from the lower (micro) level of the social system to the 
higher (macro) level. The origin of agent-based model- 
ing can be traced back to the 1940s, to the introduction 
by Von N eumann and Ulam of the not ion of cellular au- 
tomaton (|Neumannl . [19661 : lUlaml . [l96Q) . e.g., a machine 
composed of a collection of cells on a grid. Each cell can 
be found in a discrete set of states and its update oc- 
curs on discrete time steps according to the state of the 
neighboring cells. A well-known example is Conway's 
Game of Life, defined in terms of simple rules in a vir- 
tual world shaped as a 2-dimensional checkerboard. This 
kind of algorithms became v ery popular in population 
biology (|Matsuda et aLl . ll992l ). 
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The notion of agent has been very important in 
the de velopment of the concept of A rtificial Intelli- 
ffence (|McCarthvl Il959l : iMinskvl . Il96lh . which tradi- 
tionally focuses on the individual and on rule-based 
paradigms inspired by psychology. In this framework the 
term actors was used to indicate interactive objects char- 
acterized by a certain number of internal states, acting 
in parallel and exchanging messages (Hewitt, 1973). In 
computer science the notion of actor turned in that of 
agent and more emphasis has been put on the interac- 
tion level instead of autonomous actions. 

Agent-based models were primarily used for social sys- 
tems by Reynolds, who tried to model the reality of living 
bio logical agents, k nown as artificial life, a term coined 
in ( Langtonf . Il996l ). Reynolds introduced the notion of 



individual-based models, in which one investigates the 
global consequences of local interactions of members of a 
population (e.g. plants and animals in ecosystems, vehi- 
cles in traffic, people in crowds, or autonomous charac- 
ters in animation and games). In these models individual 
agents (possibly heterogeneous) interact in a given envi- 
ronment according to procedural rules tuned by charac- 
teristic parameters. One thus focuses on the features of 
each individual instead of looking at some global quantity 
averaged over the whole population. 

The Artificial Life community has been the 
first in developing agent-based models 

Steeld . 



19911: iMever and Wilsonl Il990l: 



(Maes 



Varela and Bourgind . 11995 IWeisT [T999h . but 



19951 



smce 



then agent-based simulations have become an important 
tool in other scientif ic fields a n d in p a rticular in th e study 
of social systeni s (lAxelrod . 2006t IConte et all 1997 : 
Macy and Willed. 120021; [Sc hweitzeJ, [20031; Imoldrjdge, 
2002[ ). In (|Epstein and A^t cU. 1990), by focusing on a 
bottom-up approach, the first large scale agent model, 
the Sugarscape, has been introduced to simulate and 
explore the role of social phenomena such as seasonal 
migrations, pollution, sexual reproduction, combat, 
trade and transmission of disease and culture. In 
this context it i s worth m e ntioni ng the concept of 
Brownian agent (Schweitzed. |2003[ ) which generalizes 
that of Brownian particle from statistical mechanics. 
A Brownian agent is an active particle which possesses 
internal states, can store energy and information and 
interacts with other agents through the environment. 
Again the emphasis is on the parsimony in the agent 
definition as well as on the interactions, rather than on 
the autonomous actions. Agents interact either directly 
or in an indirect way through the external environment, 
which provides a feedback about the activities of the 
other agents. Direct interactions are typically local 
in time and ruled by the underlying topology of the 
interaction network (see also Sec. III.BI) . Populations 
can be homogeneous (i.e., all agents being identical) 
or heterogeneous. Differently from physical systems, 
the interactions are usually asymmetrical since the role 
of the interacting agents can be different both for the 
actions performed and for the rules to change their 



internal states. Agent-based simulations have now 
acquired a central role in modeling complex systems 
and a large literature has been rapidly developing in 
the last few years about the internal structure of the 
agents, their activities and the multi-agent features. An 
exhaustive discussion of agent-based models is out of the 
scope of the present review, but we refer to ((_Schweitzer|, 
|2003| ) where the role of active particles is thoroughly 
discussed with many examples of applications, ranging 
from structure formation in biological systems and 
pedestrian traffic to the simulation of urban aggregation 
or opinion formation processes. 



III. OPINION DYNAMICS 
A. Introduction 

Agreement is one of the most important aspects of so- 
cial group dynamics. Everyday life presents many situa- 
tions in which it is necessary for a group to reach shared 
decisions. Agreement makes a position stronger, and am- 
plifies its impact on society. 

The dynamics of agreement /disagreement among in- 
dividuals is complex, because the individuals are. Sta- 
tistical physicists working on opinion dynamics aim at 
defining the opinion states of a population, and the el- 
ementary processes that determine transitions between 
such states. The main question is whether this is possi- 
ble and whether this approach can shed new light on the 
process of opinion formation. 

In any mathematical model, opinion has to be a vari- 
able, or a set of variables, i.e., a collection of numbers. 
This may appear too reductive, thinking about the com- 
plexity of a person and of each individual position. Ev- 
eryday life, on the contrary, indicates that people are 
sometimes confronted with a limited number of posi- 
tions on a specific issue, which often are as few as two: 
right/left, Windows/Linux, buying/selling, etc. If opin- 
ions can be represented by numbers, the challenge is to 
find an adequate set of mathematical rules to describe the 
mechanisms responsible for the evolution and changes of 
them. 

The development of opinion dynamics so far has been 
uncoordinated and based on individual attempts, where 
social mechanisms considered reasonable by the authors 
turned into mathematical rules, without a general shared 
framework and often with no reference to real sociolog- 
ical studies. The first opinion dynami cs designed b y 
a physicist was a model proposed in (Weidlichl . Il97lt ). 
The model is based on the probabilistic framework of 
sociodynamics, which was discussed in Sec. III. CI Later 
on the Ising model made i ts first appearance in opin- 
ion d ynamics ( Galam et aLl . ll98^lGalam and Moscovicil 
Il99ll ). The spin-spin coupling represents the pairwise in- 
teraction between agents, the magnetic field the cultural 
majority or propaganda. Moreover individual fields are 
introduced that determine personal preferences toward 
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either orientation. Depending on the strength of the in- 
dividual fields, the system may reach total consensus to- 
ward one of the two possible opinions, or a state where 
both opinions coexist. 

In the last decade, physicists have started to work ac- 
tively in opinion dynamics, and many models have been 
designed. We focus on the models that have received 
more attention in the physics literature, pointing out 
analogies as well as differences between them: the voter 
model (Sec. IIII.Bp . majority rule models (Sec. IIII.C|1 . 
models based on social impact theory (Sec lIILDjl . Sznajd 
fSec. IIILEj) and bounded confidence models (Sec. IIILFj) . 
In Sec. IIII.GI other models are briefly discussed. Finally, 
in Sec. IIII.Hl we review recent work that aims at an em- 
pirical validation of opinion dynamics from the analysis 
of data referring to large scale social phenomena. 



B. Voter model 

1. Regular lattices 

The voter model has been named in this way for the 
very natural interpretation of its rules in terms of opin- 
ion dynamics; for its extremely simple definition, how- 
ever, the model has been thoroughly investigated also in 
fields quite far from social dynamics, like probability the- 
ory and popul ation genetics. Voter dynam ics was first 
considered in ( Clifford and Sudburvl . I1973D as a model 
for the competition of s pecie s and named "voter model" 
in ( HoUev and LiggetU 119751 ). It has soon become pop- 
ular because, despite being a rather crude description of 
any real process, it is one of the very few non-equilibrium 
stochasti c processes tha t can be solved exactly in any di- 
mension (lRedneil . [200ll) . It can also be seen as a model 
for dimer-dimer heterogeneous catalys is in the reaction 
controlled limit ( Evans and "Ravl . ll993D . 

The definition is extremely simple: each agent is en- 
dowed with a binary variable s = ±1. At each time step 
an agent i is selected along with one of its neighbors j 
and Si = Sj, i.e., the agent takes the opinion of the neigh- 



bor. This update rule implies that agents imitate their 
neighbors. They feel the pressure of the majority of their 
peers only in an average sense: the state of the majority 
does not play a direct role and more fluctuations may be 
expected with respect to the zero-temperature Glauber 
dynamics. Bulk noise is absent in the model, so the states 
with all sites equal (consensus) are absorbing. Starting 
from a disordered initial condition, voter dynamics tends 
to increase the order of the system, as in usual coarsening 
processes ( Scheucher and Spohnl . Il988l ). The question is 
whether full consensus is reached in a system of infinite 
size. In one-dimensional lattices the dynamics is exactly 
the same of the zero-temperature Glauber dynamics. A 
look at the patterns generated in two-dimensional lattices 
(Fig. (2) indicates that domains grow but interfac es are 
very rough, at odds with usual coarsening systems ( Bravl 
I1994D . 




FIG. 2 Evolution of a two-dimensional voter model starting 
from a dropl et (top) or a fully d isordered configuration (bot- 
tom). From (jDornic et adl2(30ll ). 



Earl y studies, per formed by proba- 

bilists (IClifford and Sudbur^. Il973l: I Cox and Griffeathl . 



1X9861: iHollev and Liggettl . Il975l: iLiggettj . I1985D . have 
exploited the fact that the model can be exactly mapped 
on a model of random walkers that coalesce upon en- 
counter. This duality property allo ws to use the powerfu l 
machinery of random walk theory ( Liggettl . Il985l Il999f ) . 
We prefer to follow another derivat ion of the general 
solution on lattices (jFrac hcbourg an d Krapivskvl [l996). 
based on earlier work ifkrapivskvl . 119921 ). Considering 
a d-dimensional hypercubic lattice and denoting with 
S — {si} the state of the system, the transition rate for 
a spin k to flip is 



Wk{S) = W{s 



(5) 



where j runs over all 2d nearest neighbors and the prefac- 
tor, setting the overall temporal scale, is chosen for con- 
venience. The probability distribution function P{S, t) 
obeys the master equation 



dt 



P{S,t) = [VFfc(5'=)P(5^t) - Wk{S)P{S,t)] , (6) 



where S'^ is equal to S except for the flipped spin Sk- The 
linear structure of the rates ([5]) has the nice consequence 
that the equations for correlation functions of any order 
{sk ■ ■ ■ si) = J2s ^i^^ ■ ■ ■ Si can be closed, i.e., they 
do not de pend on higher-order functio ns and hence can 
be solved (jScheucher and Spohnl . 1 19881) . 

The equation for the one-body correlation function is 



dt 



(7) 



where Ak is the discrete Laplace operator. Summing 
over k one sees that the global magnetization (s) — 
l/-^X]fc(*fc) conserved. This conservation immedi- 
ately allows to determine the probability that a finite 
system will end up with all spins up or down (exit prob- 
ability), depending on the initial density of up spins 
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p(0) = ((s) + l)/2. This gives P„p(p(0)) = p(0) in any 
dimension. 

The two-body correlation function obeys 



^(sfcS;) = (Afe + Ai){skSi) 



(8) 



The structure of this equation, as well as of those for 
higher-order correlation functions, is similar in any di- 
mension to the equations for correlators of the one- 
dimensiona l Ising m odel with zero-temperature Glauber 
dynamics ( GlaubeJ . [1963) and can be solved analo- 
gously, via Laplace transform. In this way the asymp- 
totic behavior of the density of active interfaces na(t) — 
(sfcSfc+i))/2 is derived (,Frachebourg and Krapivskvi . 
Il99fil ^ 



^-{2-d)/2 ^ < 2 

na{t) - { l/ln(t) d = 2 
a - bt-'^/^ d>2. 



(9) 



Eq. ([9]) shows that for d < 2 the voter model undergoes 
a coarsening process leading to complete consensus. For 
d > 2 instead, it exhibits asymptotically a finite density 
of interfaces, i.e., no consensus is reached (in an infinite 
system) and domains of opposite opinions coexist indef- 
initely in time. In terms of duality the lack of order in 
high dimensions is a consequence of the transient nature 
of random walks in d > 2: diffusing active interfaces have 
a finite probability to meet and annihilate. For d = 2 the 
exact expression of the density of active interfaces for 
large times is 



la{t) - 



21n(i) +ln(256) 



(10) 



The large constant value in the denominator of Eq. (fIT 
makes the approach to the asymptotic logarithmic decay 
very slow, and explains why different la ws were hypoth- 
esized, based on numerical evide nce ( Evans and Ravi . 
119931: iMeakin and Scalapinol . Il987l) . 

Beyond the expression for the density na{t), the 
solution of Eq. ([8|) allows to write do wn a s caling 
form for the correlati o n fun ction C (r, t) (jDornid . 1 19981 : 
IScheucher and Spohnl . Il988l ). In d = 2 the solution 
violates logarithmically the standard scaling for m (see 
Sec. I II. Ap holding for usual coarsening phenomena ( Bravl 
ll994f K This violation reflects the fact that there is 
no single domain length in the system. On the con- 
trary, domains of all sizes t ", with < a < 1/2 



are simultaneously present ( Cox and Griffeathl . Il986l : 
IScheucher and Spohnl . [19881 ). 

While in an infinite system consensus is reached only 
if d < 2, in any dimension consensus is invariably 
reached asymptotically if the system is finit e. The tim e 
Tjv needed depends on the system size N dCo^ . I1989D : 
Tn ~ N"^ for d = 1, Tat ~ TVlnA^ for d = 2, while 
Tjv '--^ for d > 2. It is worth remarking that the way 
consensus is reached on finite systems has a completely 
different nature for d < 2 (where the system coherently 



tends towards order by coarsening) and for d > 2 (where 
consensus is reached only because of a large random fluc- 
tuation). 

Expression ^ for the spin-flip rates is rather 
special. How much of the voter behavior is re- 
tained if rates are modified ? A natural general- 
izatio n ( Drouffe and Godrechel Il999l : Ide Oliveira et all 
Il993l ) considers transition rates of the form Wk{S) = 
1/2[1 — Skfk(S)], where fk{S) is a local function with 
|/fe('S')| < 1. A local dynamics that is spatially symmetric 
and preserves the up-down symmetry requires fk{S) to 
be an odd function of the sum of the nearest neighbors. In 
a square lattice, the local field can assume five values, but 
only two of them are independent, /(2) = — /(— 2) = x, 
and /(4) — — /(— 4) — y. Voter dynamics corresponds 
to a: = 1/2 and y = 1, while x — y corresponds to the 
majority-vote model (Sec. IIII.C|) . y = 2a;/(l -t- x^) gives 
the transition rates of Glauber dynamics, and the case 
y — 2x corresponds to the noisy voter model (see be- 
low). The significance of the two parameters is straight- 
forward: y gauges bulk noise, i.e., the possibility that a 
spin fully surrounded by equal spins flips to the oppo- 
site position; The value y — I implies absence of such 
noise. The parameter x instead measures the amount of 
interfacial noise. Simulations and a pair approximation 
treatment show that the phase-diagram of this general- 
ized model is divided in a ferromagnetic region around 
the a; = 1,?; = 1 point (zero-temperature Glauber dy- 
namics) and a paramagnetic phase, separated by a line 
of continuous phase-transitions terminating at the voter 
model point. Changing the interfacial noise parameter 
X, while keeping y — I, one finds a jump of the order 
parameter, indicating a first-order transition. Hence the 
voter point is critical, sitting exactly at the transition 
between order and disorder driven by purely interfacial 
noise. 

More physical insight is provided by consider- 
ing a droplet of u p spi ns surrounded by negative 
spins (|Dornic et all 1200 ll ). Th e Ca h n-AU en theory 
for curvature-driven coarsening (|Bravl . Il994f l predicts 
in d = 2 a linear decay in time of the droplet 
area, the rate being proportional to surface tension. 
In the voter model instead, the interface of the 
droplet r oughens but its radius rernains s t atistically un- 
chang ed ( Dall'Asta and CastellanoL 120071 : iDornic et all 
I2OOID . showing that no surface tension is present (Fig. [2]). 

The phase diagr am of the generalized 

model ( de Oliveira et all Il993( ) described above seems 
to suggest that the voter model is rather peculiar, 
being only a point in the line of continuous transi- 
tions between the ferromagnetic and the paramagnetic 
phase. However, voter-like behavior (characterized by 
the absence of surface tension leading to logarithmic 
ordering in d = 2) can be found in other models. It has 
been argued ( Dornic et all 1200 ll ) that voter behavior is 
generically observed at order-disorder non-equilibrium 
transitions, driven by interfacial noise, between dynam- 
ically symmetric absorbing states. This symmetry may 
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be enforced either by an up-down symmetry of the local 
rules or by global conservation of the magnetization. 
The universal exponents associated to the transition are 
(3 = 0, and = 1/2 in all dimensions, while 7 — 1/2 for 
d = 1, and 7 = 1 for d > 2 with logarith mic corrections 
at the upper cri t ical d imension d = 2 ( Dornic et al\ . 
I200ll:lde Oliveiral . I2003D . 

The original voter dynamics does not in- 
clude the possibility for a spin to flip sponta- 
neously when equal to all its neighbors. The 
noisy voter model (iGranovsky and Madras! . Il995l : 
IScheucher and Spohnl. 19881) . also called linear Glauber 
model (|de Oliveiral . 12003^ includes this possibility, via 
a modification of the rates ([5]) that keeps the model 
exactly solvable. The effect of bulk noise is to destroy 
long-range order: the noisy voter model is always 
in the paramagnetic phas e of the generalized model 
of ( de Oliveira et a!] . Il993t) . so that domains form only 
up to a finite correlation length. As the strength of 
bulk noise is decreased, the length grows and the voter 
first-order transition occurs for zero noise. 

The investigation of the generalized voter univer- 
sality class and its connections with other classes of 
non-equilibrium phase transitions is a complicated and 
open issue, approached also via fiel d -theoretical meth- 
ods (iDickman and Tretvakovl Il995l : iDroz et 
iHammal et al\ . \200^ . 



2. Modifications and applications 

Being a very simple non-equilibrium dynamics with a 
nontrivial behavior, the voter model has been investi- 
gated with respect to many properties in recent years, 
including persistence, aging and correlated percolation. 
Furthermore, many modifications of the original dynam- 
ics have been proposed in order to model various types 
of phenomena or to test the robustness of the voter phe- 
nomenology. A natural extension is a voter dynamics for 
Potts variables (multitype voter model), where many of 
the res ults obtained for the Ising case are easily general- 
izable ( Sire and Maiumdail . Il995f ). 

One possible modification is the presence of quenched 
disorder, in the form of one "zealot" , i.e., an individ- 



ual that does not change its opinion (iMobilial . Hool). 
This modification breaks the conservation of magneti- 
zation: in d < 2 the zealot influences all, inducing 
general consensus with its opinion. In higher dimen- 
sions consensus is still not reached, but in the neigh- 
borhood of the zealot the stationary state is biased to- 
ward his op inion. The case of man y zealo t s has also been 
addre ssed ( Mobilia and Georgievl 120051 : iMobilia et all 

iocri). 

Anot her variant is the constrained voter 

model ( Vazquez et all l2003l ). where agents can be 
in three states (leftists, rightists, or centrists) but 
interactions involve only centrists, while extremists do 
not talk to each other. In this way a discrete analogue of 



bounded confidence is implemented. Detailed analytical 
results give the probabilities, as a function of the initial 
conditions, of ending up with full consensus in one of 
the three states or with a mixtu re of the extr e mists , 
with little chan ge between d = 1 (I Vazquez e<: al. I. I2OO3D 
and mean field ( Vazquez and Redneii 2004). A similar 
mode l with three states is the AB-model ( Castello et al\ . 
I2OO6D . Here the state of an agent evolves according to 
the following rules. At each time step one randomly 
chooses an agent i and updates its state according to 
the following transition probabilities: 

PA^AB = 1/2(75, PB^AB = l/2(7A, (H) 

PAB^B = 1/2(1 -a^), PAB^A = l/2{l-aB)ll2) 

where ct; (^=A,B,AB) are the local densities of each state 
in the neighborhood of i. The idea here is that, in order 
to go from A to B one has to pass through the inter- 
mediate state AB. At odds with the constrained voter 
model, however, here extremes do interact, since the rate 
to go from state A to AB is proportional to the density 
of neighbors in state B. This implies that consensus on 
the AB state or a frozen mixture of A and B is not pos- 
sible, the only two possible absorbing states being those 
of consensus of A or B type. 

In the original voter model the decision on the possible 
flip of a spin depe nds on one neighbor. In the "v acillat- 
ing" voter model (jLambiotte and Rednej . I2007D a site 
checks the state of two of its neighbors and flips if ei- 
ther is different from himself. This leads to a bias to- 
wards the zero magnetization state, anti-coarsening be- 
havior and consensus time scaling exponentially with 
th e system size N. A r e lated one-dimensional model 
is (jLambiotte and Rednerl . |2007| ) . 

Another modification is the introduction of memory in 
the f orm of noise reduction ( Dall'Asta and Castellanol . 
I2OO7D . Each spin has associated two counters. When an 
interaction takes place with a positive (negative) neigh- 
bor, instead of modifying the spin the positive (negative) 
counter is increased by one. The spin is updated only 
when one of the counters reaches a threshold r. This 
change induces an effective surface tension, leading to 
curvature-driven coarsening dynamics. 



A counterintuitive effect is reported in (jStark et al\ . 
I2OO7D where a voter model with "inertia" is presented: 
the local rate of change of a site decreases linearly with 
the time since the last fiip,. up to a finite satura- 
tion value. On two-dimensional lattices, Watts-Strogatz 
and fully connected networks, a sufficiently weak inertia 
makes global dynamics faster: consensus time is smaller 
than with fixed rate. 

Other variants of the original voter model have been 
devised for studying ecological problems. Recent pub- 
lications in the physics literature are the study of di- 
versity in_glaii^communities (voter model with speci- 
ation (|Zillio et all [2005)), or the investigation of fixa- 
tion in the evolution of c ompeting species (biased voter 
model (jAntal et a/.l . l2006t) ). 
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3. The voter model on networks 

Non-regular topologies have nontrivial effects on the 
ordering dynamics of the voter model. 

On a complete graph the Fokker-Planck equation for 
the probability density of the magnetization has the form 
of a one-dimensional diffusion equation with a position- 
depen dent diffusion constant, and can be solved analyt- 
ically ( Slanina and Lavickal [200^. The lack of a drift 
term is the effect of the lack of surface tension in the 
model dynamics. The average time needed to reach con- 
sensus in a finite system can be computed exactly for any 
value of the initial magnetization and scales as the size 
of the system N. The tail of the distribution can also be 
computed and has an exponential decay exp(— i/iV). 

When considering disordered topologies different ways 
of defining the voter dynamics, that are perfectly equiva- 
lent on regular lattices, give rise to noncquivalent gener- 
alizations of the voter model. When the degree distribu- 
tion is heterogeneous, the order in which a site and the 
neighbor to be copied are selected does matter, because 
high-degree nodes are more easily chosen as neighbors 
than low-degree vertices. 

The most natural generalization {direct voter model) 
is to pick up a site and make it equal to one of its neigh- 
bors. In this way one of the fundamental properties of 
the voter model , conservation of the global magnetiza - 
tion, is violated (jSuchecki et aLl . l2005atTWu et a/.l . l2004) . 
To restore conservation a link-update dynamics must be 
considered (jSuchecki et ai, 2005a): a link is selected at 
random and then one node located at a randomly cho- 
sen end is set equal to the other. If instead one chooses 
first a node and copies its variable to a randomly selected 
neighbo r one obtains the reverse (or invasion) voter dy- 
namics (ICastellanoll2005D . 

On highly heterogeneous substrates these different def- 
initions result in different behaviors. The mean consen- 
sus time Tjv has been computed in ( Sood et~ai] . l2007t 
ISood and Redneil . |2005( ) for the direct voter dynamics 
on a generic graph, by exploiting the conservation of a 
suitably defined degree-weighted density uj of up spins. 



-N^ [(1 

tJ'2 



w)ln(l - w) -)- wlnw)] , (13) 



where fik is the fc-th moment of the degree-distribution. 
For networks with scale-free distributed degree (with ex- 
ponent 7), T/v scales then as for 7 > 3 and sublin- 
early fo r 7 < 3, in good agreement with ntimerical simu- 
lations (IC astehano ajj, |2005|; ISood and Redneil . 120051 : 
ISuchecki et ai, . ,2005a). The same approach gives, for 
the other versions of voter dynamics on graphs, a lin- 
ear dependence of the consensus time on N for link- 
update dynamics (independent of the degree distribu- 
tion) and Tjv ^ N for any 7 > 2 for the reverse- 
voter dynamics, again i n good agre e ment with simula- 
tions lICastellanoL 120051: ISood eToII . I2007D . A general 
analysis of v oter-like dynamics o n generic topologies is 
presented in ( Baxter et aLl . [20081) . with particular refer- 




FIG. 3 Log-log plot of the fraction tia of active bonds be- 
tween nodes with different opinions. Empty symbols are for 
the one-dimensional case (p = 0). Filled symbols are for 
rewiring probability p — 0.05. Data are for A'^ = 200 (cir- 
cles), iV = 400 (squares), N = 800 (diamonds), N = 1600 
(triangles up) and A'^ = 3 200 (triangles left). R eprinted fig- 
ure with permission from (jCastellano et 0/1120031 ). Copyright 
2003 from EDP Sciences. 



ence to applications in population genetics and biodiver- 
sity studies. 

Another interesting effect of the topology occurs 
when v oter dynamics is con s idered on small-world net- 
works (|Watts and Strogatd . Il998[ ). After an initial 
regime equal to the one-dimensional behavior, the den- 
sity of active interfaces forms a plateau (Fig. [3]), because 
shortcuts hinder their diffusive motion. The system re- 
mains trapped in a metastable state with coexisting do- 
mains o f opposite opinions, whose typical length scales 



as l /jp ( Castellano et all 120031 : IVilone and Castellanol . 
l2004f ). p being the fraction of long-range connections. 



The lifetime of the metastable state scales with the 
linear system size L so that for finite systems consensus 
is eventually reached on a temporal scale shorter than 
on a regular one-dimensional lattice (L^)- For infinite 
systems instead, the state with coexisting opinions is ac- 
tually stable, leading to the conclusion that long-range 
connections prevent the complete ordering of the voter 
model, i n a way similar to what occurs for Glauber dy- 
namics ()Bover and Mir amontesl . 2003). A general dis- 
cussion of the interplay between topolo gy and dynam- 
ics for the voter model is presented in iSuchecki et all . 
l2005bD . A comparison between the behavior of the voter 
dy namics and the AB- model on modular networks is 
in (jCastello et a/.l . l2007t) . 
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FIG. 4 MR model. The majority opinion inside a discussion 
group (here of size three) is taken by all agents. 



C. Majority rule model 

In a population of N agents, endowed with binary opin- 
ions, a fraction p-^- of agents has opinion +1 while a frac- 
tion j9_ = 1 — opinion —1. For simplicity, suppose 
that all agents can communicate with each other, so that 
the social network of contacts is a complete graph. At 
each iteration, a group of r agents is selected at ran- 
dom (discussion group): as a consequence of the inter- 
action, all agents take the majority opinion inside the 
group (Fig. 2]). This is the basic principle of the ma- 
jority rule (MR ) model , whic h was proposed to describe 
public debates ((Galaml . [2003 ). 

The group size r is not fixed, but is selected at each 
step from a given distribution. If r is odd, there is al- 
ways a majority in favor of either opinion. If r is even, 
instead, there is the possibility of a tie, i.e., that either 
opinion is supported by exactly r/2 agents. In this case, 
one introduces a bias in favor of one of the opinions, say 
-|-1, and that opinion prevails in the group. This pre- 
scription is inspired by the principle of social inertia, for 
which people are reluctant to accept a reform if there is 
no cl ear majority in its favor ( Friedman and Friedmanl . 
Il984l ) . Majority rule with opinion bias was originally ap- 
plied within a simple model describing hierar chical voting 
in a society (iGalaml . [19861 Il990l [Tm l2000l) . 

Defined as the initial fraction of agents with the 
opinion +1, the dynamics is characterized by a threshold 
Pc such that, for > pc {p'\. < Pc), all agents will have 
opinion +1 (—1) in the long run. The time to reach 
consen sus (in number of u pdates per spin) scales like 
logA^ (jTessone et al If the group sizes are odd, 

Pcix) = 1/2, due to the symmetry of the two opinions. If 
there are groups with r even, pc < 1/2, i.e., the favored 
opinion will eventually be the dominant one, even if it is 
initially shared by a minority of agents. 

The MR model^ with a fixed group size 
r was analytically solved in the mean field 
limit ( Krapivskv and Rednerj . l2003t) . The group 
size r is odd, to keep the symmetry of the two opinions. 



^ Th e name Majority R ule model was actually coined 

in llKrapivskv and Rednerl . l2003h . S ince th i s mo del is just 
a special case of the one introduced in l|Galaml . l2003 ). we adopt 
this name since the beginning of the section. 



The solution can be derived both for a finite population 
of N agents and in the contin uum limit ol N ^ o o. Th e 
latter derivation is simpler ( Chen and RedneJ . [200I), 
and is sketched here. 

Let Sfe — ±1 be the opinion of agent fc; the av- 
erage opinion (magnetization) of the system is m = 
1 /N Sk = p+ — P-. The size of each discussion group 
is 3. At each update step, the number iV+ of agents in 
state -I- increases by one unit if the group state is + H — , 

while it decreases by one unit if the group state is H . 

One thus has: 



dN+^3{plp--p+pl) 



1), (14) 



where the factor of 3 is due to the different permutations 

of the configurations + H — and H . Eq. (fT4| can be 

rewritten as: 

dN+ N . 1 X . 

p+ = -2p+{p+ - -){p+ - 1), (15) 



N 3 



with the time increment dt — 3/iV, so that each agent 
is updated once per unit of time. The fixed points are 
determined by the condition p^ — and from Eq. (fTS]) 
we see that this happens when p+ — 0, 1/2 and 1, re- 
spectively. The point p+ = 1/2 is unstable, whereas the 
others are stable: starting from any p+ ^ 1/2, all agents 
will converge to the state of initial majority, recovering 
Galam's result. The integration of Eq. ([T5)) yields that 
the consensus time grows as log N. 

On a c?-dimensional lattice, the discussion group is 
localized around a randomly chosen lattice site. In 
one dimension, the model is not analytically solvable. 
Since the average magnetization is not conserved by 
the MR dynamics, the exit probability, i.e., the prob- 
ability that the final magnetization is -1-1, has a non- 
trivial dependence on the initial magnetization in the 
thermodynamic limit and a minority can actually win 
the contest. Consensus time grows as N^. In higher 
dimensions ( Chen and Rednen . [2005), the dynamics is 
characterized by diffusive coarsening. When the ini- 
tial magnetization is zero, the system may be trapped 
in metastable states (stripes in 2d, slabs in 3d), which 
evolve only very slowly. This leads to the existence of 
two distinct temporal scales: the most probable con- 
sensus time is short but, when metastable states ap- 
pear, the time needed is exceedingly longer. As a con- 
sequence, the average consensus time grows as a power 
of N, with a dimension-dependent exponent. When 
the initial magnetization is non-zero, metastable states 
quickly disappear. A crude coarse-graining argument re- 
produces qualitatively the occurrences of metastable con- 
figurations for any d. Numerical simulations show that 
the MR model in four dimensions does not reproduce 
the results of the mean field limit, so the upper critical 
dimension of the MR model is larger than four. The 
MR dynamics was also inves tigated on netwo rks with 
strong degree heterogeneities ( Lambiott^ . I2OOII) and on 
networks with community structure, i.e., graphs consist- 
ing of groups of nodes with a comparatively large density 
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of internal links with respe ct to the density of connections 
between differen t group s (|Lambiotte and AusloosL [2OO7I: 
iLambiotte et aLl . l2007atl. The MR model was studied on 
small world lattices ( Li et al\ . l2006l ) as well. 

The MR model has be en extended to mu l ti-stat e opin- 
ions and plurality rule ( Chen and Redneil . |2005[ ). The 
number of opinion states and the size of the interaction 
groups arc denoted with s and G, respectively. In the 
mean field limit, the system reaches consensus for any 
choice of s and G, in a time that scales like log N , as in 
the 2-state MR model. On a square lattice, if the number 
of states s is too large, there are no groups with a major- 
ity, so the system does not evolve, otherwise the evolution 
is based on diffusive coarsening, similarly to that of the 2- 
state MR model. Again, two different timescales emerge 
when s is small, due to the existence of metastable states. 
When s and G approach a threshold, there is only one do- 
main that grows and invades all sites, so there is only one 
time scale. The plurality rule is a special extension of the 
MR rule when there are more than two opinion states: in 
this case, all agents of a group take the opinion with the 
most representatives in the group. The evolution leads 
to consensus for any s and G, because all interaction 
groups are active (there is always a relative majority); 
when the opinions reduce to two, the dynamics becomes 
identical to that of the 2-state MR model, so there will 
be metastable states and two different timescales. 

Modifications of the MR model include : a model where 
agents can move in space ( Galam et al\ . l2002t IStauffed . 
l2002al ): a dynamics where eac h agent interacts wi th a 
variable number of neighbor s (iTessone et al\. 20041) : an 
extension to three opinions (|Gekle et all l2005li ~ the in- 
troduction of a probability to favor a particular opinion, 
that could v ary arn o ng diff erent individuals and/or so- 
cial groups (|Galaml . laOOSa); the presence of "contrari- 
ans", i.e., agents that initially take the majority opin- 
ion in a group discussion, but that right after t h e dis- 
cussion switch to the opposite opinion ( Galaml . 120041 : 
IStauffer and Sa Martin s. 2004); th e presence of one-sided 
contrarians and unsettled agents ( Borghesi and Galaml . 



I2OO6D : the pr esence of inflexible agents , that always stay 
by their side ( Galam and Jacobi . l2007t ) . 



We now discuss some variants of the ma- 
jority r ule. In the m ajority -minority (MM) 

model ( Mobiha and Redneil . l2003l ). one accounts 
for the possibility that minorities take over: in a 
discussion group the majority opinion prevails with a 
probability p, whereas with a probability 1 — p it is the 
minority opinion that dominates. For a discussion group 
of three agents, the magnetization m changes by an 
amount 2p — 4(1 — p) at each interaction, which means 
that, for p — pc — 2/3, m does not change on average, 
like in the voter model. In the mean field limit, the 
model can be solved analytically: the exit probability 
turns out to be a step function for p > pc, (i.e., the 
system will evolve towards consensus around the initial 
majority opinion), whereas it equals 1/2 for p < pc, 
which means that the system is driven towards zero 



magnetization. 

Another interesting mo del based on m ajority rule is 
the majority- vote model ( Liggettl . Il985l ). At each up- 
date step, with a probability 1 — q a spin takes the 
sign of the majority of its neighbors, with a probabil- 
ity q it takes the minority spin state. If there is a 
tie, the spin is flipped with probability 1/2. The pa- 
rameter q is the so-called noise parameter. We stress 
that a single spin is updated at each time step, at 
variance with the MR model. For q = the model 
coincides with the Ising m odel with zero-temperature 
Glauber kinetics (jGlaubeii . Il963f) . On a regular lattice, 
the majority- vote model presents a phase transition from 
an ordered to a disorde red state at a crit ical value qc 
of the noise parameter ( de Oliveiral Il992l ). The criti- 
cal exponents of the transition are in the Ising univer- 
sality class. Recent studies showed that the majority- 
vote model also generates an order-disorder phase transi- 
tion on small- wo rld lattices (Campos et al. 2003.) and on 
random graphs ( Lima et aU 2008rPereira and Moreiral 
I2OO5D . 

In a recent model, an agent is convinced if there is 
at least a fraction p of its n eighbors sharing the same 
opinion ( Klimek et all l2007fl . This model interpolates 
between the rule where majority dominates {p = 1/2) 
and the unanimity rule {p = 1), where an agent is in- 
fluenced by its neighbors only if they all have the same 
opinion ( Lambiotte et a/.l . l2007bl ). Another model based 
similar to majority-vote is p resented an d studied on di- 
rected networks in ([Sanchez et al 



D. Social impact theory 



Th e psychological theory of social impact ( Latanel 
Il98ll ) describes how individuals feel the presence of their 
peers and how they in turn influence other individuals. 
The impact of a social group on a subject depends on the 
number of the individuals in the group, on their convinc- 
ing power, and on the distance from the subject, where 
the distance may refer both to spatial proximity or to the 
closeness in an abstract space of personal relationshigs. 
The original cellular automata intr o duced in (jLatang . 
ll98lD and refined in (jNowak et all Il990l ) represent a 
class of dynamic models of statistical mechanics, that are 
exact ly solvable in the mean field limit ( Lewenstein et all . 
I1992D . 

The starting point is a population of N individuals. 
Each individual i is characterized by an opinion ai — ±1 
and by two real-valued parameters, that estimate the 
strength of its action on the others: persuasiveness pi 
and supportiveness Si, that describe the capability to 
convince someone to change or to keep its opinion, re- 
spectively. These parameters are assumed to be random 
numbers, and introduce a disorder that is responsible for 
the complex dynamics of the model. The distance of a 
pair of agents i and j is dij. In the simplest version, 
the total impact li that an individual i experiences from 
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his/her social environment is 



N 



N 



^(1 

3 = 1 '1 



(16) 



where a > 2 expresses how fast the impact decreases 
with the distance da between two individuals. The first 



term of Eq. ()16|) expresses the persuasive impact, i.e. the 
pressure exerted by the agents with opposite opinions, 
which tend to enforce an opinion change; the second term 
instead is the supportive impact, i.e. the pressure exerted 
by the agents with the same opinion of z, which favor the 
status quo. In both cases, the impact of each agent on i 
is proportional to its persuasiveness/supportiveness. 
The opinion dynamics is expressed by the rule 



ai{t + 1) = -~sgn[ai{t)Ii{t) + 



(17) 



where hi is a random field representing all sources other 
than social impact that may affect the opinion (e.g. mass 
media) . According to Eq. p?]) , a spin flips if the pressure 
in favor of the opinion change overcomes the pressure to 
keep the current opinion {li > for vanishing hi). 

For a system of fully connected agents, and without 
individual field s, the model presents infi nitely many sta- 
tionary states ( Lewenstein et all Il992l ). The order pa- 
rameter of the dynamics is a complex fu nction of one 
variable, hke in spin glasses (jMezard et all . ll987) . 

In general, in the absence of individual fields, the dy- 
namics leads to the dominance of one opinion on the 
other, but not to complete consensus. If the initial mag- 
netization is about zero, the system converges to config- 
urations characterized by a large majority of spins in the 
same opinion state, and by stable domains of spins in 
the minority opinion state. In the presence of individ- 
ual fields, these minority domains become metastable: 
they remain stationary for a very long time, then they 
suddenly shrink to smaller clusters, which again persist 
for a very long time, before shrinking again, and so on 
("staircase dynamics"). 

The dynamics can be modified to account for other 
proc esses related t o social behavior, such as learn- 
ing ( Kohrin j . Il996l ). the response of a population to 
the simultaneous acti on of a st r ong leader and exter- 
nal influenc e (Holyst aLl . [2000l : lKacperski and Holvstl . 
Il996l . Il997l [l999l l206ol ) and the mitigation of social 
impact due to the coexistence of d ifferen t individuals 
in a group ( Bordogna and Albanol . |2007[ ). For a re- 
view of_statisticalniechanical models of social impact, 
see (|Holvst et aZ.l . l200ll ). 

Social impact theory neglects a number of realistic fea- 
tures of social interaction: the existence of a memory 
of the individuals, which reflects the past experience; a 
finite velocity for the exchange of information between 
agents; a physical space, where agents have the possibil- 
ity to move. An important extension of social impact the- 
ory that includ es those features is based on active Brow- 
nian particles ( Schweitzeil l2003t ISchweitzer and Holvstl . 



I2OOOI) . that are Brownian particles endowed with some 
internal energy depot that allows them to move and to 
perform several tasks as well. The interaction is due 
to a scalar opinion field, expressing the social impact 
of all agents/opinions at each point in space; the par- 
ticles/agents act as sources of the field and are in turn 
affected by it, both in opinion and in space. Each agent i 
is labeled by its opinion ai = ±1 and its personal strength 
Si. The field of opinion cr, at position r and time t, is in- 
dicated with hai^Tt). The transition probability rates 
w{a'\a), for an agent to pass from opinion a to opinion 
ct', with (7 7^ (t', are defined as: 



w{a'\a) = rje^p{[h^,{r,t) - h^{r,t)]/T}, 



(18) 



where T is a social temperature. The dynamics is ex- 
pressed by two sets of equations: one set describes the 
spatio-temporal change of the opinion field 



d 



N 



— /iCT(r, t) SiS^^„^S{r - Ti) - jh„{r, t) + DhAh„{r, t), 

i—l 

(19) 

the other set presents reaction-diffusion equations for the 
density na-{r, t) of individuals with opinion a, at position 
r and time t 



d 



r{v,t) = -VK(r,t)aV/i<,(r,i)] +D„A7i^(r,t) 
-^[w{cT'\a)n„{r,t) - w{cT\cj')n„,{v,t)]. (20) 

a' ^(7 



In the equations above, N is the number of agents, I/7 
is the average lifetime of the field, Dh is the diffusion 
constant for information exchange, I?„ the spatial diffu- 
sion coefficient of the individuals, a measures the agents' 
response to the field. The three terms on the right-hand 
side of Eq. (fT9)) represent the microscopic density of the 
agents' strength, the relaxation of the field (modeling the 
memory effects) and its diffusion in space, respectively. 
The first two terms on the right-hand side of Eq. ([^0]) 
indicate the change in the agents' density due to their 
motion in space, whereas the sum expresses the balance 
between the gain and loss of individuals due to opinion 
changes. 

Eqs. and ([^D]) are coupled: depending on the lo- 
cal intensity of the field supporting either opinion, an 
agent can change its opinion, or migrate towards loca- 
tions where its opinion has a larger support. Opinion 
changes and migrations have a non-linear feedback on 
the communication field, which in turn affects the agents, 
and so on. The model presents three phases, depending 
on the values of the parameters: a paramagnetic phase, 
where both opinions have the same probability (1/2) of 
being selected at every place (high-temperature, high- 
diffusion), a ferromagnetic phase, with more agents in fa- 
vor of one opinion over the other (low-temperature, low- 
diffusion), and a phase in which either opinion prevails 
in spatially separated domains (segregation). 
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E. Sznajd model 

In the previous section we have seen that the im- 
pact exerted by a social group on an individual in- 
creases with the size of the group. We would not pay 
attention to a single guy staring at a blank wall; in- 
stead, if a group of people stares at that wall, we may 
be tempted to do the same. Convincing somebody is 
easier for two or more people than for a single indi- 
vidual. This is the ba s ic p rinciple behind the Sznajd 
model (jStauffeil . l2003al: [sTn aid-Weron. 20 05bl). In its 
original version ( Sznaid-Weron and Szna^ 120001 ) . that 
we call Sznajd B, agents occupy the sites of a linear chain, 
and have binary opinions, denoted by Ising spin variables. 
A pair of neighboring agents i and i -I- 1 determines the 
opinions of their two nearest neighbors i — 1 and i + 2, 
according to these rules: 
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if Si = Si+i, then Si_i = Si 
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(21) 



FIG. 5 Sznajd model. In the most common version of the 
model (Sznajd A), a pair of neighboring agents with the same 
opinion convince all their neighbors (top), while they have no 
influence if they disagree (bottom) . 



if s,j 7^ Si+i, then s,j_i = s^+i and s.^ 
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(22) 



So, if the agents of the pair share the same opinion, they 
successfully impose their opinion on their neighbors. If, 
instead, the two agents disagree, each agent imposes its 
opinion on the other agent's neighbor. 

Opinions are updated in a random sequential order. 
Starting from a totally random initial configuration, 
where both opinions are equally distributed, two types 
of stationary states are found, corresponding to con- 
sensus, with all spins up (m = 1) or all spins down 
(to = —1), and to a stalemate, with the same number of 
up and down spins in antiferromagnetic order (to = 0). 
The latter state is a consequence of rule ((22|) . that fa- 
vors antiferromagnetic configurations, and has a prob- 
ability 1/2 to be reached. Each of the two (ferromag- 
netic) consensus states occurs with a probability 1/4. 
The values of the probability can be easily deduced from 
the up-down symmetry of the model. The relaxation 
time of the system into one of the possible attractors 
has a log-normal distribution ( Behera and Schweitzeii 
I2003D . The number of agents that never changed opin- 
ion first decays as a power law of time, and then it 
reaches a constant but finite value, at o dds with the Ising 
model ( Stauffer and de Oliveiral [20021 ). The exit proba- 
bility has been calculated analy tically for both random 
and correlated initial cond itions ( Lambiotte and Redneil . 
120071 : ISlanina et o/.l . l2007tl . 

Since the very introduction of the Sznajd model, it has 
been argued that a distinctive feature of its dynamics is 
the fact that the information flows from the initial pair 
of agents to their neighbors, at variance with the other 
opinion dynamics models, in which instead agents are in- 
fluenced by their neighbors. Because of that the Sznajd 
model was supposed to describe how opinions spread in 
a soc iety. On the other hand, in ( Behera and Schweitzeii 
I2003D it has been shown that, in one dimension, the di- 
rection of the information flow is actually irrelevant, and 



that the Sznajd B dynamics is equivalent to a voter dy- 
namics. The only difference with the classic voter model 
is that an agent is not influenced by its nearest neighbors 
but by its next-to-nearest neighbors. Indeed, the dy- 
namics of Sznajd B on a linear chain can be summarized 
by the simple sentence "just follow your next-to-nearest 
neighbor" . The fact that in Sznajd a pair of agents is 
updated at a time, whereas in the voter model the dy- 
namics affects a single spin, introduces a factor of two in 
the average relaxation time of the equivalent voter dy- 
namics; all other features are exactly the same, from the 
probability to hit the attractors to the distributions of 
decision and relaxation times. Therefore, Sznajd B does 
not respect the principle of social validation which moti- 
vated its introduction, as each spin is influenced only by 
a single spin, not by a pair. 

Sznajd rule (j22p is unrealistic and was soon replaced 
by alternative recipes in subsequent studies. In the most 
popular alternative, that we call Sznajd A, only the ferro- 
magnetic rule (HD) holds, so the neighbors of a disagree- 
ing agents' pair maintain their opinions. Extensions of 
the Sznajd model to different substrates usually adopt 
this prescription and we shall stick to it unless stated 
otherwise. On the square lattice, for instance, a pair of 
neighboring agents affect the opinions of their six neigh- 
bors only if they agree (Fig. [5]). In this case, the exit 
probability is a step function with threshold at to = 0: 
if the initial magnetization m < 0, the system always 
attains consensus with to = — l;ifTO>0 initially, the 
steady state is consensus with m = 1. The distribu- 
tion of the times required to reach complete consensus 
is broad, b ut not a log- normal like for Sznajd B in one 
dimension ( Stauffer et aLl . [2000l ). We stress that Sznajd 
B in one dimension has no phase transition in the exit 
probability, due to the coexistence of ferro- and antifer- 
romagnetic stationary states. 

The fixed points of Sznajd A dynamics hold if one 
changes the size of the pool of persuading agents. The 
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only exception is represented by th e so-called Ochrom bel 
simplification of the Sznajd model ( Ochrombe]ll200il) . in 
which a single agent imposes its opinion on all its neigh- 
bors. 

The results m entioned above were derived by computer 
simulations. In ( Slanina and Lavickal l20d3[ ) an exact so- 
lution for a Sznajd-like dynamics on a complete graph 
has been given. Here a pair of randomly selected agents 
i and j interacts with a third agent fc, also taken at ran- 



dom. If Si = s,-, then Sk 



Sj, otherwise nothing 



happens. The evolution equation for the probability den- 
sity P{m, t) that the system has magnetization m at time 
t reads: 

^^P{m,t) = -^[{\~m^)mP{m,t)]. (23) 

Eq. (j23p is derived in the thermodynamic limit and it 
represents a pure drift of the magnetization. The general 
solution is: 



P(m,i) = [(1 -m2)m]-V(e 



(24) 



where the function / depends on the initial conditions. 
If P(m, t = Q) = 5{ni — too)i i-e., the system starts with 
a fixed value mp of the magnetization, P(rn, t) is a S- 
function at any moment of the evolution; the center is 
pushed by the drift towards the extremes -f 1 if mo > 
or —1 if mo < 0, which are reached asymptotically. So, 
the initial magnetization toq determines the final state 
of the system, which is consensus, and there is a phase 
transition when mo changes sign. Eq. (|23p also allows 
to derive the behavior of the tail of the distribution of 
the times to reach the stationary states of the dynamics, 
which turns out to be exponential. 

Some effort has been devoted to find a proper Hamil- 
tonian for mulation of Sznajd dynamics I Sznaid-Weronl . 
I2OO2I [2004 . 2005a). It turns out that the rules of the 
model are equivalent to the minimization of a local func- 
tion of spin-spin interactions, the so-called disagreement 
function. On a linear chain of spins, the disagreement 
function for spin i reads: 



Ei = -JlSiSi+l - J2SiS 



i+2, 



(25) 



where Ji and J2 are coupling constants, whose values 
determine the type of dynamics, and i + 1, i + 2 are the 
right nearest and next-to- nearest neighbors of i. Here, 
spin i takes the value that minimizes Ei . The function Ei 
and its minimization defines the Two-Component (TC) 
model (Sznaid-Weroni . ,20021) . We remark that, when 
J1J2 > 0, the two terms of Ei are equivalent, so only 
one can be kept. Sznajd B dynamics is recovered for 
— J2 < Ji < J2, J2 > 0, but the model has a much richer 
behavior. Based on the values of the pair of parameters 
Ji and J2, one distinguishes four phases, delimited by the 
bisectors Ji ± J2 = 0. Besides the known ferro- and anti- 
ferromagnetic attractors, a new stationary configuration 
emerges, with pairs of aligned spins whose signs alternate 



(... + H hH ...). The TC mode l has b een extended 

to the square lattice (ISznaid-WeronL I2004D. and can be 
exactl y solved in the mean field limit ( Sznaid-Weronl . 
l2005aD . In general, we stress that the model is not equiv- 
alent to a Hamiltonian model at zero temperature, be- 
cause it is not possible to define a global energy for the 
system. The sum of the disagreement function Ei over 
all spins does not play the role of the energy: the local 
minim ization of Ej can lea d to an increase of its global 
value (ISznaid-Weronl . l2004D . 

Sznajd dynamics turns out to be a special cas e of the 
genera l sequential probabilistic model (GPM) (jCalaml . 
l2005bl ). Here, opinions are Ising spins: the proportions 
of both opinions at time t are p{t) {+) and 1 — p{t) (— ). 
In the mean field limit, a random group of fc agents is 
selected, with j agents with opinion + and k — j with 
opinion — . The opinion dynamics of the GPM enforces 
consensus among the agents of the group, which adopt 
opinion -I- with a suitably defined probability m^j and 
opinion — with probability 1 — m^ j. The probability 
p(t + 1) to find an agent sharing opinion + after the 
update is 



p(i+l) = ^m,,,p(t)^"[l-p(i)]'=- 



fc! 



(26) 



The size fc of the random group along with the local prob- 
abilities {mfej} completely define the dynamics of the 
GPM. A phase diagram can be derived as a function of 
the local probabilities. Only two different phases arc ob- 
tained, corresponding to consensus and coexistence of the 
two opinions in equal proportions. The phase transition 
occurs at those values of the {rukj} for which magnetiza- 
tion is on average conserved: here the model has a voter 
dynamics. With suitable choices of the set {ruk j} the 
GPM reproduces the MP behavior of all known models 
with binary opinions: voter, majority rule, Sznajd, the 
majority- minority model, etc.. 

We now briefly review the modifications of the Sznajd 
model. The dynamics has b een studied on many different 
topol ogies: regular lattices (Chang, '2001';' Stauffer et all 
\200(t ) , comple te graphs (Slan ina and Layicka. ,20031 , ran- 
dom graphs ( Rodrigues and da F. CostSTioOSh, s mall- 



scale- free networks (Bernardes et all 



world networks (lElgazzar . 2003[~He et al. 



2002 



200i) and 
Bonnckoh 



2003t iRodrigues and da F. Costal 120051: ISousal 12005 



Sousa and Sanchez . 2006( ). The S znajd model 011 scale 
free networks was recently studied ( Gonzalez et a/.l . [2006l ) 
within a real space renormalization framework. On any 
graph, if only Sznajd's ferromagnetic rule (|2ip holds, 
the system undergoes a sharp dynamic phase transi- 
tion from a state with all spins down to a state will all 
spins up. If the graph is not fixed, but in evolution, 
like a growing network, the transi tion becomes a smoot h 
crossover between the two phases ( Gonzalez et all . \2004) . 
The ph ase transition holds as well if one introduces di- 
lution ( Moreira et Ml, 120011). if the number of opin ion 
states is larger than two (jSlanina and Lavickal [200! ■ if 
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the influence of the active pa i r of ag ents extends beyond 
their neighborhood (ISchulzel . [2003bh . so it is a very ro- 
bust feature of the Szna.i d model, although it disappears 
when one includes noise (Stauff er et al. . 2000l) or ant ifer- 
romagnetic rules ( Chan g. 2001: iSznaid-Weronl . 12004 ). 

If the random sequential updating so far adopted is 
replaced by synchronous updating, i.e., if at each it- 
eration all agents of the configurations are paired off 
and act simultaneously on their neighbors, it may hap- 
pen that an agent is induced to choose opposite opin- 
ions by different neighboring pairs. In this case the 
agent is "frustrated" and mai ntains it s opiii i on. Such 
frustration hinders consensus (jStauffed . [200I ITu et all . 
I2005D . due to the emergence of stable clusters where 
both opinions coex ist. This problem can be limi ted if 
noise is introduced ( Sabatelli and Richmondl . [20041 ). or if 
agents have memory, so that, in case of conflicting ad- 
vice, they follow the most frequent opin ion they had in 
the past ( Sabatelh and Richmond! . [20031 ). 

When the possible opinion states are q > 2, one can 
introduce bounded confidence, i.e., the realistic princi- 
ple that only people with similar opinions can have an 
influence on each other. If we assume that two opin- 
ions are similar if their values differ by at most one unit, 
and that a pair of agents with the same opinion can con- 
vince only neighbors of similar opinions, the Sznajd dy- 
namics always leads to complete consensus for g = 3, 
whereas for g > 3 it is very likely that at least two opin - 
ions survive in the final stationary state ( Stauffeil . l2002bt ). 
Bounded confidence allows for an extension of the Sznajd 
model to real- valued opinions ( Fortunatol . l2005b[ ) . Other 
studies focused on the dynamics of clusters of agents 
with regular opi nion patterns, ferromagnetic and /or an- 
tiferromagnetic |Schneider and Hirtreitei , 2005bl ) j dam - 
age spreading (jKhetsch . 2005 : Roehner et all \2004) . 
the comb ination of Szna.id with other convincing 
stra tegies (ISousa and Sanched . 20061). contraria n behav- 
ior (|de la Lama et all l2005l : Iwio et all l2006t) . the ef- 
fect on the dynamics of agents biased t owards the 
global majority and/or minority op inion (|Schneideij . 
120041 : [Schneider and Hirtreiteil . [2005af ). Recently exten- 
sions of Szna .id B to higher dimensions have bee n consid- 
ered as well (iKondrat and Sznai d-WeronLl2008h . 

The Sznajd model has found applications in dif- 
ferent areas. In politics, it has b een used to de- 
scrib e voting behavi o r in elections ( Bernardes et all 
l2002t iGonzalez et'ol] . \2004 : we shall discuss this is- 
sue in Sec. IIII.HI Moreover, it was applied to study 
the interaction of economic and personal attitudes of 
individuals, which evolve a ccording to different rules 
but in a coupled manner ( Sznaid-Weron and Sznaidl . 
I2OO5I ). Sznajd dynamics has also been adopted to 
model the co mpetition of different products in an 
open market ( Sznaid-Weron and Weronl l2003l ). The 
effects of aging, diffusion and a mult i -layered society 
have been considered as well ( Schulzd . l2003al : ISchulzd . 
|2004( ). Sznajd dynamics has been adapted in a model 
that describes the spread of opinions among a group 



of traders ( Sznaid-Weron and Weronl . |2002| ). Finally, 
Sznajd-like rules have bee n employed to generate a new 
class of complex networks ( da Fontoura Costal . l2005l ). 



F. Bounded confidence models 

1. Continuous opinions 

In the models we have so far investigated opinion is a 
discrete variable, which represents a reasonable descrip- 
tion in several instances. However, there are cases in 
which the position of an individual can vary smoothly 
from one extreme to the other of the range of possible 
choices. As an example, one could think of the politi- 
cal orientation of an individual, that is not restricted to 
the choices of extreme Right/Left, but it includes all the 
options in between, which may be indicated by the geo- 
metric position of the seat of a deputy in the Parliament. 

Continuous opinions invalidate some of the concepts 
adopted in models with discrete choices, like the con- 
cepts of majority of an opinion and equality of opinions, 
so they require a different framework. Indeed, continu- 
ous opinion dynamics has historically followed an alter- 
native path. The first studies were carried out by applied 
mathematicians and aimed at identifying the conditions 
under which a panel of experts would reach a common de- 
cision (IChatteriee and Senetal . llQTTt ICohen et all \l98d 
IStond . ll96lh " 

The initial state is usually a population of N agents 
with randomly assigned opinions, represented by real 
numbers within some interval. In contrast to discrete 
opinion dynamics, here all agents usually start with dif- 
ferent opinions, and the possible scenarios are more com- 
plex, with opinion clusters emerging in the final station- 
ary state. The opinion clusters could be one (consensus), 
two (polarization) or more (fragmentation). In princi- 
ple, each agent can interact with every other agent, no 
matter what their opinions are. In practice, there is a 
real discussion only if the opinions of the people involved 
are sufficiently close to each other. This realistic aspect 
of human communications is called bounded confidence 
(BC); in the literature it is expressed by introducing a 
real number e, the uncertainty or tolerance, such that an 
agent, with opinion x, only interacts with those of its 
peers whose opinion lies in the interval ]x — e, x -\- e[. 

In this section we discu ss the most popular B C models, 
i.e., the Deffuant model (iDeffuant et all\200(ji] and that 
of Hegselmann-Krause (jHegselmann andKrausfJ, '2002^. 
BC models have been recently reviewed in (|Lorenz . .2007.) . 



2. Deffuant model 

Let us consider a population of N agents, represented 
by the nodes of a graph, where agents may discuss with 
each other if the corresponding nodes are connected. 
Each agent i is initially given an opinion Xi, randomly 
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chosen in the interval [0, 1]. The dynamics is based on 
random binary encounters, i.e., at each time step, a ran- 
domly selected agent discusses with one of its neighbors 
on the social graph, also chosen at random. Let i and j 
be the pair of interacting agents at time with opinions 
Xi{t) and Xj{t), respectively. DefFuant dynamics is sum- 
marized as follows: if the difference of the opinions Xi{t) 
and Xj{t) exceeds the threshold e, nothing happens; if, 
instead, \xi{t) — Xj{t)\ < e, then: 

x,{t + 1) = x,{t) + ti[xj{t) - x,{t)l (27) 

xj{t + 1) = Xj{t) + ii[x,{t) - Xj{t)]. (28) 

The parameter /i is the so-called convergence param- 
eter, and its value lies in the interval [0,1/2]. Deffuant 
model is based on a compromise strategy: after a con- 
structive debate, the positions of the interacting agents 
get closer to each other, by the relative amount fx. If 
yU = 1/2, the two agents will converge to the average of 
their opinions before the discussion. For any value of 
e and /i, the average opinion of the agents' pair is the 
same before and after the interaction, so the global av- 
erage opinion (1/2) of the population is an invariant of 
DefFuant dynamics. 

The evolution is due to the instability of the initial 
uniform configuration near the boundary of the opinion 
space. Such instability propagates towards the middle of 
the opinion space, giving rise to patches with an increas- 
ing density of agents, that will become the final opinion 
clusters. Once each cluster is sufficiently far from the 
others, so that the difference of opinions for agents in 
distinct clusters exceeds the threshold, only agents inside 
the same cluster may interact, and the dynamics leads to 
the convergence of the opinions of all agents in the cluster 
to the same value. Therefore, the final opinion configura- 
tion is a succession of Dirac's delta functions. In general, 
the number and size of the clusters depend on the thresh- 
old e, whereas the parameter /i affects the convergence 
time of the dynamics. However, when is small, the fina l 
cluster configuration als o depends on ^ ( Laguna et al\ . 
120041 : iPorfiri aLl . l2007D . 

On complete graphs, regular lattices, random graphs 
and scale-free networks, for e > — 1/2, all agents 
share t he same opinion 1/2, so there is complet e con- 
sensus (|Fortunatol . |2004 iLorenz and Urbid , |2007t ). This 



may be a general property of Deffuant model, indepen- 
dently of the underlying social graph. If e is small, more 
clusters emerge (Fig. IH]). 

Monte Carlo simulations reveal that the number Uc of 
clusters in the final configuration can be approximated 
by the expression l/(2e). This can be understood if we 
consider that, at stationarity, agents belonging to dif- 
ferent opinion clusters cannot interact with each other, 
which means that the opinion of each cluster must differ 
by at least e from the opinions of its neighboring clusters. 
In this way, within an interval of length 2e centered at 
a cluster, there cannot be other clusters, and the ratio 
l/(2e) is a fair estimate for ric- 
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FIG. 6 Deffuant model. Opinion profile of a population of 500 
agents during its time evolution, for e = 0.25. The population 
is fully mixed, i.e., everyone may interact with everybody else. 
The dynamics leads to a polarization of the population in two 
factions. 



Most results on DefFuant dynamics are derived through 
numerical simulations, as the model is not analytically 
solvable. However, in the special case of a fully mixed 
population, where everybody interacts with everybody 
else, it is possible to write the general rate eq uation 
governing the opinion dynamics (Bcn-Nai m et aZ.1 . 12003). 
For this purpose, one neglects individual agents and fo- 
cuses on the evolution of the opinion population P{x^ t), 
where P(cc, t)dx is the probability that an agent has opin- 
ion in the interval [x, x-^dx\. The interaction threshold is 
e = 1, but the opinion range is the interval [—A, A]; this 
choice is equivalent to the usual setting of the Deffuant 
model, if e = 1/2A. For simplicity, ^ — 1/2. The rate 
equation then reads: 



d_ 

dt 



Pix,t) 



/ / dxidx2P{xi,t)P{x2,t) 

J \xi—X2\<l J 

5(x~'^^^)-5{x~x,)\. (29) 



Eq. conserves the norm Mq = P{x, t)dx and the 
average opinion. The question is to find the asymptotic 
state Poo{x) — P{x, t — !■ oo), starting from the flat initial 
distribution P{x,t = 0) = 1, for x e [-A, A]. If A < 
1/2, all agents interact and Eq. ([29|) is integrable. In this 
case, it is possible to show that all agents approach the 
central opinion and Poo (a:) = Mo(5(x). 

If A > 1/2, the equation is no longer analytically solv- 
able. The asymptotic distribution is a linear combination 
of delta functions, i.e.. 



Poo(x) 



1=1 



miS{x — Xi). 



(30) 
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The cluster masses nii must obey the conditions J^i "^i = 
Mq and J2i TTT'i^i = 0; the latter comes from the con- 
servation of the average opinion. Numerical solutions of 
Eq. reveal that there are only three types of clusters: 
major (mass > 1), minor (mass < 10~^) and a central 
cluster located at a; = 0. These clusters are generated 
by a periodic sequence of bifurcations, consisting in the 
nucleation and annihilation of clusters. 

On any graph topology, Deffuant dynamics is always 
characterized by a stationary state with one or more opin- 
ion clusters. However, as the interaction range of an 
agent is restricted to its topological neighborhood, more 
opinion clusters emerge for low values of the uncertainty. 
Opinion homogenization involves only agents in the same 
cluster: in this way, if two clusters are geometrically sepa- 
rated, there will be no communication between the corre- 
sponding agents and the final opinions will be in general 
different in each cluster, even if their opinions are com- 
patible, so that they would converge to the same opinion 
on a complete graph. The result is an increased fragmen- 
tation of the agents' population. On scale-free networks, 
the number of surviving opinions in the stationary state 
is proportional to the number of agents of the popula - 
tion, for fixed e ( Staufi^er and Mever-Ortmannd . 12004 ). 
In particular, nodes with few connections have a sizeable 
probability to be excl uded from the dy namics and to keep 
their opinion forever (IWeisbuchl . [200l . The result holds 
for both static and evolving networks ( Sousal . [2004 ). 

Deffuant model can be defined as well if opinion s 
are not continuous but discretized ( Stauffer et al ]. |2004) . 
Here the opinion s of any agent can take one of Q values, 
s = 1, 2, Q. Opinions Si and Sj ^ Si are compatible if 
|si ~ Sjl ^ L, where L is an integer. The rules are the 



same as in Eqs. (|27)) and ((28|) . still with a real- valued 
convergence parameter /i, but the shift of the opinions 
is rounded to the nearest integer. In the limit L — s- oo 
and Q oo, with the ratio e = L/Q kept constant, one 
recovers the results of the original model with continu- 
ous opinions. If L = 1, on a complete graph consensus is 
the only stationary state if Q = 2, i.e., for binary opin- 
ions^. Instead, for Q > 2, complete consensus is never 
attained, but there are at least two opinion clusters in the 
final configuration. On scale-free networks the number 
of surviving opinions in the stationary state is a scaling 
function of Q and the population size N. 

Simple modifications of Deffuant model yield 
ric h dynam i cs. I f agents have in di vidua l values of 
e (|Lorenzl . l2008bl : IWeisbuch efaj] . |2002| ). the dy- 
namics is dominated by the ag ents with large uncer- 
tainties. In several papers (lAldashev and Carletti 
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2006 



Amblard and Deffuantl . |2004 iDeffuant 



iDeflFuant aU |2004 iDeffuant e< aU 12002 : 



We remark that, for L = 1, it is impossible for the two interact- 
ing opinions to shift towards each other, as only integer opinion 
values are allowed; so, as a result of the discussion, one agent 
takes the opinion of the other. 



IWeisbuch et al. I, iool), the uncertainties are also 
affected by the dynamics. In addition, they are also cou- 
pled to the opinions, based on the principle that a small 
uncertainty also implies more confidence and a higher 
probability to affect other opinions. These models are 
able to explain how extremal positions, initially shared 
by a minority of people, may eventually dominate in 
society. In (Bcn-Naim, 2005) a model in which Deffuant 
compromise strategy is combined with spontaneous 
changes of the agents' opinions, has been studied. The 
latter phenomenon is described as a diffusion process in 
the opinion space, which affects cluster formation and 
evolution, with large clusters steadily overtaking small 
ones. Other extensions of Deffuant dynamics include 
the introduction of an external periodic perturbation 
affecting all agents at onc e , to simulate the effect of 
propaganda ( Carletti et all . l2006t) . and the study of a 
more complex opinion dynamics where the interaction of 
pairs of agents depends not only on the compatibility of 
their opinions , but also on the coevolving mutual affinit y 
of the agents (|Bagnon et all . l2007t ICarletti et aLl . l2007t) . 
This coupling provides a natural and endogenous way 
of determining the number of opinion clusters and their 
positions. 



3. Hegselmann-Krause model 



Th e model proposed in ("Hegse lmann and Krausel 
l2002f ) (HK) is quite similar to that of Deffuant. Opinions 
take real values in an interval, say [0, 1], and an agent i, 
with opinion Xi , interacts with neighboring agents whose 
opinions lie in the range [xi — e,Xi + e], where e is the 
uncertainty. The difference is given by the update rule: 
agent i does not interact with one of its compatible neigh- 
bors, like in Deffuant, but with all its compatible neigh- 
bors at once. Deffuant's prescription is suitable to de- 
scribe the opinion dynamics of large populations, where 
people meet in small groups, like pairs. In contrast, HK 
rule is intended to describe formal meetings, where there 
is an effective interaction involving many people at the 
same time. 

On a generic graph, HK update rule for the opinion of 
agent i at time t reads: 



Xi{t - 



^ T,jj\x,(t)-x,{t)\<e "-V^jii) ^^^^ 



{t)-'Xj{t)\<e 



where is the adjacency matrix of the graph. So, agent 
i takes the average opinion of its compatible neighbors. 
The model is fully determined by the uncertainty e, un- 
like Deffuant dynamics, for which one needs to specify 
as well the convergence parameter fi. The need to cal- 
culate opinion averages of groups of agents that may 
be rather large makes computer simulations of the HK 
model rather lengthy as compared to Deffuant's. This 
may explain why the HK model has not been studied by 
many authors. 
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The dynamics develops just like in DefFuant, and leads 
to the same pattern of stationary states, with the num- 
ber of final opinion clusters decreasing if e increases. In 
particular, for e above some threshold (-c, there can only 
be one cluster. On a complete graph, the final configura- 
tions are symmetric with respect to the central opinion 
1/2, because the average op inion of the system is con- 
served by the dynamics (jFortunato et ai . 2005), as in 
DefFuant. The time to reach the stationary state diverges 
in correspondence to the bifurcation thresholds of opin- 
ion clusters, due to the prese nce of isolated agents ly ing 
between consecutive clusters (jFortunato et aZ.I . |2005( ). 

The threshold for complete consensus Ec can only take 
one of two values, depending on the behavior of the av- 
erage degree (fc) of the underlyi ng social graph w hen the 
number of nodes N grows large (lFortunatol . l2005aD . If (/c) 
is constant in the limit of large N, as for example in lat- 
tices, ec = £1 = 1/2. Instead, if (fc) oo when N ^ oo, 
as for example in complete graphs, ec = £2 ~ 0.2. We 
have seen instead that, for Deffuant, ec = 1/2 on any 
graph. 

Th e extensio n of HK to discretized opin- 
ions ( Fortunatol . |2004[ ) is essentially a voter model 
with bounded confidence: an agent picks at random the 
opinion of a compatible neighbor. For three opinion 
values and uncertainty one, the model reduc es to the 
constrained voter model ( Vazquez et aLl . [20031 ) . 

Other developments include: the use of al- 
ternative recipes to average t he o pinions 



Eq. (|3T|) (jHegselmann and Krausd. I2005f ) ; an anal- 
ysis of damage spreading ( Fortunatol . l2005h : the 
introduction of a general framework where the size of 
the groups o f interacting agent s varie s from 2 (Deffuant) 
to N (HK) (jUrbig and Lorend . l200i ): the reformulation 
of Deff uant and HK dynamics as interactive Markov 
chains (|Lorenzl . [2005bl : II orenj . |2006[ ) ; a nalytical results 
on the stability of EC opinion dynamics ( Loren j . l2005al ) 
and their ability to preserv e the relative ordering of the 
opinions ([Hendrick 2 l2007D . 



G. Other models 

The opinion dynamics models described so far are 
based on elementary mechanisms, which explain their 
success and the many investigations they have stimu- 
lated. Such models, however, do not exhaust the wide 
field of opinion dynamics. The last years witnessed a 
real explosion of new models, based on similar concepts 
as the classical models or on entirely new principles. Here 
we briefly survey these alternative models. 

The basic models we have seen are essentially deter- 
ministic, i.e., the final state of the system after an in- 
teraction is always well defined. Randomness can be in- 
troduced, in the form of a social temperature or pure 
noise, but it is not a fundamental feature. Most models 
of last generation, instead, focus on the importance of 
randomness in the process of opinion formation. Ran- 



domness is a necessary ingredient of social interactions: 
both our individual attitudes and the social infiuence of 
our peers may vary in a non-predictable way. Besides, 
the infiuence of external factors like mass media, propa- 
ganda, etc., is also hardly predictable. In this respect, 
opinion dynamics is a stochastic process. 

An interesting variant of Ising dyna r nics h as been re- 
cently explored in (Jiang ef^aZT l2007l . l2008l ). Here an 
agent surrounded by a majority of neighbors with equal 
opinion will fiip its opinion with the usual Metropolis 
probability exp{—AE/T), where AE is the increase of 
the Ising energy due to the flip and T the (social) temper- 
ature. When the majority of the neighbors disagrees with 
the agent's opinion, the latter flips with a probability q, 
which accounts for the possibility that agents keep their 
opinion in spite of social pressure ( "inflexible agents" ) . 
For q = 1 one recovers the usual Ising dynamics. For 
q < 1 instead, the model obeys a non-equilibrium dynam- 
ics, since detailed balance is violated. Simulations of the 
model on the lattice reveal a non-trivial phase diagram: 
the ferromagnetic transition is continuous for q larger 
than a critical value qc, and discontinuous for q < qc- 
The critical exponents of the continuous transitions dif- 
fer both from the mean fleld and the Ising exponents. 
On small-world lattices with a large density of shortcuts 
ther e is no mag i ietiza tion transition if q is sufficiently 
low (| Jiang et aU |2008[ ). 

In ( Bartolozzi et all l2005l ) a model with binary opin- 
ions, evolving according to a heat bath dynamics, is pro- 
posed. The opinion field acting on a spin is given by a 
linear combination with random weights of a term pro- 
portional to the average opinion of its nearest neighbors 
on the social network, with a term proportional to the av- 
erage opinion of the whole network. When the stochastic 
noise exceeds a threshold, the time evolution of the aver- 
age opinion of the system is characterized by large inter- 
mittent fiuctuations; a comparison with the time series 
of the Dow- Jones index at New York's Stock Exchange 
reveals striking similarities. 

In a recent model ( Kuperman and Zanettd . l2002f ). 
opinions are affected by three processes: social imitation, 
occurring via majority rule, fashion, expressed by an ex- 
ternal modulation acting on all agents, and individual 
uncertai nty, expressed by random noise. Stochastic res- 
onance ( Gammaitoni et all . 1 19981 ) was observed: a suit- 
able amount of noise leads to a strong amplification of the 
response of the system to the external modulation. The 
phenomenon occurs as well if one v aries the size of the 
system for a fixed amount of noise ( Tessone and Tor"ail . 
[2005) : here the best response to the external solicitation 
is achieved for an optimal population size (system size 
stochastic resonance) . 

Kinetic models of opinion dynamics were proposed 
m ([Toscanil . |2006| ). Interactions are binary, and the 
opinions of the interacting pair of agents vary according 
to a compromise strategy a la Deffuant, combined with 
the possi bility of opinion diffusion, following the origi- 
nal idea (|Ben-Naiml . 120051 ) discussed in Sec. IIII.F.2I The 
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importance of diffusion in the process is expressed by a 
random weight. The dynamics can be easily reformu- 
lated in terms of Fokker-Planck equations, from which it 
is possible to deduce the asymptotic opinion configura- 
tions of the model. Fokker-Planck equations have also 
been employed to study a d ynamics similar to th at of 
the constrained voter model ( Vazquez et all 12003 ^ , but 
in the presence of a soc ial temperature, inducin g sponta- 
neous opinion changes ( de la Lama et aLl . l2006l ). 



In (jMartinsl . l2007l HqoI) binary and continuous opin- 
ions are combined. The model is based on the simple 
idea that declared opinions, or actions, are only a pro- 
jection of the actual opinion of the agents. Two persons 
can express the same preference but their individual cer- 
titudes towards that preference may be different. So one 
has to distinguish between the internal opinion, which is 
expressed by a probability, and the external opinion, or 
action, which is a binary variable. The internal opinion 
is a measure of how extreme a position is. Agents vary 
both their actions and their internal opinions based on 
the observation of the actions of their peers. In this way 
one can monitor how the convictions of individuals are 
coupled to their actions. Clusters of agents with the same 
external opinions display a characteristic pattern, where 
the agents are very convinced of their choices (extrem- 
ists) if they are well inside the cluster, whereas they are 
more open-minded if they lie at the boundary. 

Synchronization has also been used to explain 
consen sus formation. A variant of the Kuramoto 
model (|Kuramotolll975t ). where the phases of the oscilla- 
tors are replaced by unbounded real numbers, represent- 
ing the opinions, displays a phase transition from an in- 
coherent phase (anarchy) , to a s y nchronized pha s e (con - 
sen sus) (iPluchino et all. 120061: iPluchino et'dl . 120051 ). 
In ( Di Mare and Latoral [2007 ) it was shown that sev- 
eral opinion dynamics models can be reformulated in the 
context of strategic game theory. 

Some mo dels focus on s pecific aspects of opinion dy- 
namics. In ( Indekeul . 12004 ) it has been pointed out that 
the influence of network hubs in opinion dynamics is 
overestimated, because it is unlikely that a hub-agent 
devotes much time to all its social contacts. If each 
agent puts the same time in its social relationships, this 
time will be distributed among all its social contacts; 
so the effective strength of the interaction between two 
neighboring agents will be the smaller, the larger the de- 
grees of the agents. If the spin-spin couplings are renor- 
malized according to this principle, the Ising model on 
scale- free networks always has a ferromagnetic threshold, 
whereas it is known that, with uniform couplings, net- 
works with in finite degree variance are magnetized at an y 
temperature ( Aleksieiuk et a/.l . l2002HLeone et al. 1. 120021) . 
The issue of how opinion dynamics is influenced by 
the hierarchical struct ure in societies/organizations has 
also been inv e stigat ed ( Grabowski and Kosinskil . l2006bl : 



iLaguna et al. . 2005t). Other authors investigated fash- 



ion ( Nakavama and Nakamural 12004 ), the interplay be- 
tween opinions and personal taste ( Bagnoli et all . 12004 ) 



and t he effect of opinion surveys on the outcome of elec- 
tions ijAlves et a/.l . |2002| ). 

It is worth mentioning how the close formal similari- 
ties between the fields of opinion and language dynamics 
leads to the idea that models proposed in the framework 
of language dynamics could suitably apply also in mod- 
eling opinion formation. One exampl e is repre sented by 
a variant of the naming game (Baronchelli et al, . ,2007.) . 
as defined in Sec. IVl 



H. Empirical data 

One of the main contributions of the physical approach 
to opinion dynamics should be to focus on the quantita- 
tive aspects of the phenomenon of consensus formation, 
besides addressing the mere qualitative question of when 
and how people agree/disagree. What is needed is then 
a quantitative phenomenology of opinion dynamics, to 
define the phenomenon in a more objective way, posing 
severe constraints on models. Sociological investigations 
have been so far strongly limited by the impossibility of 
studying processes involving large groups of individuals. 
However, the current availability of large datasets and of 
computers able to handle them makes for the first time 
such empirical analysis possible. 

Elections are among the largest scale phenomena in- 
volving people and their opinions. The number of vot- 
ers is of the order of millions for most countries, and 
it can easily reach hundreds of millions in countries like 
Brazil, India and the USA. A great deal of data is nowa- 
days publicly available in electronic form. The first em- 
pirical investigation s carried out by physicis ts concerned 
Brazilian elections ( Costa Filho et all . 19991 ). The study 
focused on the distribution of the fraction v of votes re- 
ceived by a candidate. Datasets referring to the federal 
elections in 1998 revealed the existence of a characteris- 
tic pattern for the histogram -P(i^), with a central portion 
following the hyperbolic decay l/v, and an exponential 
cutoff for large values of ly. Interestingly, datasets cor- 
responding to candidates to the office of state deputy in 
several Brazilian states revealed an analogous pattern. 
A successive analysis on data referring to state and fed- 
eral elections in 2002 confirmed the results for the elec- 
tions in 1998, in spite of a change in the p olitical rules 
that constrained alliances between parties (iFilho et all . 
|2003[ ). Indian data displayed a similar pattern for P(j^) 
across different states, altho ugh discrepancies were also 
found (^Gonzal ez et 'ail. l2004) . Data on Indonesian elec- 
tions are consistent with a power law decay of Piv), 
with e xponent close to one, but are too noisy to be re- 
liable (|Situngkiij . [200l . Claims that Mexican elections 
also obey a similar pattern are not clearly supported by 
the data ( Morales-Matamoros et "aZl. l2006D . 

The peculiar pattern of P{i') was interpreted as the 
result of a multiplicative process, which yields a log- 
norma l distribution for v , due t o the Central Limit The- 
orem (ICosta Filho et all . Il999tl . The 1/ V behavior can 



23 



indeed be reproduced by a log-normal function, in the 
limit where the latter has a large variance. A micro- 
scopic m odel based on Szna.id op inion dynamics was pro- 
posed in teer nardes et aLl . l2002l ). Here, the graph of per- 
sonal contacts between voters is a scale-free network d la 
Barabasi- Albert; candidates are initially the only nodes 
of the network in a definite opinion state, a suitably mod- 
ified Sznajd dynamics spreads the candidates' opinions 
to all nodes of the network. The model reproduces the 
empirical curve P(i^) derived from Brazilian elections. 
The same mechanism yields, on different social gr aphs, 
like pseudo-fractal networks ( Gonzalez et all \2004) and 
a m odified Bara basi- Albert network with high cluster- 
ing (jSousal . 120051 ^ ■ a good agreement with empirical data. 
The weakness of this model, however, is that a non-trivial 
distribution is only a transient in the evolution of the sys- 
tem. For long times the population will always converge 
to the only stable state of Sznajd dynamics, where every 
voter picks the same candidate, and the corresponding 
distribution is a J-function. All studies stopped the mod- 
ified Sznajd dynamics after a certain, carefully chosen, 
time. A recent model based on simple opinion spread- 
ing yields a distribution similar to the Brazilian curve, if 
the underlying social graph is an Erdos-Renyi network, 
whereas on scale-free networks the same dynamics fails 
to reproduce the data ( Travieso and da Fontoura Costal . 
[2006) . 

The power law decay in the central region of -P(t^), 
observed in data sets relative to different countries and 
years, could suggest that this pattern is a universal fea- 
ture of the distribution. But this is unlikely because can- 
didates' scores strongly depend on the performance of 
their parties, which is determined by a much more com- 
plex dynamics. Ind eed, munic i pal el ection data display 
a different pattern ( Lyra e^"all . |200^ . Instead, the per- 
formances of candidates of the same party can be objec- 
tively compared. This can be done in proportional elec - 
tions with open lists ( Fortunato and Castellancl l2007t ). 
In this case the country is divided into constituencies, 
and each party presents a list of candidates in each con- 
stituency. There are three relevant variables: the number 
of votes V received by a candidate, the number Q of can- 
didates presented by the party in the corresponding list 
and the total number N of votes received by the party 
list. Therefore, the distribution of the number of votes 
received by a candidate should be a function of three vari- 
ables, P{v,Q,N). It turns out instead that P{v,Q,N) 
is a scaling function of the single variable vQ/N, with a 
log-normal shape, and, remarkably, this function is the 
same in different countries and years (Fig. [7|) . This find- 
ing justifies a simple microscopic description of voting be- 
havior, using the tools and methods of statistical physics. 
A model based on word-of-mouth spread i ng, si milar to 
that of ( Travieso and da Fontoura Costal . [20061 ). is able 
to reproduce the data. 

Other studies disclose a correlation between the scores 
of a part y and the number of its members in German 
elections ( Schneider and Hirtreiteil . l2005cl ) and a polar- 
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FIG. 7 Distribution of electoral performance for candidates in 
proportional elections held in Italy, Poland and Finland. The 
remarkable overlap shows that th e curve is a universal fea- 
ture of the voting process. From ([Fortunato and Castellanol . 
120071 ). 
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FIG. 8 Relation between the maximal speed of change and 
the duration of the change for birth rates and the num- 
ber of mobile phones in several European countries. The 
linear regression of the data points in double logarithmic 
scale is consistent with the universal behavior predicted 
by the Rando m Field Is ing Model at zero temperature. 
From (,Michard and Bouchaud . .200511 . 



ization of the distribution of votes around two main can- 
didat es in Brazilian elections for mayor ( Araripe et all 
I2OO6I ) . An empirical investigation of the relation between 



party size and temporal correlations in electoral results 
is described in ( Andresen et aLl . 12007^. 

In ( Michard and Bouchaudl [2005i) it was suggested 
that extreme events like booms of products/fashions, fi- 
nancial crashes, crowd panic, etc., are determined by 
a combination of effects, including the personal atti- 
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tude of the agents, the public information, which affects 
all agents, and social pressure, represented by the mu- 
tual interaction between the agents. This can be for- 
mally described within the framework of the Random 
Field Ising Model at zero temperature, which success- 
fully describes hysteresis in random magnets and other 
physical ph enomena, like the occurrence of crackling 
noise dSethna et all . l200ll ). Here, opinions are binary, 
attitudes are real- valued numbers in ] — cx), -l-oo[, corre- 
sponding to the random fields, the public information is 
a global field F{t), slowly increasing with the time t, and 
the interaction term is the sum of Ising-like couplings 
between pairs of agents. The order parameter O of the 
system is the average opinion of the population. By in- 
creasing the field F, O displays a sharp variation, due to 
large groups of agents that simultaneously switch opin- 
ion. The evolution of the speed of change dO/dF as a 
function of F follows a universal bell-shaped curve in the 
transition region, with a characteristic relation between 
the height h of the peak and its width w: h ^ u;"^/^. 
This relation was indeed observed in empirical data on 
extreme events, such as the dramatic drop of birth rates 
in different European countries in the last decades, the 
rapid diffusion of mobile phones in Europe in the late 
'90s, and the decrease of the clapping intensity at the 
end of applauses (Fig. [5]). 

For the future, more data are needed. Several phe- 
nomena of consensus formation could be empirically an- 
alyzed, for instance spreading of fads and innovations, 
sales dynamics, etc.. 



from both social scientists and physicists. 

The origin of its success among social scientists is in 
the inclusion of two mechanisms that are believed to be 
fundamental in the understanding of the dynamics of cul- 
tural assimilation (and diversity): social influence and 
homophily. The first is the tendency of individuals to be- 
come more similar when they interact. The second is the 
tendency of likes to attract each other, so that they inter- 
act more frequently. These two ingredients were generally 
expected by social scientists to generate a self-reinforcing 
dynamics leading to a global convergence to a single cul- 
ture. It turns out instead that the model predicts in some 
cases the persistence of diversity. 

From the point of view of statistical physicists, the Ax- 
elrod model is a simple and natural "vectorial" general- 
ization of models of opinion dynamics that gives rise to a 
very rich and nontrivial phenomenology, with some gen- 
uinely novel behavior. The model is defined as follows. 
Individuals are located on the nodes of a network (or on 
the sites of a regular lattice) and are endowed with F 
integer variables (cti, . . . ^up) that can assume q values, 
CT/ = 0, 1,...,(7 — 1. The variables are called cultural fea- 
tures and q is the number of the possible traits allowed 
per feature. They are supposed to model the different 
"beliefs, attitudes and behavior" of individuals. In an 
elementary dynamic step, an individual i and one of his 
neighbors j are selected and the overlap between them 



3 = jll ^-f 



(32) 



IV. CULTURAL DYNAMICS 

In the previous section we have reviewed the active 
field of opinion dynamics. In parallel, there has been 
in recent years a growing interest for the related field of 
cultural dynamics. The border between the two fields 
is not sharp and the distinction is not clear-cut. The 
general attitude is to consider opinion as a scalar variable, 
while the more faceted culture of an individual is modeled 
as a vector of variables, whose dynamics is inextricably 
coupled. This definition is largely arbitrary, but we will 
adopt it in the review. 

The typical questions asked with respect to cultural 
influence are similar to those related to the dynamics 
of opinions: what are the microscopic mechanisms that 
drive the formation of cultural domains? What is the 
ultimate fate of diversity? Is it bound to persist or all 
differences eventually disappear in the long run? What 
is the role of the social network structure? 



A. Axelrod model 

A prominent role in the investigation of cultural dy- 
namics has_been played by a model introduced by Axel- 
rod in (|Axelrodl . fl997h . that has attracted a lot of interest 



is computed, where 5ij is Kronecker's delta. With prob- 
ability cuij the interaction takes place: one of the fea- 
tures for which traits are different (cr/(i) ^ cr/(j)) is se- 
lected and the trait of the neighbor is set equal to <Jf{i)- 
Otherwise nothing happens. It is immediately clear that 
the dynamics tends to make interacting individuals more 
similar, but the interaction is more likely for neighbors 
already sharing many traits (homophily) and it becomes 
impossible when no trait is the same. There are two sta- 
ble configurations for a pair of neighbors: when they are 
exactly equal, so that they belong to the same cultural 
region or when they are completely different, i.e., they 
sit at the border between cultural regions. 

Starting from a disordered initial condition (for exam- 
ple with uniform random distribution of the traits) the 
evolution on any finite system leads unavoidably to one of 
the many absorbing states, which belong to two classes: 
the ordered states, in which all individuals have the 
same set of variables, or the other, more numerous, frozen 
states with coexistence of different cultural regions. 

It turns out that which of the two classes is reached 
depends o n the number of possibl e traits q in the initial 
condition ( Castellano et all . |2000| ) . For small q individ- 



uals share many traits with their neighbors, interactions 
are possible and quickly full consensus is achieved. For 
large q instead, very few individuals share traits. Few 
interactions occur, leading to the formation of small cul- 
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FIG. 9 Axelrod model. Behavior of the order parame- 
ter {Smax) / L'^ VS. q for three different system sizes and 
F = 1 0. In the inset the same quantity is reported for F — 2. 
From (|Castellano et all l200d ) . 



tural domains that are not able to grow: a disordered 
frozen state. On regular lattices, the two regimes are 
separated by a phase transition at a critical value qc, de- 
pending on F (Fig. [9]). 

Several order parameters can be defined to character- 
ize the transition. One of them is the average fraction 
{Smax) /N of the system occupied by the largest cultural 
region. N is the number of individuals in the system. 
In the ordered phase this fraction is finite (in the limit 
N ^ oo), while in the disordered phase cultural domains 
are of finite size, so that {Smax )IN ^ 1/N. Another 



(dis)order parameter often used (| Gonzalez- Avella et al\ . 
I2OO5) is g = (Ng) /N , where Ng is the number of different 
domains in the final state. In the ordered phase g — > 0, 
while it is finite in the disordered phase. 

In two dimensions the nature of the transition depends 
on the value of F. For F — 2 there is a continuous change 
in the order parameter at qc, while for F > 2 the tran- 
sition is discontinuous (Fig. Correspondingly the 
distribution of the size s of cultural domains at the tran- 
sition is a power law with P{s) ~ s""^ exponent smaller 
than 2 (r « 1.6) for = 2 while the exponent is larger 
than 2 (t « 2.6 ) for any F > 2. I n one-dimensional 
systems instead (jKlemm et all l2003c( ) , the transition is 
continuous for all values of F. 

It is worth remarking that, upon interaction, the over- 



FIG. 10 Plot of the density of active links na{t) for 
= 10, L = 150 and (top to bottom) q = 
1, 100, 200, 230, 240, 250, 270, 300, 320, 400, 500, 10000. The 
inset reports the dependence of the freezing time tco on L 
for F = 10 and q = 100 < qc- The bold line has slope 2. 
From (jCastellano ef aZ.ll2000ll . 



lap between two neighbors always increases by but 
the change of a trait in an individual can make it more 
dissimilar with respect to his other neighbors. Hence, 
when the number of neighbors is larger than 2, each in- 
teraction can, somewhat paradoxically, result in an in- 
crease of the general level of disorder in the system. This 
competition is at the origin of the nontrivial temporal be- 
havior of the model in d = 2, illustrated in Fig. [101 below 
the transition but close to it (g < qc) the density of active 
links (connecting sites with overlap different from and 
1) has a highly non monotonic behavior. 

Most investigations of the Axelrod model are based 
on numerical simulations of the model dynamics. Ana- 
lytical app roaches are just a few. A simple mean field 
treatment (ICastellano et "aLl. l2000l:IVaia uez and Rednej . 
120071 : IVilone et all I2002D consists in writing down rate 
equations for the densities Pm of bonds of type m, i.e., 
connecting individuals with m equal and F — m differ- 
ent features. The natural order parameter in this case is 
the steady state number of active links Ua = X]m=i 
that is zero in the disordered phase, while it is finite in 
the ordered phase. This approach gives a discontinuous 
transition for any F. In the particular case of -F = 2 the 
mean field equations can be studied analytically in de- 
tail (fVazqucz an d Rcdncr. 2007), providing insight into 



^ Since q is discrete calling the transition "continuous" is a slight 
abuse of language. We will adopt it because the transition is 
associated with the divergence of a length, as in usual transitions. 



the non-monotonic dynamic behavior for <? < (?c and 
showing that the approach to the steady state is gov- 
erned by a timescale diverging as \q — qc\~^^'^ ■ Some 
information about the behavior of the Axelrod model for 
F = 2 and g = 2 is obtained also by a mapping to the 
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FIG. 11 Order parameter {Sn-iayi)/N as a function of the ef- 
fective noise rate r' = r(l — 1/q) for different values of q. 
Simulations have b een run in system s of size iV = 50^ with 
F = 10. From (,Klemm et am2003al 'l. 



constrained voter model ( Vazquez and Redneij . l20(3 ) dis- 
cussed in Sec. IIII.BI 



B. Variants of Axelrod model 

In his seminal paper, Axelrod himself mentioned many 
possible variants of his model, to be studied in order to in- 
vestigate the effect of additional ingredients as the topol- 
ogy of the interactions, random noise, the effect of mass 
media and many others. Over the years this program has 
been followed by many researchers. 

The possibility of one individual to change sponta- 
neously one of his traits, independently of his neighbor- 
hood, is denoted as "cultural drift" in social science and 
corresponds to the addition of flipping ev ents driven by 
random noise. In fKlem m et al I l2003aD it is demon- 
strated that the inclusion of noise at rate r has a pro- 
found influence on the model, resulting in a noise-induced 
order-disorder transition, practically independent of the 
value of the parameter q (Fig. [TT|). 

For small noise the state of the system is monocultural 
for any q, because disordered configurations are unsta- 
ble with respect to the perturbation introduced by the 
noise: the random variation of a trait unfreezes in some 
cases the boundary between two domains leading to the 
disappearance of one of them in favor of the other. How- 
ever, when the noise rate is large, the disappearance of 
domains is compensated by the rapid formation of new 
ones, so that the steady state is disordered. The thresh- 
old between the two behaviors is set by the inverse of 
the average relaxation time for a perturbation T(N), so 
that the transition occurs for rcT{N) = 0(1). An ap- 
proximate evaluation of the relaxation in = 2 gives 
T = N\ti{N), in good agree ment with simulation s, while 
T ^ N on one dimension ( Klemm et "all. I2005D . bmce 
T{N) diverges with N , the conclusion is that, no matter 
how small the rate of cultural drift is, in the thermody- 
namic limit the system remains always disordered for any 
9- 



The discovery of the fragility of the Axelrod model 
with respect to the presence of noise immediately raises 
the question "What is the simplest modification of the 
original model that preserves the exi stence of a transi- 
tion in the presence of noise?". In ( Kupermanl . |2006| ) 
two modified Axelrod-like dynamics have been intro- 
duced, where the interaction between individuals is also 
influenced by which trait is adopted by the majority of 
agents in the local neighborhood. Similar ingredients are 
present in two ot her variants of th e Axelrod model re- 
cently proposed ( Flache and Macvl . [2007). A convinc- 
ing illustration that these modifications lead to a robust 
phenomenology with respect to the addition of (at least 
weak) noise is still lacking. 

Another variant of the original definition of the model 
is the introduction of a threshold such that, if the over- 
lap is smaller than a certain va lue 9, no interaction takes 
place ( Flache and Macvl . l2007t ) . Unsurprisingly no qual- 
itative change occurs, except for a reduction of the or- 
dered region of the phase diagram ( De Sanctis and Gallal . 
l2007f ). Another possibility named "interaction noise", is 
that for u smaller than the threshold, the interaction 
takes place with probability 5. This kind of noise favors 
ordering but aga in does not lead t o drastic changes of the 
model behavior ( De Sanctis and Galla. 20(53) ■ 

In order to understand the effect of complex in- 
teraction topologies on its behavior, the Axelrod 
model has been studied on small-world and scale-free 
graphs ( Klemm et aLl . I2b03b( ) . In the first case, the tran- 
sition between consensus and a disordered multicultural 
phase is still observed, for values of the control parame- 
ter qc that grow as a function of the rewiring parameter 
p. Since the WS network for p = 1 is a random net- 
work (and then practically an infinite-dimensional sys- 
tem) this is consistent with the observa tion of the transi- 
tion also in the mean field approa ches ( Castellano et al\ . 
120001 : IVazg uez and Rednerj . [2iD07tl . The scale-free nature 
of the BA network dramatically changes the picture. For 
a given network of size N a somewhat smeared-out tran- 
sition is found for a value gc, with bistability of the 
order parameter, the signature of a first-order transi- 
tion. However numerical simulations show that the tran- 
sition threshold grows with iV as ~ ]V° '^^, so that 
in the thermodynamic limit the transition disappears 
and only ordered states are possible. This is similar to 
what occurs for the Ising model on scale-free networks, 
wher e the transition te mperature diverges with system 
size ( Leone et al\ . |2002| ). 

Another natural modification of the original Axelrod 
model concerns the effect of media, represented by some 
external field or global coupling in the system. One 
possible way to implement an external field consists in 
defining a mass media cultural message as a set of fixed 
varia bles M — (/ii, /i2, . . . , /if) (jGonzalez-Avella et al\ . 
I2005D . With probability B the selected individual in- 
teracts with the external field M exactly as if it were a 
neighbor. With probability 1 — B the individual selects 
instead one of his actual neighbors. Rather unexpect- 
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edly the external field turns out to favor t he multicultura l 
phase , in agreement with early findings ( Shibanai et all 
l200lh . The order-disorder transition point is shifted to 
smaller values of the control parameter qc{B). For B 
larger than a threshold such that qc{B*) — only the 
disordered phase is present: a strong external field favors 
the alignment of some individuals with it, but it simul- 
taneously induces a decoupling from individuals too far 
from it. 

Similar conclusions are drawn when a global cou- 
pling or a local non-uniform coupling are consid- 
ered (jGonzalez-Avella et aZ.1 . 120061) . In all cases the or- 
dered region of the phase diagram is reduced with respect 
to the case of zero field and it shrinks to zero beyond a 
certain field strength B* . Interestingly, for q > qc{B ^ Q) 
a vanishing field has the opposite effect, leading to an 
ordered monocultural state. The limit £? ^ is there- 
fore discontinuous. The same type of behavior is also 
found for indirect mass-media feedback, i.e., when sites 
accept the change of a trait only with probability R, if 
the new value of the trait i s not the same of the major- 
ity (jGonzalez-Avella et all . l2007f ) . 

In the Axelrod model the numerical value of traits is 
just a label: nothing changes if two neighbors have traits 
that differ by 1 or by 9 — 1. In order to model situa- 
tions whe re this difference actuall y matters, it has been 
proposed ( Flache and Macvl . |2006[ ) to consider some fea- 



tures to be "metric", i.e., such that the contribution to 
the overlap of a given feature is [1 — Acr//(g — l)]/F, 
where Aery is the difference between the trait values. In 
this way the Axelrod model becomes similar to the vecto- 
rial version of the Deffuant model. Although a systematic 
investigation has not been performed, it is clear that this 
variation favors the reaching of consensus, because only 
maximal trait difference (Act/ = q—1) totally forbids the 
interaction. A related variation with "me tric" features is 
described in (iDe Sanctis and Gallal . [2007h . 

Other recent works deal with a version of the 
Axelr od model with both an external field and 
noise ( Mazzitello et aLl . l2007l ) and one where individuals 
above a fixed threshold do not interact ( Parravano et all . 
I2OO6D . 



C. Other multidimensional models 

At odds with the detailed exploration of the behavior 
of the Axelrod model, much less attention has been paid 
to other types of dynamics for vectors of opinions. 

In t he original paper on the Deffuant 
model (jPeffuant et all I2OOOI ). a generalization to 



vectorial opinions is introduced, considering in this case 
binary variables instead of continuous ones. This gives 
a model similar to the Axelrod model with q = 2 traits 
per feature, with the difference that the probability of 
interaction between two agents as a function of their 
overlap is a step function at a bounded confidence 
threshold d. In mean field a transition between full 



consensus for large threshold and fragmentation for 
small d is found. 

A similar model is studied in ( Laguna et aLl . fJoOSt ) . In 
this case, when two agents are sufficiently close to inter- 
act, each pair of different variables may become equal 
with a probability /i. Again a transition between con- 
sensus and fragmentation is found as a function of the 
bounded confidence threshold, but its properties change 
depending on whether ^ — 1 or fi < 1. 

A generalization of continuous opinions (the HK 
mo del) to the vecto r ial (2- dimensional) case is reported 
in ( Fortunato et all l2005l ) for a square opinion space, 
with both opinions ranging between and 1, and square 
or circular confidence ranges. Assuming homogeneous 
mixing and solving the rate equations, it turns out that 
no drastic change occurs with respect to the ordinary 
HK model. The consensus threshold is practically the 
same. When there is no consensus the position of coex- 
isting clusters is determined by the shape of the opinion 
space. An extension of Deffuant an d HK models to vecto - 
rial opinions has been proposed in ( Lorend . l2006l l2008at ) . 
Here opinions sit on a hypercubic space or on a simplex, 
i.e., the components of the opinion vectors sum up to 
one. It turns out that consensus is easier to attain if the 
opinion space is a simplex rather than hypercubic. 

Other vectorial models are considered in Section 1X1 on 
the coevolution of networks and states 



V. LANGUAGE DYNAMICS 



La nguage dynamics is an emerging field ( Wichmannl . 
I2OO8D that focuses on all processes related to emergence, 
change, evolution, interactions and extinction of lan- 
guages. In this section we shall review some of the main 
directions in which this field is evolving. 

Models for language dynamics and evolution can be 
roughly divided in two main categories: sociobiological 
and sociocultural approaches. This distinc tion somehow 
paralle l s the debate nature versus nurture (lGaltonl . fl87l 
iRidlevl . [200! which concerns the relative importance of 
an individual's innate qualities (" nature" ) with respect 
to personal experiences ("nurture") in determining or 
causing individual differences in physical and behavioral 
traits. 

The sociobiolo g ical approach ( Hurfordl . Il989l : 
iPinker and Blooiiil Il990l ) postulates that successful 
communicators, enjoying a selective advantage, are 
more likely to reproduce than worse communicators. 
Successful communication contributes thus to biological 
fitness: i.e., good communicators leave more offspring. 
The most developed branch of research in this area 
is represented by the evolutionary approaches. Here 
the main hypothesis is that communication strategies 
(which are model-dep endent) are inna te, in the spirit of 
the nativist approach ( Chomskvl . 119651 ). and transmitted 
genetically across generations. Thus if one of them is 
better than the others, in an evolutionary time span it 



28 



will displace all rivals, possibly becoming the unique 
strategy of the population. The term strategy acquires 
a precise meaning in the context of each particular 
model. For instance, it can be a strategy for acquir- 
ing the lexicon of a language, i.e., a function from 
samplings of observed be haviors to acquired commu- 



nicative behavior patterns (jHurford", '1989^; 'Nowa k_et all 
Il999bl : lOliphanti 1 19971 : [oliphant and Batali, 19961) . or 



it ca n simply coincide w i th th e lexicon of the par- 
ents ( Nowak and Krakaueii [l999l ) . but other possibilities 
exist (|Steelsl . l2Q05l l 

On the other hand, in sociocultural approaches lan- 
guage is seen as a complex dynamical system that evolves 
and self-o rganizes, continuously sh aped and reshaped by 
its users ( Steels and Baillij . l200Cl( ). Here good strate- 
gies do not provide higher reproductive success but only 
better communication abilities. Agents can select better 
strategies exploiting cultural choice and direct feedback 
in communications. Moreover, innovations can be intro- 
duced due to the inventing ability of the agents. Thus, 
the study of the self-organization and evolution of lan- 
guage and meaning has led to the idea that a community 
of language users can be seen as a complex dynamical 
system which collectively solves the problem of develop- 
ing a shared communication system. In this perspective, 
which has been adopted by the novel field of semiotic 
dynamics, the theoretical tools developed in statistical 
physics and complex systems science acquire a central 
role for the study of the self-generated structures of lan- 
guage systems. 



A. Evolutionary approaches 



Ac c ording t o the s ociobiological a pproac h (iHurford 
19891: iNowakl 120061: iNowak et all Il999bl : lOliphant . 
1997 : lOliphant and Bat's!! . 1996| l. evolution is the main 



responsible both for the origin a nd the emergence 
of n atural language in humans ( Pinker and Blooml . 
Il990l ). Consequently, natural selection is the fun- 
damental driving force to be introd uced i n mod- 
els. Evolutionary game theory (ISmithl . I1982D 
was formulated with the aim of adapting classi- 
cal game theory (Ivon Neumann and Morgensternl Il947l : 



lOsborne and Rubinsteinl . 19941 ) to deal with evolution- 
ary issues, such as the possibility for agents to adapt, 
learn and evolve. The approach is phenotypic, and the 
fitness of a certain phenotype is, roughly speaking, pro- 
portional to its diffusion in the population. Strategies 
of classical game theory are substituted by traits (ge- 
netic or cultural), that are inherited, p ossibly with mu- 
tations. The search for Nash equilibria ( Nashl . [T95(]| ) be- 
comes the quest for evolutionary stable strategies. A 
strategy is stable if a group adopting it cannot be in- 
vaded by another group adopting a different strategy. Fi- 
nally, a fundamental assumption is that the payoff from 
a game is interpreted as the fitness of the agents involved 
in the game. The Evolutionary Language Game (ELG) 



(jNowak and Krakaueii Il999l : iNowak et all . Il999bl ) aims 
at modeling the emergence of language resorting to evo- 
lution ary game theory and to the concept of language 
game (IWittgensteinl . Il953allbl). For a recent experimen- 
tal paper we refer to ( Lieberman et a/.l . l2007l ). 



1. Evolutionary language game 

In this section we analyze in some detail how the prob- 
lem of the evolution of a common vocab ulary [or mor e 
generally a common set of conventions ( Lewis! . Il969t) . 
syntactic or grammatical rules] is addressed in the frame- 
work of evolutionar y game theory. The formalism we 
use is mutuated by ( Nowak et aLl . ll999bl ). but the basic 
struct ure of the game was already included in the seminal 
paper (lHurfordLll989D about the evolution of Saussurean 
signs (Ide Saussurd.ll9ld) 

A population of agents (possibly early hominids) lives 
in an environment with n objects. Each individual is able 
of produce a repertoire of m words (sounds or signals, in 
the original terminology) to be associated with objects. 
Individuals are characterized by two matrices P and Q, 
which together form a language L. The production ma- 
trix P is a n X rii matrix whose entry, py , denotes the 
probability of using word j when seeing object z, while 
the comprehension matrix Q is a mxn matrix, whose en- 
try, Qji, denotes the probability for a hearer to associate 
sound j with object i, with the following normalization 
conditions on the rows J2]LiPij — 1 ^^^d J27=i Iji — 1- 

A pair of matrices P and Q identifies a language L. 
Imagine then two individuals I\ and I2 speaking lan- 
guages L\ (defined by P\ and Qi) and Li (defined by P2 
and Qi)- The typical communication between the two 
involves the speaker, say /i, associating the signal j to 
the object i with probability pij . The hearer I2 infers ob- 



If one sums over all 



ject i with probability X^Jli Plj'lji 
the possible objects, one gets a measure of the ability, for 

/i, to convey information to h: J2i=i J2T=i Pif l-H ■ ^ 
symmetrized form of this expression defines the so-called 
payoff function, i.e., the reward obtained by two individ- 
uals speaking languages Li and I/2 when they communi- 
cate: 



F{Li,L2 



iT.T.i^^+^'ff)- (33) 

1=1 j=l 



From the definition of the payoff it is evident that each 
agent is treated once as hearer and once as speaker and 
they both receive a reward for successful communication. 

The crucial point of the model is the definition of the 
matrices P and Q which have to be initialized in some 
way at the beginning of the simulation. In principle there 
is no reason why P and Q should be correlated. On the 
other hand the best possible payoff is obtained by choos- 
ing P as a binary matrix having at least one 1 in every 
column (if n > m) or in every row (if n < m) and Q as a 
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binary matrix with qji = 1 , if pij is the largest entry in a 
column of P. If ri = m the maximum payoff is obtained 
for P having one 1 in every row and column and Q being 
the transposed matrix of P. In general the maximum 
payoff is given by Fmax — min{m,n}. It is also worth 
noting that the presence of two completely uncorre- 
lated matrices for the producti on, P, and comprehension, 
Q, rnodes, already present in (iHurfordl . il989 : Olip hant , 
119971: lOliphant and Bat alii . [l99^ could lead to patho- 
logica l situations as remarked in ( Komarova and Nivogil . 
I2004D . where a single matrix is adopted for both tasks. 

In a typical situation one simulates a population of 
N individuals speaking N different languages Li to Ln 
(by randomly choosing the matrices Pk and Qk, for 
k = 1, N). In each round of the game, every individ- 
ual communicates with every other individual, and the 
accumulated payoffs are summed up, e.g. the payoff re- 
ceived by individual k is given by Fk = X]il=i ^i^k, Li), 
with I k. As already mentioned, the payoff is inter- 
preted as fitness. In a parental learning scheme each 
individual will produce an offspring (without sexual re- 
production) with the probability fk — Fk / X); ^i- this 
way each individual gives rise on average to one offspring 
for the next generation and the population size remains 
constant. The individuals of the new generation learn the 
language of their parents by constructing an association 
matrix A, whose element records how many times 
the individual has observed its parent associating object 
i and signal j in K different samplings. The production 
and comprehension matrices P and Q are easily derived 
from the association matrix A as: 



P^3 



1=1 



^ji/ ^ "y- 
1=1 



(34) 



The form of the matrix A clearly depends on K. In the 
limit K oo the offspring reproduces the production 
matrix of its parent and A = P. For finite values of 
K, learning occurs with incomplete information and this 
triggers mutations occurring in the reproduction process. 

An important observation is in order. In such a scheme 
the language of an individual, i.e., the pair (P, Q), deter- 
mines its fitness and, as a consequence, the reproduction 
rate. On the other hand what is inherited is not directly 
the language but a mechanism to learn the language 
which is language specific, i.e., a language acquisition 
device in the spirit of the nativist approach ( Chomsky! . 
I1965D . Therefore the traits transmitted to the progeny 
can be different from the language itself. 

This evolutionary scheme leads the population to con- 
verge to a common language, i.e., a pair of {P,Q) ma- 
trices shared by all individuals. The common language 
is not necessarily optimal and the system can often 
get stuck in sub-optimal absorbing states where syn- 
onymy (two or more signals associated to the same ob- 
ject) or homonymy (the same signal used for two or 
more objects) are present. The convergence proper- 
ties to an absorbing state depend on the population 



size N as well as on K, but no systematic analysis 
has been performed in this direction. Another inter- 
esting direction is related to the underlying topology 
of the game. What described so far corresponds to 
a fully connected topology, where each agent interacts 
with the whole population. It is certainly of interest ex- 
ploring different topological structures, more closely re- 
lated to the structure of social networks, as discussed 



in (iHauert et all [20 05: S zab o and Fathl . l2007t) (see also 



the recent (Tavares et aLi. i2007h on the model proposed 



in (Komarova et all 
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The model can then be e nriched by adding a probabi l- 
ity of errors in perception ( Nowak and Krakaueii Il999l ). 
i.e., by introducing a probability Uij of misinterpreting 
signal i as signal j. The terms Uij are possibly defined 
in terms of similarities between signals. The maximum 
payoff for two individuals speaking the same language is 
now reduced, hence the maximum capacity of informa- 
tion t ransfer. This result is r eferred to as linguistic erro r 
limit ( Nowak and Krakaueii Il999l : iNowak et al. I,[l999i): 
the number of distinguishable signals in a protolanguage, 
and therefore the number of objects that can be accu- 
rately described by this language, is limited. Increas- 
ing the number of signals would not increase the capac- 
ity of information transfer [it is worth mentioning here 
the interesting parallel between the formalism of evo- 
lutio nary language game with that of information the- 
ory (Plotki^^^d[Now^, |2003)]- A possible way out i s 
that of combining signals into words ( Smith et aZ.1 . 12003V 
opening the way to a potentially unlimited number of ob- 
jects to refer to. In this case it is shown that the fitness 
function can overcome the error limit, inc reasing expo- 



nentially with the length o f words ( Nowak and Krakaue^ 
119991: iNowak et aZ] . ll999aD . This is considered one of the 
possible ways in which evolution has selected higher or- 
der structur es in lan guage, e.g. syntax and g rammar. 
We refer to (|Nowakl . 120061: INowak and Krakaueii \l99^ 
for details about the higher stages in the evolution of 
language. 



2. Quasispecies-like approach 

The model described in the previous section can be cast 
in the framework of a determi nistic dyn amical system 
(see for a recent discussion (Traulsen et aLl . ,2005) and ref- 
erences therein). We consider again the association ma- 
trix A, a, n X m matrix whose entries, a^, are non-zero if 
there is an association between the object i and the signal 
j. In this case we consider a binary matrix with entries 
taking either the value or 1. The possible number of 
matrices A is then M = 2"™. This matrix is a l so den oted 
as the lexical matrix ( Komarova and Nowakl l200ll ). In 
a population of N individuals denote now with Xk the 
fraction of individuals with the association matrix Ak, 
with X)fc=i = 1- One can define the evolution of clS 
given by the following equation: 
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Hk 



i)Xk , 1 = l,..,M 



(35) 



where /; is the fitness of individuals with the asso- 
ciation matrix Ai (from now on individual I), fi = 
'^f, F(Ai, Ak)xk, with the assumption that Xk is the 
probability to speak with an individual k; (j) defines the 
average fitness of the population, cj) — fixi, while Qik 
denotes the probability that someone learning from an in- 
dividual with Ai will end up with Ak- The second term 
on the right-hand side keeps the population size constant. 

Eqs. ((35|) repre s ent a particular case of the quasis pecies 
equations (|E igCTl I1971I : [Eigen and Schusteij . Il979l ). The 
quasispecies model is a description of the process of Dar- 
winian evolution of self-replicating entities within the 
framework of physical chemistry. These equations pro- 
vide a qualitative understanding of the evolutionary pro- 
cesses of self-replicating macromolecules such as RNA or 
DNA or simple asexual organisms such as bacteria or 
viruses. Quantitative predictions based on this model 
are difficult, because the parameters that serve as input 
are hard to obtain from actual biological systems. In the 
specific case in which mutation is absent, i.e., Qij = 
if i ^ J, one recovers the so-c alled r e plicato r equations 
of evolutionary game theory (|Smithl . I1982D . which, it 
is worth recalling, are equivale nt to the Lotka-Volterra 
equat ions in M — 1 dimensions ()Hofbauer and Sigmundl . 
I1998D . 



B. Semiotic Dynamics approach 

Semiotic dynamics looks at language as an evolving 
system where new words and grammatical constructions 
may be invented or acquired, new meanings may arise, 
the relation between language and meaning may shift 
(e.g. if a word adopts a new meaning), the relation be- 
tween meanings and the world may shift (e.g. if new 
perceptually grounded categories are introduced). 



1. The Naming Game 

The Naming Game (NG) possibly represents the sim- 
plest example of the complex processes leading pro- 
gressively to the establishment of human-like languages. 
It was expressly conceived to explore the role of self- 
organ ization in the evolution of language ( Steelsl . Il995l 
11996') and it has acquired, since then, a paradigmatic 
role in the whole fie ld of semiotic dynamics. The original 
paper (|Steelsl . [T995h . focuses mainly on the formation of 
vocabularies, i.e., a set of mappings between words and 
meanings (for instance physical objects). In this con- 
text, each agent develops its own vocabulary in a ran- 
dom private fashion. But agents are forced to align their 
vocabularies, through successive conversation, in order 
to obtain the benefit of cooperating through communi- 
cation. Thus, a globally shared vocabulary emerges, or 



should emerge, as the result of local adjustments of in- 
dividual word-meaning association. The communication 
evolves through successive conversations, i.e., events that 
involve a certain number of agents (two, in practical im- 
plementations) and meanings. It is worth remarking that 
here conversations are particular cases of language garn es, 
which, as already pointed out in (I Wittgensteinl . [l953al lbl) . 
can be used to describe linguistic behavior, even if they 
can include also non-linguistic behavior, such as pointing. 

This original seminal idea triggered a series of con- 
tributions along the same lines and many variants have 
been proposed over the years. It is par ticularly interest - 
ing to mention the work proposed in ( Ke et ai\ . |2002( ). 
that focuses on an imitation model which simulates how 
a common vocabulary is formed by agents imitating each 
other, either using a mere random strategy or a strategy 
in which imitation follows the majority (which implies 
non-local information for the agents). A further contri- 
bution of this paper is the introduction of an interaction 
model which uses a probabilistic representation of the 
vocabulary. The probabilistic scheme is formally simi- 
lar to the framework of evolutionary game theory seen 
in Sec. IV. A. 11 since to each agent a production and a 
comprehension matrices are associated. Differently from 
the approach of ELG, here the matrices are dynamically 
transformed according to the social learning process and 
the cultural trans mission rule. A simi lar approach has 
been proposed in ( Lenaerts et a/.l . l2005l ). 

In the next section we shall present a minimal ver- 
sion of the NG which results in a drastic simplification 
of the model definition, while keeping the same overall 
phenomenology. This version of the NG is suitable for 
massive numerical simulations and analytical approaches. 
Moreover the extreme simplicity allows for a direct com- 
parison with other models introduced in other frame- 
works of statistical physics as well as in other disciplines. 



a. The Minimal Naming Game Th e simplest version of 
the NG ( Baronchelli et all l2006b[ ) is played by a popu- 
lation of N agents, on a fully connected network, trying 
to bootstrap a common name for a given object. Each 
player is characterized by an inventory of word-object 
associations he/she knows. All agents have empty inven- 
tories at time t — Q. At each time step [t — 1,2,..), two 
players are picked at random and one of them plays as 
speaker and the other as hearer. Their interaction obeys 
the rules described in Fig. [T^ 



b. Macroscopic analysis The first property of interest is 
the time evolution of the total number of words owned by 
the population Nw(t), of the number of different words 
Nd{t), and of the success rate S{t) (Fig. fTS]). 

After a transient period, the system undergoes spon- 
taneously a disorder-order transition to an asymptotic 
state where global consensus emerges, i.e., every agent 
has the same word for the same object. It is remark- 
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FIG. 12 Naming Game. Examples of the dynamics of the 
inventories in a failed (top) and a successful (bottom) game. 
The speaker selects the word highlighted. If the hearer does 
not possess that word he includes it in his inventory (top). 
Otherwise both agents erase their inventories only keeping 
the winning word (bottom). 
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FIG. 13 Naming Game, a) Total number of words present in 
the system, Niu{t); b) Number of different words, Nd{t); c) 
Success rate S{t), i.e., probability of observing a successful in- 
teraction at time t. The inset shows the linear behavior of S{t) 
at small times. The system reaches the final absorbing state, 
described by N^{t) = N, Nd{t) = 1 and S (t) = 1, in which a 
global agreement has been reached. From (jBaronchelli et all 
l2006bl ). 



able that this happens starting from completely empty 
inventories for each agent. The asymptotic state is one 
where a word invented during the time evolution takes 
over with respect to the other competing words and im- 
poses itself as the leading word. In this sense the sys- 
tem spontaneously selects one of the many possible co- 
herent asymptotic states and the transition can thus be 
seen as a symmetry breaking transition. The dynamics 
of the Naming Game is characterized by the following 
scaling behavior for the convergence time tconv, the time 
and the height of the peak of Nu,(t), namely tmax and 
]\[max _ N'^(t^^^y It turns out that all these quantities 



Nmax N'^ and tdiff = [tconv - Unax) ~ , with 

exponents a — (3 — ^ — 5— 1.^ (with a subtle fea- 
ture around the disorder-order transition where an ad- 
ditional timescale emerges). The values of those expo- 
nents can be understood throu gh simple scaling argu- 
ments ( Baronchelli et a/] . l2006bf )^. 



2. Symmetry breaking: a controlled case 

Consider now a simpler case in which there arc only 
two words at the beginning of the process, say A and B, so 
that the population can be divided into three classes: the 
fraction of agents with only A, ua, the fraction of those 
with only the word B, ub, and finally the fraction of 
agents with both words, uab- Describing the mean-field 
time evolution of the three species is straightforward: 
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(36) 



-2nAnB - 2nAB - ("A + nB)nAB- 



The system of differential equations (|36p is determin- 
istic. It presents three fixed points in which the sys- 
tem can collapse depending on the initial conditions. If 
nA{t = 0) > nB{t = 0) [nB{t = 0) > ^^(t = 0)], at the 
end of the evolution we will have the stable fixed point 
UA — 1 [nB = 1] and, consequently ub = uab — 
[nA = nAB = 0]. If, on the other hand, we start 
from nA{t = 0) = nB{t = 0), the equations lead to 
71 A — nB ^ '^.UAB = 0.4. The latter situation is clearly 
unstable, since any external perturbation would make the 
system fall in one of the two stable fixed points. 

Eqs. ([55)1 however, are not only a useful example to 
clarify the nature of the symmetry breaking process. In 
fact, they also describe the interaction among two differ- 
ent populations that converged separately to two distinct 
conventions. In this perspective, Eqs. (j36p predict that 
the larger population will impose its conventions. In the 
absence of fluctuations, this is true even if the difference 
is very small: B will dominate if nB(i — 0) = 0.5 -I- e 
and UAit = 0) = 0.5 — e, for any < e < 0.5 and 
nAB{t = 0) = 0. Data from simulations show that the 
probability of success of the convention of the minor- 
ity group ua decreases as the system size increases, go- 
ing to zero in the thermodynamic limit {N — > cx)). A 
similar approach has been propo sed to model the com- 
petit ion between two languages ( Abrams and Strogatj . 
I2003D . We discuss this point in Sec. ITD] Here it is 
worth remarking the formal similarities between model- 
ing the competition between synonyms in a NG frame- 
work and the competition between languages: in both 
cases a synonym or a language are represented by a sin- 
gle feature, e.g. the characters A or B, for instance. 



follow power law behaviors: t„ 



AT" f 



Here the time is the number of binary interactions 



in Eqs. (|36p . The similarity has been made more ev- 
ident by the subsequent variants of the model intro- 
duced in ( Abrams and Strogat j . I20d3h to include explic- 



itly th£_20Ssibilit;^_of bilingual i ndividuals. In partic ular 
in (jMinett and Wand . [2008: Wang and Minettl . l2005[ ) de- 



terministic models for the competition of two languages 
ha ve been prop osed, which include bilingual individuals. 
In (jCastello et ai, 2006) a modified version of the voter 
model (see Sec. IIII.BI) including bilingual individuals has 
been proposed, the so-called AB-model. In a fully con- 
nected network and in the limit of infinite population 
size, the AB-model can be described by coupled differ- 
ential equations for the fractions of individuals speaking 
language A, B or AB, that are, up to a constant normal- 
ization factor in the timescale, identical to Eqs. (l36|) . 



3. The role of the interaction topology 



As already mentioned in Sec. III. Bl social networks play 
an important role in determining the dynamics and out- 
come of language change. The first investigation of the 
role of topology was proposed in 2004, at the 5th Con - 
ference on Language evolution, Leipzig (! Ke et aLl . l2008l ). 
Since then many approaches focused on adapting known 
models on topologies of increasing complexity: regular 
lattices, random graphs, scale-free graphs, etc. 

The NG, as described above, is not unambiguously 
defined on general networks. As already observed in 
Sec. llLBl and in Sec. UlLBl when the de gree distribution 
is heterogeneous, the order in which an agent and one 
of its neighbors are selected does matter, because high- 
degree nodes are more easily chosen as neighbors than 
low-degree nodes. Several variants of the NG on generic 
networks can be defined. In the direct NG [reverse NG) 
a randomly chosen speaker (hearer) selects (again ran- 
domly) a hearer (speaker) among its neighbors. In a 
neutral strategy one selects an edge and assigns the role 
of spe aker and hearer with eq ual probability to the two 
nodes (jDalPAsta et aLl . l2006bf ). 

On low-dimensional lattices consensus is reached 
throug h a coarsening phenomenon (|Baronchelli et all . 
l2006ah with a competition among the homogeneous clus- 
ters corresponding to different conve ntions, driven by 
the curvature of the interfaces (|Bravl .[l994). A scaling 
of the convergence time as 0{N'^~^^^'^) has been con- 
jectured, where c? < 4 is the lattice dimension. Low- 
dimensional lattices require more time to reach con- 
sensus compared to a fully connected graph, but a 
lower use of memory. A si milar analysis h as be en per- 
formed for the ^B-model (|Castell6 et all l2006f) . The 
effe ct of a sm all-world topology has been investigated 
in (iDall'Asta et al. . 2006a) in the framework of the NG 
and in (|Castell6 et ali l2006l ) for the AB-model. Two 
different regimes are observed. For times shorter than 
a crossover time, t cross = 0{N/p'^), one observes the 
usual coarsening phenomena as long as the clusters are 
one-dimensional, i.e., as long as the typical cluster size 
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FIG. 14 AB-model. Time evolution of the average density 
(p) of bilingual individuals in small- world networks for differ- 
ent values of the rewiring parameter p. From left to right: 
p = 1.0,0.1,0.05,0.01,0.0. The inset shows the dependence 
of the characteristic lifetime r on the rewiring parameter p. 
The d ashed line correspon ds to a power law fit r ~ p""'^®. 
From (|Castell6 et aZ.I. |2006| ) . 



is smaller than 1/p. For times much larger than t cross, 
the dynamics is dominated by the existence of short- 
cuts and enters a mean field like behavior. The con- 
vergence time is thus expected to scale as N^^^ and not 
as (as in c? = 1). Small- world topology allows to 
combine advantages from both finite-dimensional lattices 
and fully connected networks: on the one hand, only a 
finite memory per node is needed, unlike the ©(TV^/^) 
in fully connected graphs; on the other hand the con- 
vergence time is exp ected to be much sh orter than in 
finite dimensions. In ( Castello et aZ.I . [2006t) . the dynam- 
ics of the yl_B-model on a two-dimensional small world 
network has been studied. Also in this case a dynamic 
stage of coarsening is observed, followed by a fast decay 
to the A or B absorbing states caused by a finite size 
fluctuation fFig. [Ti)) . The NG has been studied on com- 
plex networks as well. Here the convergence time tconv 
scales as N^, with c:^ 1.4 ± 0.1, for both Erdos-Renyi 
(ER) (Erdos and Rcnvi, 1959,1960") and Barabasi- Albert 
(BA) ( Barabas i and Albert, 1999. ) networks. The scal- 
ing laws observed for the convergence time are general 
robust features not affected by further topologic al de- 
tails fo airAsta and BaronchcUi, 2006; DairA sta et all 
l2006br ). Finally it is worth mentioning how the Nam- 
ing Games with local broadcasts on random geomet- 



ric graphs have been investigated in ( Lu et all l2008l ) 
as a model for agreement dynamics in large-scale au- 
tonomously operating wireless sensor networks. 
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C. Other models 

A variant of the NG has been introduced with the 
aim of mimicking the mechanisms leading to opinion 
and convention forma t ion in a population of individu- 
als ( Baronchelli et all . l2007f ). In particular a new pa- 
rameter, (3, has been added mimicking an irresolute at- 
titude of the agents in making decisions (/? = 1 cor- 
responds to the NG). The parameter /3 is simply the 
probability that, in a successful interaction, both the 
speaker and the hearer update their memories erasing 
all opinions except the one involved in the interaction 
(see Fig. [12]). This negotiation process displays a non- 
equilibrium phase transition from an absorbing state in 
which all agents reach a consensus to an active (not frozen 
as in the Axelrod model ( Axelrodl . ll997( )) stationary state 
characterized either by polarization or fragmentation in 
clusters of agents with different opinions. Interestingly, 
the model displays the non-equilibrium phase transition 
also on heterogeneous networks, in contrast with other 
opinion dynamics models, lik e for instance the Axelrod 
model (jKlemm et all l2003bl ). for which the transition 
disappears for heterogeneous networks in the thermody- 
namic limit. 

A hybrid approach, combining vertical and hori- 
zontal tran smission of cultu ral traits, has been pro- 
posed in (|Ke et all |2002|) . while an evolutionary 



ver sion of the Naming Game has been introduced 
in ( Lipowski and Lipowskal [2008t) . 

A very interes ting approach to language chan ge has 
been proposed in ( Baxter et a/.l . l200"6l : ICroftll200Cll ). based 
on an utterance selection model. Another interesting 
approach to language dynam i cs, th e so-called Iterated 



Learning Model (ILM) (iKirbvl. 
evolution and learning (jNivog: 



2001 V, focuses on cultural 



20061 ) and explores how 



the mappings from meanings to signals are transmitted 
from generation to generation. In this framework several 
results have been obtained concerning the e mergence of a 
lingu istic structure, e.g. compositionality ( Smith et all 
I2003D . 



D. Language competition 

Models of language evolution usually focus on a sin- 
gle population, which is supposed to be isolated from the 
rest of the world. However, real populations are not iso- 
lated, but they keep interacting with each other. These 
steady interactions between peoples play a big role in the 
evolution of languages. 

At present, there are about 6,500 languages in the 
world, with a very uneven geographic distribution. Most 
of these languages have very few speak ers, and are threat- 
ened by extinction ( Sutherland! . |2003( ) . Indeed, it is plau- 
sible that in the future increasing numbers of people 
will be pushed to adopt a common language by socio- 
economic factors, leading to the survival of a few ma- 
jor linguistic groups, and to the extinction of all other 
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FIG. 15 Distribution of language sizes. The i-axis repre- 
sents the number of individuals speaking a language. In 
the upper diagram, on the y axis the number of languages 
spoken by x individuals is reported. In the lower dia- 
gram the number of languages is divided by the size of the 
corresponding bins, obtaining the probability distribution, 
well fitted by a log-normal function (continuous line). The 
total number of languages sampled is 6142. (Data from 
http: //www . ethnologue . com/ ). 



languages. According to some estimates, up to 90% of 
present languages might disappear by the end of the 21st 
century (Krauss . ,1992.) . The histogram of the language 
sizes, shown in Fig. llSi has a regular shape, which closely 
resembles a log- normal distribution^. Several models of 
language competition have been proposed with the aim 
of reproducing such distribution. However, we stress that 
the observed distribution of language sizes may not be a 
stable feature of language diversity, as there is no reason 
to believe that it has kept its shape over the past cen- 
turies and that it will keep it in the future. Morevoer, 
it has been suggested that the current histog ram may 
be a c onsequence of pure demographic growth ( Zanettd . 
l2007bD . 

Modeling language competition means studying the 
interaction between languages spoken by adults. Lan- 
guage evolution shares several features with the evolu- 
tion of biological species. Like species, a language can 
split into several languages, it can mutate, by modifying 
words/expressions over time, it can face extinction. Such 
similarities have fostered the application of models used 
to describe biological evolution in a language competi- 
tion context. The models can be divided in two cate- 
gories: macroscopic models, where only average proper- 
ties of the system are considered, are based on differential 



Similarfy, in l|Austoos and Petronil . l2007h . tfie distribution of the 
number of adfierents to religions has been investigated. 
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equations; microscopic models, where the state of each 
individual is monitored in time, are based on computer 
simulations. 



1. Macroscopic models 

The first macroscopic model of language competition 
was a dynamic model proposed in ( Abrams and Strogatj . 
I2OO3I ) (AS), describing how two languages, A and B, 
compete for speakers. The languages do not evolve in 
time; the attractiveness of each language increases with 
its number of speakers and perceived status, which ex- 
presses the social and economic benefits deriving from 
speaking that language. We indicate with x and with 
< s < 1 the fraction of speakers and the status of 
A, respectively. Accordingly, language B has a fraction 
y = 1—x of speakers, and its status is l—s. The dynamics 
is given by the simple rate equation 



dx 

— = c(l 
dt ^ 



x)x''s-cx{l~x)''{l-s), (37) 



where a and c are parameters which, along with s, fix 
the model dynamics^ . Eq. ([57)) expresses the balance be- 
tween the rates of people switching from language B to 
A and from A to B. The dynamics has only two sta- 
ble fixed points, corresponding to a; = and x — I. 
There is a third fixed point, corresponding to x = 1/2, 
s — 1/2, when the two languages are equivalent, but 
it is unstable, as confirmed by numerical simulations 
of a microscopic version of the AS m odel on different 
graph topologies ( Stauffer et all l2007l ). Therefore, the 
AS model predicts the dominance of one of the two lan- 
guages and the consequent extinction of the other. The 
dominant language is the one with the initial majority of 
speakers and/or higher status. Comparisons with empir- 
ical data reveal that the model is able to reproduce the 
decrease in time of the number of speakers of various en- 
dangered languages (Fig. . The AS model is minimal 
and neglects important aspects of sociolinguistic interac- 
tion. In actual situations of language competition, the 
interaction between two languages A and B often occurs 
thr ough speakers who are p roficient in both languages. 
In ( Mira and Paredei . 12005 1 bilingual speakers were in- 
troduced in the AS model. A parameter k expresses the 
similarity of the two competing languages A and B and 
is related to the probability for monolingual speakers to 
turn bilingual. For each choice of the AS parameters 
a, s, there is a critical value fcmi„(a,s) such that, for 
k > kmin{o-, s), the system reaches a steady state char- 
acterized by the coexistence of one group of monolingual 
speakers with a group of bilinguals. Monolingual speak- 
ers of the endangered language are bound to disappear. 



We remark that the parameter c is an overall multiplicative con- 
stant of the right-hand side of Eq. I|37|l and can be absorbed in 
the time unit, without affecting the dynamics. 
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FIG. 16 . Dynamics of language extinction according to the 
model of Abrams and Strogatz. The four panels show the 
comparison of the model with real data on the proportion 
of speakers over time for (a) Scottish Gaelic in Sutherland, 
Scotland, (b) Quechua in Huanuco, Peru, (c) Welsh in Mon- 
mouthshire, Wales and (d) Welsh in all of Wales. In (d) both 
historical data (squares) and the results of a recent cens us 
(triangles) are plotted. From ([Abrams and Strogatj . [20031 ') . 



but the survival of the language is ensured by bilingual- 
ism, provided A and B are similar enough. The model 
describes well historical data on the time evolution of the 
proportions of spe akers of Galician and Ca stillian Span- 
ish in Galicia. In ( Minett and Waii3 . l2008f ) a more com- 
plex modification of the AS model, incorporating bilin- 
gualism and language transmission between adults and 
from adults to children, was proposed. The model pre- 
dicts the same extinction scenario of the AS model, un- 
less special strategies of intervention are adopted when 
the number of speakers of the endangered language de- 
creases below a threshold. Effective intervention strate- 
gies are the enhancement of the status of the endangered 
language and the enforcement of monolingual education 

of the chil dren. 

In (Patriarca and Leppanenl |2004| ) the effect of popu- 
lation density is introduced, by turning the rate equation 
of the AS model into a reaction-diffusion equation. Here 
people can move on a plane, divided in two regions; in 
each regions one language has a higher status than the 
other one. The system converges to a stable configuration 
where both languages survive, although they are mostly 
concentrated in the zones where they are favored. In a 
recent work ( Pinasco and Romanellil . I2OO6D it was shown 



that language coexistence in the same region is possible, 
if one accounts for the population dynamics of the two 
linguistic communities, instead of considering the whole 
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population fixed, fike in the AS model. The dynamics 
is now ruled by a set of generalized Lotka-Volterra equa- 
tions, and presents a non-trivial fixed point when the 
rate of growth of the population of speakers of the en- 
dangered language compensates the rate of conversion to 
the dominant language. 



2. Microscopic models 

Many microscopic models of language competition rep- 
resent language as a set of F independent features {F 
usually goes from 8 to 64), with each feature taking one 
out of Q values. This is also the representation of culture 
in the Axelrod model (see Sec. IIV.A[) : indeed, language 
diversity is an aspect of cultural diversity. If Q = 2, lan- 
guage is a bit-string, a representation used for biological 
species (|Eigeij . Il97ll) . For a recent review of la nguage 
competition simulations see ( Schulze et al ], 120081). 

In the Schulze model ( Schulze and Stauffeil . l2005al ). 
the language of each individual evolves according to three 
mechanisms, corresponding to random changes, trans- 
fer of words from one language to another and learn- 
ing of a new language. There are three parameters: p, 
q and r. With probability p, a randomly chosen fea- 
ture of an agent's language is modified: with probabil- 
ity q, the new value is that of the corresponding fea- 
ture of a randomly picked individual, otherwise a value 



taken at random. Finally, there is a probability (1 — x) r 
that an agent switches to a language spoken by a frac- 
tion X of the population. If agents are the nodes of a 
network, the language of an individual can only be af- 
fected by its neighbors. Simulations show that there is a 
sharp transition between a phase in which most people 
speak the same language (low p), and a phase in which 
no language dominates (high p) and the distribution of 
language sizes is roughly log-normal, like the empirical 
distribution fFig. I15p . The agreement with the data im- 
proves by sampling t he evolving model dist ribution over 
a large time interval ( Stauffer et a/.l . [2006c[ ) . An analyt- 
ical formu lation of the Sc hulze model has been recently 
proposed (jZanettd . [2007al) . 

We notice that here languages have no intrinsic fit- 
ness, i.e., they are all equivalent for the dynamics, at 
variance with biological species and with the macroscopic 
models of the previous section, where the different status 
of languages is responsible for their survival/extinction. 
The possible dominance of one language is determined 
by initial fluctuations, that make a linguistic community 
slightly larger than the others and more likely to capture 
speakers fleeing from other communities. 

Several modifications of the Schulze model have 
been proposed. Agen t s can age, reproduce and 
die ( Schulze and Stauffeil . l2005al ): they can move on 
the sites of a lattice, forming linguistic communities 
that are spatially localiz ed ( Schulze and Stauffer, 20053); 
they can be bilingual (jSctiulze et ail I2008D . To avoid 
the dominance of a single language, it is enough to 



stop the flight from an endangered language when 
the number of its speakers dec reases below a thresh- 
old ( Schulze and Stauffeil . [200^. A hnguistic taxon- 
omy can be introduced, by classifying languages into 
families based on simila r ities of the corresponding bit 
strings ( Wichmann et a/l |2006( ). This enables to control 
the dynamics of both individual languages and of their 
families. 

The model in ( de Oliveira et all . l2006bl) describes the 
colonization of a territory by a population that eventu- 
ally splits into different linguistic communities. Language 
is represented by a number, so it has no internal struc- 
ture. The expansion starts from the central site of a 
square lattice, with some initial population size. Free 
sites are occupied by a neighboring population with a 
probability proportional to the number of people speak- 
ing that language, which is a measure of the fitness of 
that population. The language of a group conquering a 
new site mutates with a probability that is inversely pro- 
portional to its fitness. The simulation stops when all 
lattice sites have been occupied. The resulting linguis- 
tic diversity displays similar features as those observed 
in real linguistic diversity, like the distribution of lan- 
guage sizes (Fig. [T5|) . The agreement improves by in- 
trodu cing an upper bound for the fitness of a popula- 
tion ( de Oliveira et al. . 2006al). or by repres enting lan- 
guages as bit strings (|de Oliveira et 



Social Impact Theory (see Se c . IIII.DD was applied to 
model language change (jNettlel . Il999al lbl). Here, there 
are two languages and agents are induced to join the 
linguistic majority because it exerts a great social pres- 
sure. Language mixing, for which a new language may 
originate from the merging of two l anguages, was im- 
plemented in ( Kosmidis et all l2005t ). In this model, 
the biological fitness of the agents may increase if they 
learn words of the other language. The model ac- 
counts for the emergence of bilingualism in a community 
wh ere people initially speak only one of two languages. 
In ( Schwammlel . l2005h there are two languages and agents 



move on a lattice, are subjected to biological aging and 
can reproduce. People may grow bilingual; bilinguals 
may forget one of the two languages, if it is minoritar- 
ian in their spatial surroundings. As a result, if the two 
linguistic communities are spatially separated, they can 
coexist for a long time, before the dynamics will lead to 
the dominance of one of them. Bilingual agents are also 
present i n the modified versio n of the voter model pro- 
posed in ( Castello et aLl . I2OO6I ). discussed in Sec. IIII.Bl 



VI. CROWD BEHAVIOR 

Collective motion is very common in nature. Flocks 
of birds, fish schools, swarms of i nsects are among the 
most spectacular manifestations ( Parrish and Hamneil . 
Il997f ). Humans display similar behavior in many in- 
stances: pedestrian motion, panic, vehicular traffic, etc.. 

The origin of collective motion has represented a puz- 



36 



zle for many years. One has the impression that each 
individual knows exactly what all its peers are doing in 
the group and acts accordingly. It is plausible instead 
that an individual has a clear perception of what hap- 
pens in its neighborhood, ignoring what most of its peers 
are doing. We are then faced again with a phenomenon 
where local interactions determine the emergence of a 
global property of the system, in this case collective mo- 
tion. Therefore it is not surprising that in the last years 
physicists have worked in this field. In this section we 
shall give a brief account of the most important results 
on crowd behavior. For a review o f the studies on ve- 
hicular traffic w e refer to (jHelbintd . I2OOII : iKerneii |2004 
lNagataniLl2003 ). 



A. Flocking 

To study the collective motion of large groups 

of organisms, the concept of self-propelled p a rticles 

(SPP) has been introduced (ICzirok and Vicsekl . Il999t 



ICzirok and Vicsekl . I2OOOI : IVicsek et all Il995^ . SPP are 



particles driven by an intrinsic force, produced by an 
energy depot which is internal to the particles, just 
like it happens in r eal or ganisms. In the original 
model ( Vicsek et all Il995l ). N particles move on a 
squared surface with periodic boundary conditions. A 
particle i is characterized by its position Xi and veloc- 
ity Vj. For simplicity, it is assumed that the velocities 
of the particles equal in modulus the constant vq at any 
moment of the system's evolution. This is where the self- 
propelled feature of the particles sets in. The interaction 
is expressed by a simple rule: at each time step a particle 
i assumes the average direction of motion of the particles 
lying within a local neighborhood S{i), with a random 
perturbation, i.e. 



,(i + At) = (0(<))s(.)+C, 



(38) 



where 9i is the angle indicating the velocity direction of 
particle i, At the time step and the noise ^ is a ran- 
dom variable with a uniform distribution in the range 
[—77/2,77/2]. Noise is a realistic ingredient of the model, 
relaxing the strict prescription of alignment to the aver- 
age direction of motion of the neighboring peers. 

The model resembles the classical ferromagnetic XY 
model, where the velocities of the particles play the 
role of the spins and the noise that of tempera- 
ture ( Binnev et all . [l"992f ) . As in the XY model, spins 
tend to align. In the limit 7;o ^ the SPP model 
resembles Monte Carlo simulations of diluted XY fer- 
romagnets. However, when vq 7^ 0, it is a non- 
equilibrium model involving transport, as particles con- 
tinuously move on the plane. 

Initially, the directions of the velocity vectors are ran- 
domly assigned, so that there is no organized flow of par- 
ticles in the system. After some relaxation time, the 
system reaches a steady state. The order parameter is 
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FI G. 17 Velocity fields of particles in the model introduced 
in (jVicsek et al\ . Il995f l. A phase transition from a disor- 
dered (a) to an ordered (d) state of collective motion is 
observ ed by decreasing t he amount of noise in the system. 
From (jVicsek et ali 119951 ) . 



the magnitude $ of the average velocity of the system 

j 

Collective motion arises when $ > 0. The level of noise 
77 is the control parameter, just like temperature in spin 
models. Numerical simulations indicate that, in the ther- 
modynamic limit of infinite system size, $ varies with 77 
as follows 



$(77) 



riaip)-V 







for 77 < 77c (p) 
for 77 > 77c (p) 



(40) 



where p is the density of the particles. Eq. (|40|) indicates 
that the system undergoes a continuous kinetic transi- 
tion from a phase, where the initial rotational symmetry 
of the system is preserved, to a phase where it is bro- 
ken due to the emergence of collective motion (Fig. [17]) . 
The flow is generated for values of the noise lower than 
a threshold 77c (p), which depends on the particle density. 
Like in equilibrium phase transitions, the fluctuations of 
the order parameter diverge at the critical threshold. Nu- 
merical estimates yield (3 = 0.42±0.03, different from the 
mean field value 1/2. The critical line 77c(p) p'^, with 
K = 0.45 ± 0.05, and differs from the corresponding line 
for diluted ferromagnets: here, when 7/ ^ 0, the critical 
density vanishes; for diluted ferromagnets, instead, where 
the absence of noise corresponds to a vanishing tempera- 
ture, the critical density approaches the finite percolation 
threshold. 

These results hold for any non-vanishing value of the 
velocity vq. It is remarkable that the model has an 
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ordering phase transition in two dimensions, as it is 
known that the equiUbrium XY model does not or- 
der, due to the formation of Kost erlitz-Thouless vor- 
tices ( Kosterlitz and Th ouless'. 'l973h. Simulations show 



that, for vq w 0, vortices indeed form in the first steps 
of the evoluti on, but they are uii stable and disappear in 
the long run (I Czirok et al. I. I1997D . 

Similar SPP models were also int roduced in one and 
three dimensions. In one dimension ( Czirok et aLl . [l999l ) 
the interaction rule has to be modified: since there is 
only one direction of motion, imposing a constant par- 
ticle's velocity would yield a trivial dynamics, without 
contact to real situations. In a realistic dynamics, parti- 
cles change direction of motion after slowing down. This 
can be accomplished with a simple modification of rule 
psp . The resulting model displays the same kinetic tran- 
sition as in two dimensions, but with different exponents 
{(3 = 0.60±0.05, K « 1/4). Collective motion emerges in 
three dimensions as well dCzirok et at, 1999). Here the 
kinetic phase transition is less surprising, as in three di- 
mensions even equilibrium continuous spin models have 
a ferromagnetic phase transition. The critical exponents 
are found to be consistent with the mean field exponents 
in equilibrium systems. In the limit vq the model be- 
comes analo gous to the three-diin ensional diluted Heisen- 
berg model ( Binnev et aZ.l .ll992^. However, for 3D SPP 
collective motion can be generated for any value of the 
particle density p, whereas the static model cannot or- 
der for densities lower than the percolation threshold 

Pperc ~ 1 • 

Equilibrium systems cannot display an ordering phase 
transitio n in 2D because of the Mermin- Wagner theo- 
rem t Mermin and Wagnerl . [l96B), so the kinetic transi- 
tion observed in the SPP model is due to its intrinsic 
dynamical, non-equilibrium character. To explain the 
peculiar featur es of the model continuum theories have 
been proposed (jCzirok et all \l99^ . iToner and f^ . Il995l 
IT998) . where the relevant variables are the velocity and 
density fields v(x, i) and p{'x.,t). 

In a first approach, continuum hydrodynamic equa- 
tions of motion were used (, Toner and T u. 1995, 1998,;). 
The equation terms are dictated by the strict require- 
ment that rotational symmetry must be preserved 

afV + (vV)v = av-/3|v|2v- VP-|-i:>LV(Vv) 
dtp + V{pv) = 0, 

(41) 

where P is the analogue of pressure, a, /3 are coefficients 
ruling the symmetry breaking, the -Dl,i.2 are diffusion 
constants and S, is uncorrelated Gaussian random noise. 
The first equation contains some additional terms with 
respect to standard Navier-Stokes equations for a simple 
fluid, because the SPP model does not conserve the mo- 
mentum by construction, so Galilean invariance is bro- 
ken. The second equation expresses the conservation of 
the mass. An analysis of these equations based on the 
dynamic renormalization group ( Forster et aLl ■ [l977^ re- 
veals that the system orders in 2D, and that in more than 



four dimensions it behaves just like the corresponding 
equilibrium model, with equal exponents. On the other 
hand, there should be no phase transition in one dimen- 
sion, in contrast to numerical findings. For this reason, 
a different continuum theory was derived directly from 
the m aster equation of the ID SPP model ( Czirok et all 
I1999I ). Its equations read 



dtU = /([/) + M'a2[^ + ai^^^^^-f c 
dtp = -VQd^{pU) + Ddlp, 



(42) 



where U{x,t) is the velocity field, f{U) a self-propulsion 
term, C the noise. The crucial feature of Eqs. (|^^ is 
the existence of the non-linear term {dxU){dxp)/ p- A 
linear stability analysis of the equations reveals that there 
is an ordered phase if the noise is sufficiently low, as 
the domain walls separating groups of particles travelling 
in opposite directions are unstable. Numerical solutions 
of Eqs. (|42|) show that the continuum theory and the 
discrete ID SPP model belong to the same universality 
class. 



B. Pedestrian behavior 

Pedestri an behavior has been empir ically studied since 
the 1950s (|Hankin and WrightL Il958l). The first physi- 
cal modeling was proposed in (Henderso^, Il97l[ ) where 
it was conjectured that pedestrian fiows are similar to 
gases or fluids and measurements of pedestrian flows 
were compared with Navier-Stokes equations. However, 
realistic macroscopic models should account for effects 
like maneuvers to avoid collisions, for which energy and 
momentum are not conserved. Moreover, they should 
consider the "granular" structure of pedestrian fiows, 
as each pedestrian occupies a volume that cannot be 
penetrated by others. Therefore, micros copic models 
have recently attracted rn u ch at tention ( Galeal . 120031 : 
ISchreckenberg and Sharrnal . l200lh . One distinguishes 
two main approaches: models based on cellular automata 
(CA) and the social force model. 

In CA models of pedestrian dynamics fBlue and Adlei 
1998 , 2000; Burstedd_ej^;., 2001; Fukui and Ishibashj 



I999I: iManiccaml. . 
Muramatsu and Nagatanil! 



20031: iMuramatsu et all 
l2000l ). thM and 



|199£; 
space 

are discretized. The pedestrian area is divided into 
cells, which can be either empty or occupied by a 
single agent or an obstacle. A pedestrian can move 
to an empty neighboring cell at each time step. The 
motion of a single pedestrian is a biased random walk, 
where the bias is represented by a field residing on the 
space cells, which determines the transition rates of the 
agent towards neighboring cells, much like it happens 
in chemotaxis. CA models are computationally very 
efficient, but they do not describe well the complex 
phenomenology observed in real pedestrian dynamics, 
mostly because of space discretization, which constrains 
the directions of traffic fiows. Therefore, models where 
agents can move in continuous space are more likely 
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to be successful. Among them, the soc ial force model 
introduced by Helb i ng and coworkers (iHelbins . Il994 



mtroduced by ±ieib i ng and coworkers lltieiDms . Iiijjj^ 
iHelbing et aU . l200l l2000al : iHelbing and Molnaii UMl 

had a big impact: the main reason is that the actual 
forces between agents are computed, which allows more 
quantitative predictions as compared to CA models. 

The social force model is based on the concept that be- 
havioral changes of individuals are driven by an external 
social force, which affects the motivation of the individual 
and determines its actions. Pedestrians have a particular 
destination and a preferred walking speed. The motion of 
a pedestrian is determined by its tendency to maintain its 
speed and direction of motion and the perturbations due 
to the presence of other pedestrians and physical barriers 
(walls) . 

The general equation of motion for a pedestrian i is 



■^(pers) 



,(soc) 



(43) 



where rrii is the mass of i and its velocity. The ten- 
dency to maintain the desired walking speed vf^"^ is ex- 
pressed by the first force term F,-^'^'^*'' , which reads 



(0) 



(44) 



where ti is the reaction time of i. The social force r|'*°'^\ 
which describes the influence of the environment, is a su- 
perposition of two terms, expressing the interaction with 
other pedestrians and with barriers, respectively. The 
interaction with other pedestrians is described by a re- 
pulsive potential, expressing the need to avoid collisions, 
and by an attractive potential, expressing the tendency to 
come closer to persons/objects that the pedestrian finds 
interesting. A typical choice for the repulsive force F^y*^^' 



between two pedestrians i and j is 



ij 



Ai exp 



n 



I] , 



(45) 



where Ai and Bi are constants, indicating the strength 
and the range of the interaction, respectively, Vij is the 
sum of the radii of the pedestrians, which are modelled as 
disks, dij the distance between their centers of mass and 
Uij the normalized vector pointing from j to i. Pedes- 
trians try to keep a certain distance from walls/ obstacles 
as well, and the repulsion is often described by the same 
term as in Eq. (|45)) . where one considers the distance of 
i from the nearest point of the walls. The noise ^i{t) is 
added to account for non-predictable individual behav- 
ior. 

This simple framework predicts realistic scenarios, like 
the formation of ordered lanes of pedestrians who in- 
tend to walk in the same direction and the alternation 
of streams of pedestrians that try to pass a narrow door 
into opposite directions. The existence of lanes reduces 
the risk of collisions and represents a more efficient config- 
uration for the system. On the other hand, this is a spon- 
taneously emerging property of the system, as the agents 



are not explicitly instructed by the model to behave this 
way. The repulsion between pedestrians moving towards 
each other implies that the pedestrians shift a little aside 
to avoid the collision: in this way small groups of people 
moving in the same direction are formed. These groups 
are stable, due to the minimal interactions between peo- 
ple of each group, and attract other pedestrians who are 
moving in the same direction. 

In (jHelbing et al\ . l2000bD it is shown that a nontrivial 
non-equilibrium phase transition is induced by noise: a 
system of particles driven in opposite directions inside a 
two-dimensional periodic strip can get jammed in crys- 
tallized configurations if the level of noise exceeds some 
critical threshold {freezing by heating), in contrast to the 
expectation that more noise corresponds to more disor- 
der in the system. This can explain how jams can arise in 
situations of great collective excitation, like panic. Sur- 
prisingly, the crystallized state has a higher energy than 
the disordered state corresponding to particles flowing 
along the corridor, so it is metastable. 

The model introduced in ( Helbing et all . l2000b( ) has 
been adapted to simulate situations in which people in- 
side a room are reached by a sudden alarming informa- 
tion (e.g. fire alarm) and t ry to run away throug h one 
of the exits (escape panic) ( Helbing et all l2000a^ . Ad- 
ditional force terms are considered to account for real- 
istic features of panicking crowds, like the impossibility 
of excessive body compression and of tangential motion 
of people about each other and along the walls. The 
model describes phenomena observed in real panic situ- 
ations: for example, people attempting to leave a room 
through a single narrow exit, generate intermittent clog- 
ging of the exit, so that people are unable to flow con- 
tinuously out of the room, but groups of individuals of 
various sizes escape in an irregular succession (Fig. 118b .) . 
This bursty behavior was actually observed in an em- 
pirical stud^;_on_jnice_ate to exit out of a water 
pool (Salo ma et all |2003[ ). Moreover, due to the fric- 
tion of people in contact, the time to empty the room is 
minimal in correspondence to some optimal value of the 
individual speed: for higher speeds, the total escape time 
increases {faster is slower effect). Placing columns near 
the exits improves the situation, even if it seems against 
intuition. Another situation deals with people trying to 
escape from a smoky room, i.e., a room whose exits are 
not visible unless one happens to stand close to them 
fFig. [T8b). In this case, the agents do not have a prefer- 
ential direction of motion, as they have first to find the 
exits. The question is whether it is more effective for the 
individuals to act on their own or to rely on the action of 
the people close to them. The process is modeled by in- 
troducing a panic parameter, that expresses the relative 
importance of independent action and herding behavior 
[where herding is simulated by a term analogous to the 
alignment rule in (|vicsekiLaII,[i99l)]. It turns out that 
the optimal chances of survival are attained when each 
individual adopts a mixed strategy, based both on per- 
sonal initiative and on herding. In fact, through Individ- 
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FIG. 19 Applause dynamics. Time evolution of the sound in- 
tensity of applauses in theaters and concert halls. The rhyth- 
mic applause occurs in the time window between 12s and 25s. 
a) Global noise intensity, measured by a microphone (aver- 
aged over a narrow time window of 0.2s). b) Local noise in- 
tensity, c) Average noise intensity over a moving time window 
of 3s. d) Order parame ter q. e) Clapping period. Adapted 
from l|Neda et ali |2000| ') . 



FIG. 18 Panic behavior. (Top) Escape from a room with a 
single exit. The exit is clogged by the people, who can leave 
the room only from time to time in bunches of individuals. 
(Bottom) Escape from a smoky room. The time to empty the 
room is minimal if people maintain their self-c ontrol and look 
at wh at the others are doing. Adapted from l|Helbing et ali . 
I2000al '). 



ualistic behavior some lucky ones find quickly the exits 
and are followed by the others because of imitation. A 
detailed account of evacuation dynamics can be found in 
the recent review ( Schadschneider et a/.l . l2008l ). 

Other studies concern the statistical features of 
mo bility patterns of ind ividuals in physical space. 
In ( Brockmann et aZl |2006| ) the scaling properties of hu- 
man travels have been investigated by tracking the world- 
wide dispersal of bank notes through bill-tracking web- 
sites. It turns out that the distribution of traveling dis- 
tances decays algebraically, and is well reproduced within 
a two-parameter continuous time random walk model. 
Moreover also in this case a power-law distribution for 
the interevent times between two consecutive sighting of 
the banknotes has been reported. These studies high- 
light the importance of the web as a platform for social 
oriented experiments (see also Sect. IVIlI)) . Very recently 
the mobility patterns of individuals have been investi- 
gated by tra cking the geographica l location of mobile 
phone users ( Gonzalez et al\ . |2008| ). Also in this case 



the distribution of displacements over all users is well 
approximated by a truncated power-law and analyzed in 
terms of truncated Levi flights, i.e. random walks with 
power-law distributed step sizes. 



C. Applause dynamics 

Applause represents another remarkable example of so- 
cial self-organization. At the end of a good performance 
the audience, after an initial uncoordinated phase, often 
produces a synchronized clapping, where everybody claps 
at the same time and with the same frequency. Synchro- 
nization occurs in many biological and sociological pro- 
cesses (Strogatz and Stewart, 1993; Strogatz, 1994), from 
the flashing of Southeast Asian fireflies to the chirping of 
crickets, from oscillating chemical reactions to menstrual 
cycles of women living together for long times. 

Rhythmic applauses have been explored in detail, both 
empirically a nd theoretic a lly. In the first pioneerin g in- 
vestigations (|Neda et all l2000l: iNeda et all l2000t) ap- 
plauses were recorded after several good theater and 
opera performances, with microphones placed both at 
some distance from the audience and close to randomly 
selected spectators. The intensity of noise, heard at some 
distance from the audience or by one of the spectators 
(Fig. [TOb . b) shows that the signal becomes periodic dur- 
ing the phase of the rhythmic applause, and that the av- 
erage intensity of the sound decreases (Fig. [TOb). Spec- 
tators usually start with a high frequency of clapping, 
which reduces to approximately a half in the synchro- 
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nized phase (Fig. Moreover, an order parameter 

for synchronization was introduced and computed, i.e. 
the maximum of the normahzed correlation between the 
signal s(t) and a harmonic function 



max 



/;+^s(t)sin(2Vr- 



It-T 



(46) 



where $ and T vary in [0, 2tt] and between two convenient 
values Tmin and Tmax, respectively. The time evolution 
of the order parameter q{t) clearly indicates the emer- 
gence of synchronization (Fig. [TUH). 

The dynamics of rhythmic applause has bee n explained 
in th e framework of the Kuramoto model ( Kuramotol 
Il975[ ). The audience is described by a system of non- 
identical N rotators. Each rotator i is characterized by 
its phase 0i and natural frequency tOi, which is assigned 
according to a probability distribution g{Lo). The dy- 
namics is expressed by a system of A'^ coupled differential 
equations 



dt 



K 

N 



N 



(47) 



D. Mexican wave 

We conclude with another striking example of coher- 
ent collective motion, i.e., the Mexican wave, also called 
La Ola, which is the wave created by spectators in 
football stadia when they rapidly leap from the seats 
with their arms up and successively sit down while a 
nei ghboring sect i on of people starts the same s equence. 
In (jFarkas all 120021 : iFarkas and Vicsekl . [looi) an em- 
pirical study of this peculiar phenomenon has been re- 
ported and simple models to describe it proposed. These 
mod els were inspired by the literature on excitable me - 
dia (|Bub et all 120021 : iCreenberg and Hastingl lUli), 
where each unit of the system can switch from an in- 
active to an active state if the density of active units in 
their neighborhood exceeds a critical threshold. The in- 
fluence of a neighbor on an excitable subject decreases 
with its distance from the subject and is higher if the 
neighbor sits on the side where the wave comes from. 
The total influence of the neighbors is compared with 
the activation threshold of the spectator, which is uni- 
formly distributed in some range of values. It turns out 
that a group of spectators must exceed a critical mass 
in order to initiate the process. The models are able to 
reproduce size, form, velocity and stability of real waves. 



where K is the interaction coupling. Each rotator in- 
teracts with all its peers, trying to minimize the phase 
difference with the other rotators. In the thermodynamic 
limit N oo there is a critical coupling Kc such that, 
ii K > Kc, there is a synchronized phase in which the 
rotators are partially aligned. If the distribution of natu- 
ral frequencies g{i-tj) is Gaussian, with dispersion D, one 
finds that = i.e. the critical coupling is 

proportional to the dispersion of the natural frequencies. 
So, if D is small, synchronization is likely to occur. The 
clapping frequency of the rhythmic applause is indeed 
small, and so is its dispersion, as confirme d by experi- 
ments performed on individual spectators ( Neda et al\ . 
[2OOQ ) . On the other hand, the frequencies of the enthu- 
siastic clapping at the beginning of the applause have a 
much higher dispersion, which hinders synchronization. 

From Fig. [19] we see that the phase of rhythmic ap- 
plause is not stable, as it disappears after a while, al- 
though it may be reached again at a later time. This is 
due to the fact that the noise intensity decreases during 
the rhythmic applause, because of the lower clapping fre- 
quency (Fig. [int). Spectators may find that the reduced 
intensity of the sound does not adequately represent their 
appreciation of the performance and decide to increase 
the clapping frequency in order to produce more noise. 
However, this brings the system again in the phase of 
enthusiastic clapping, and synchronization is lost. 

In a more realistic model, spectators are represented as 
two-mode stochastic oscillators, and are only driven by 
the goal of producing some desir ed global level of n oise, 
with or without synchronization ( Neda et aLl . l2003l ). 



VII. FORMATION OF HIERARCHIES 

Hierarchical organization is a peculiar feature of many 
animal species, from inse cts to f ishes, from b i rds to mam- 
mals, includin g hum ans (|Alled . Il942t IChasd . Il980l : iGuhll . 
119681 : IWilso^ . Il97lh . Individuals usually have a well de- 
fined rank inside their group, and the rank essentially 
determines their role in the community. Highly-ranked 
individuals have easier access to resources, they have bet- 
ter chances to reproduce, etc.. Hierarchies also allow for 
an efficient distribution of different tasks within a society, 
leading to a specialization of the individuals. 

The origin of hierarchical structures in animal and hu- 
man societies is still an open issue and has stimulated 
a lot of activity in the past decades. The problem is 
to understand why and how from individuals with ini- 
tial identical status, inequalities emerge. For instance, 
one wonders how, in human societies, a strongly elitar- 
ian wealth distribution could arise starting from a society 
where people initially own an equal share of resources. A 
possible explanation is that hierarchies are produced by 
intrinsic attributes of the individuals, e.g. differences in 
weight or size (for animals), and ta lent or charisma (for 
humans). However, already in 1951 ( Landaul . ll951al lbl). it 
was pointed out that intrinsic factors alone could not be 
responsible for the hierarchies observed in animal com- 
munities, and that the interactions between individuals 
play a crucial role in the establishment of dominance re- 
lationships. The hypothesis that hierarchy formation is a 
self-organization phenomenon due to social dynarn i cs has 
meanwhile become the most widespread ( Chasel Il982t 
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IChase et all l2002t iFrancid . Il988f) . Here we discuss the 
most important results in this area. 



A. The Bonabeau model 



Dominance relationships seem to be determined by 
the outcome of fights between individuals. Laboratory 
experiments on various species hi nt at the ex i stence 
of a positive feedback me c hanism (IChase et al. ^ 1994t 



iHogeweg and Hespeil . Il983l : iTheraulaz et al\ . Il995 ). ac- 
cording to which individuals who won more fights have 
an enhanced probability to win future fights as compared 
to those who were less successful (winner/loser effects). 
Memory effects are also important: observations on cock- 
roaches show that insects removed from a group and rein- 
serted after s ome time do not rega in immediately their 
original rank ( Dugatkin et aLl . ll994f) . Based on these em- 
pirical findings, Bonabeau et al. proposed a simple model 
to explain the eme rgence of hierarchies f rom an initial 
egalitarian society ( Bonabeau et aLl . ll995[ ). 

In the Bonabeau model, agents occupy the sites of a 
two-dimensional square lattice with linear dimension L. 
Each site can host only one agent and the density of the 
agents on the lattice is p, which is the control parameter 
of the system. Every agent performs a random walk on 
the lattice, moving to a randomly selected neighboring 
site at each iteration. If the site is free, the agent occupies 
it. If the site is hosting another agent, a fight arouses 
between the two, and the winner gets the right to occupy 
the site. In this way, if the winner is the attacking agent, 
the two competitors switch their positions at the end of 
the fight, otherwise they keep their original positions. 

The outcome of the fight depends on the relative 
strength of the two opponents. The strength F of an 
agent grows with the number of fights it wins. Agent i is 
stronger than agent j if Ft > Fj . The fight is a stochastic 
process, in which the stronger agent has better chances 
to prevail, but it is not bound to win. The probability 
Qij that agent i defeats agent j is expressed by a Fermi 
function: 



(48) 



l + exp{-77[F,-^,]}' 



where is a sort of inverse temperature, measuring the 
level of randomness in the dynamics: for r] oo, the 
stronger agent always wins, for — s- both agents win 
with the same probability 1/2. When an agent wins/loses 
a fight, its strength is increased/decreased by one unit. 
Memory effects are expressed by a relaxation process of 
the strengths of all agents, which decrease by a constant 
fraction fi at each time step. 

In an egalitarian society, all agents have equal strength. 
A broad distribution of strength would indicate the ex- 
istence of hierarchies in the system and is reflected in 
the distribution of the winning probabilities Qij in the 
stationary state. So, the variance 



can be used as order pa rameter for the sys- 
tem ( Sousa and Stauffe^ . 120001 ). For an egalitarian so- 
ciety, (7 = 0; a hierarchical society is characterized by a 
strictly positive value of the variance a. 

In simulations of the Bonabeau model, agents are ini- 
tially distributed at random on the lattice, the strengths 
of all agents are usually initialized to zero (egalitarian 
society) and one iteration consists of one sweep over all 
agents, with each agent performing a diffusion/fighting 
step. The main result is that there is a critical density 
Pc(?7, m) for the agents such that, for p < Pc{r], /i), society 
is egalitarian, where as for p > Priv, iS ] a hie rarchical or- 
ganization is created ( Bonabeau et aLl . ll995f ). This is due 
to the interplay between competition and relaxation: at 
low densities, fights are rare, and the dynamics is domi- 
nated by relaxation, which keeps the individual strengths 
around zero; at large densities, instead, the grow of the 
strengths of some individuals is not counterbalanced by 
relaxation, and social differences emerge. 



1. Mean field solution 

The model can be analytically solved in the mean field 
limit, in which spatial correlations are absent and the 
density p expresses the probability for two agents to meet 
and fight. The evolution equation for the strength of an 
agent Fi is 



dFi 



N 



sinh?7(Fi — Fj) 
+ coshrj{Fi — Fj 



(50) 

where N is the number of agents. From Eq. ([50]) one 
derives that the mean strength {F) decays exponentially 
to zero, so the stationary states are all characterized by 
{F) = 0. To check for the stability of the special solution 
Fi = 0,Vi, corresponding to an egalitarian society, one 
computes the eigenvalues of the Jacobian matrix 



ms,,-P^. (51) 



(49) 



It turns out that there are only two different eigenval- 
ues, namely —p < and prj/2 — p (in the large- TV limit), 
which is negative only when p < pc{r],p) = 2p/rj. We 
conclude that, for p > Pcij], p), the egalitarian solution is 
linearly unstable, and a hierarchy emerges. For 77 < 2 the 
hierarchy transition is discontinuous, because a subcriti- 
cal bifurcation occurs at pc{ri,p). In this case there are 
several metastable states, along with the egalitarian one, 
so it is possible to observe hierarchical structure even for 
p < pc{ri,p). Indeed, simulations reveal that the sta- 
tionary stat e is very sensitive to th e choice of the initial 
conditions ([Bonabeau et all . Il995h . For 77 > 2, the bi- 
furcation is critical, so the transition is continuous and 
characterized by critical exponents and critical slowing 
down. 
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2. Modifications of the Bonabeau model 

The most popu lar modification of the Bonabeau model 
( Staufferl . l2003lJ ) consists in replacing 77 in Eq. ([48]) with 
the standard deviation a of Eq. By doing so, the 

probabilities {Qij} are calculated using the variance of 
their distribution, which changes in time, so there is 
a feedback mechanism between the running hierarchical 
structure of society and the dominance relationships be- 
tween agents. 

For this model, analytical work in the mean field 
limit ( Lacasa and Luqug . [loO^) showed that the egali- 
tarian fixed point is stable at all densities, at odds with 
simulation results, that support the existence of a phas e 
transition to a hierarchical system ( Malarz et aZ.l . l200^ . 
The apparent discrepancy is due to the fact that, above 
a critical density, a saddle-node bifurcation takes place. 
Both the egalitarian and the hierarchical fixed points 
are stable, and represent possible endpoints of the dy- 
namics, depending on the initial condi tions. The mode l 
has been simulated on regular lattices ( Stauffeil . l2003b( ). 
compl ete graphs (M alarz et at, 200^ and scale-free net- 
works ( Gallosl . l2005( ). The phase transition holds in every 
case, although on scale-free networks the critical density 
may tend to zero in the thermodynamic limit of infinite 
agents. On the lattice, the model yields a society equally 
divided into leaders and followers, which is not realistic. 
If the variation of the strength is larger for a losing agent 
than for a winning agent, instead, the fraction of agents 
that turn into lead ers decreases rapidly with the a mount 
of the asymmetry ( Stauffer and Sa Martinsl [2003^ . 



Some authors proposed modific ations of the moving 
rule for the agents. In (|Odagaki and Tsuiiguchil . l2006t 



iTsuiiguchi and Odagakil |2007[ ) two particular situations 
have been investigated, corresponding to what is called 
a timid and a challenging society, respectively. Timid 
agents do not look for fights, but try to move to a free 
neighboring site. If there are none, they pick a fight with 
the weakest neighbor. Two phase transitions were ob- 
served by increasing the population density: a continu- 
ous one, corresponding to the emergence of a middle class 
of agents, who are fairly successful, and a discontinuous 
one, corresponding to the birth of a class of winners, who 
win most of their fights. In a challenging society, agents 
look for fights, and choose the strongest neighbor as op- 
ponent. Hierarchies already emerge at low values of the 
population density; in addition, since strong agents have 
a big attractiveness, spatial correlations arise with the 
formation of small domains of agents at low and inter- 
mediate densities. 



B. The Advancement- Decline model 

A simple model based on the interplay between ad- 
vancement and decline, s imilar to the Bonabeau mode l, 
has been proposed in ( Ben-Naim and Redneii l2005t ). 
Agents have an integer-valued strength, and interact 



pairwise. The advancement dynamics is deterministic: 
the strength of the stronger competitor increases by one 
unit. If both agents have equal strength, both advance. 
The memory effect of the Bonabeau model now consists 
of a declining process, in that the strength of each indi- 
vidual decreases by one unit at rate r, as long as it is 
positive. The parameter r fixes the balance of advance- 
ment and decline. 

The model is solvable in the mean field limit. We call 
fj{t) the fraction of agents with strength j at time t, and 

Fk — j=o /j • '^^^ dynamics of the cumulative distri- 
bution Fk is described by the non-linear master equation 



dFk 
dt 



r(Ffe+i-Ffe)-f (52) 



where the first two terms express the contribution of the 
decline process, the other two that of advancement. At 
the beginning all agents have zero strength, i.e. /o(0) = 
1, corresponding to i^fc(O) = l,Vfc > (the reflecting wall 
at fc = conies from imposing F-i{t) = 0,Vi). 

In the limit of continuous strength, Eq. (|52p becomes 
a partial differential equation. 



dF{k,t) 
di 



= [r-F{k,t)] 



dF{k,t) 
dk ■ 



(53) 



In the limit of t — > 00 and finite strength, and since 
F(fc,oo) < 1 by definition, Eq. (|53p leaves only two pos- 
sibilities: 



F{k, 00) 
dF{k, 00) /dk 



r, 
0, 



for r < 1; 
for r > 1. 



(54) 



In this way, for r > 1, all agents have finite strength 
in the stationary state. For r < 1, instead, we have a 
hierarchical society, with a fraction r of agents having 
finite strength (lower class), coexisting with a fraction 
1 — r of agents whose strength diverges with t (middle 
class). 

A dimensional analysis of Eq. reveals that the 

strength of the middle class agents increases linearly 
with t. One can then check for solutions of the form 
Fk ~ ^{x), with X — k/t and the boundary condition 
<i>(oo) = 1. Eq. ([53]) becomes 



x^'{x) = mx) ~r)<^>'{x), 



whose solution is 



$(a;) 



X - 
1, 



a; < 1 
x> 1 



(55) 



(56) 



From Eq. (j56p one deduces that the maximum strength 
of the system is (1 — r)t, i.e. it grows linearly with time, 
as expected. In the limit of finite strength {x —^ 0) one 
has $(0) — r, which recovers the result on the fraction 

of lower class agents^ 

In a later paper (jBen-Naim et al ]. l2006bD . the model 
has been generalized by introducing a stochastic ad- 
vancement dynamics, in which the stronger competitor 
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of a pair of interacting agents wins with a probabil- 
ity p. This model yields a richer phase diagram. In 
some region of the parameter space, a new egalitarian 
class emerges, in which the strength distribution of the 
agents is strongly peaked and moves with constant veloc- 
ity, like a traveling wave. The model has been success- 
fully apphcd to describe the dynamics of sport competi- 
tions (Ben-Naim et ai, 2007). Moreover, it has inspired a 
generalization to competitive games involving more th an 
two players at the same time ()Ben-Naim et aLl . l2006al) . 



VIII. HUMAN DYNAMICS 
A. Empirical data 

One of the key questions in social dynamics concerns 
the behavior of single individuals, namely how an indi- 
vidual chooses a convention, takes a decision, schedules 
his tasks and more generally decides to perform a given 
action. Most of these questions are obviously very diffi- 
cult to address, due to the psychological and social fac- 
tors involved. Nevertheless in the last few years several 
studies have tried to quantitatively address these ques- 
tions, mainly relying on the availability of data through 
the web. A first valuable source of data have been the 
logs of email exchanges. In particular the stru cture of 
emai l networks has been first studied in (Ebel et all . 
I2OO2I : iNewman et'd] . l2002f ). focusing on the spreading 
of informatic viruses. The emergence of coherent, self- 
org anized, structures in email traffic has been reported 
m (jEckmann et all\2004) . using an information-theoretic 
approach based on the analysis of synchronization among 
triplets of users. It has been highlighted how nontrivial 
dynamic structures emerge as a consequence of time cor- 
relations when users act in a synchronized manner. The 
observed probability distribution of the response time t 
until a message is answered, features a broad distribu- 
tion, roughly approximated by a 1/t power law behavior 

(Fig. EOl) . T he sam e kind of data have been analyzed 

in ((Johansenl . l2004t) and a generalized response time dis- 
tribution ~ 1/t for human behavior in the absence of 
deadlines has been suggested. 

A different scaling, as far as the distribution of re- 
sponse times is concerned , has been observed in mail 
correspondence patt erns ( Oliveira and Barabasil 120051 : 
I Vazquez et aLl . [2006l ). In the correspondence of Einstein, 
Darwin and Freud it has been found that the probability 
that a letter will be replied in r days is well approximated 
by a power-law t~", with a — 3/2. 

From all these empirical data a picture emerges where 
human dynamics is characterized by bursts of events sep- 
arated by short interevent time s followed by long periods 
of inactivity. This burstiness ( KleinberS . mm can be 



10^ 



quantified by looking for instance at the distribution of 
interevent times that typically displays heavy-tailed non - 
Poissonian statistics. A recent study (l Alfi et all . l2007t ) 
investigated how people react to deadlines for conference 
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FIG. 20 Distribution of the response time until an email mes- 
sage is answered. (Inset) The same distribution is measured 
in ticks, i.e., units of messages sent in the system. Binning 
is logarithmic. The solid lines follow At~ ^ and are meant as 
guides for the eye. From l|Eckmann et all |2004| ) 



registration. 



B. Queuing models 

Let us now focus on the explanation of the above men- 
tioned pheno menology. The very sam e database collected 
an d used in j Eckmann et al\ . \2004) has been analyzed 
m (|Barabasil [2005l), where the origin of bursts and heavy 
tails in the probability distribution of the response time 
to an email has been explained as a consequence of a 
decision-based queuing process. The model is defined as 
follows. Each human agent has a list with L tasks, each 
task being assigned with an a priori priority parameter 
Xi (for i G [1, chosen from a distribution p{x). At 

each time step the agent selects the task with the high- 
est priority with probability p and executes it, while with 
probability 1 — p a randomly selected task is performed. 
The executed task is then removed from the list and re- 
placed with another one with priority again randomly 
extracted from p{x). Computer simulations of the model 
showed that for the deterministic protocol {p ^ 1) the 
probability distribution of the times spent by the tasks 
on the list features a power-law tail with exponent a — 1. 
The exa ct solutio n of th e Barabasi model for two tasks 
2) (|Vazaued . |2005() confirmed the 1/t behavior 



with an exponential cut-off over a characteristic time 
(In (2/(1 +p)))~^, which diverges for p ^ 1. The limit 
p — > 1 is subtle because not only the characteristic time 
diverges but the prefactor vanishes. 

This issue has been addressed 

in ( Gabrielli and Caldarellil . 120071 ). where an exact 
probabilistic description of the Barabasi model for L = 2 
has been given in the extremal limit, p = 1, i.e., when 
only the most recently added task is executed. In this 
limit the Barabasi model ca n be exactly mapped into 
Invasion Percolation (IP) (jWilkinson and Willemsenl . 
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I1983D i n d = 1 . Usin g the so-called run time statistics 
(RTS) ([Marsilil Il994[ ). originally introduced to study 
the IP problem in d = 2, it has been found that the 
exact waiting time distribution for a task scales a s 
r~^, unhke the result l/r found in ( Vazquej . l2005l ). 



which is valid for the stationary state when < p < 1. 
This behavior disappears in the limit p — *■ 1, since the 
prefactor vanishes. In summary, for < p < 1 finite 
time deviations from stationarity relax exponentially 
fast and the dynamics is well described by the stationary 
state. On the other hand, for p — > 1 the stationary state 
becomes trivial and finite time deviations relax so slowly 
that the task list dynamics has to be described as an 
intrinsically nonstationary dynamics. 

In (I Vazquez et the case where limitations on 

the number of tasks an indi vidual can handle at any time 
is discussed (ICobhaml . [T95l . Here tasks are, for instance, 
letters to be replied. One assumes that letters arrive at 
rate A following a Poisson process with exponential ar- 
rival time distribution, while responses are written at rate 
fi. Each letter is assigned a discrete priority parameter 
X — 1,2, ...,r upon arrival and the unanswered letter with 
highest priority will be always chosen for a reply. This 
model predicts a power law behavior for the waiting time 
distribution of the individual tasks, with an exponent 
equal to 3/2, with an exponential cut-off characterized 
by a characteristic time tq given hy tq = 1 / {fi{l — ^/pY), 
where p = X/ fJ.. The ratio p is called traffic intensity 
and it represents the natural control parameter for the 
model. For p = 1 one has a critical regime where the 
queue length performs a 1 — d random walk with a lower 
bound at zero. This allows to explain the origin of the 
exponent 3/2 in terms of the exponent of the return time 
probability for a random walker. 

Conditions for the emergence of scaling in the inter- 
event time distribution have been addressed in (Hidalgci, 
hood ). A generalization of the queuing model introducing 
a time-dependent priority for the tasks in the list (aging) 
has been discussed in (Blanchard and Honglcr, 2006). A 
further generalization of the queue model with contin- 
uo us valued priorities has been introduced and studied 
in ( Grinstein and Linskeii [2006) . In this case a new task 
is added to the queue with probability A and assigned 
a priority a; (0 < a; < 1) chosen from a probability dis- 
tribution p{x). With probability ^ the highest priority 
task is executed. Two asymptotic waiting time distri- 
butions have been found analytically. For A = /i < 1 
one obtains a power law (with exponent 3/2) for the 
waiting time distribution. For X < p < 1 a character- 
istic time To = (^ — A)^/(4p(l — A)) > emerges such 
that, for r ^ To, the waiting time distribution goes as 
P(t) - e-^/^OT-^/^, while for 1< T < To P(t) - t'^/^. 
The asymptotic behaviour of P{t) changes with the 
introduction of a cost for switching between different 
classes of tasks ( Grinstein and Linsked . [2006 ). 

Finally, in ( Bedogne and Rodgersl . |2007[ ) a continu- 
ous model has been introduced, which adopts an infi- 
nite queuing list where each element is assigned a real 



non-negative number representing its priority of execu- 
tion. The priorities are drawn from a given probability 
distribution function p{x). At each time step a task is 
selected with a selection probability function, which is an 
increasing function of the priority: 



U{x,t) 



yp{x,t) 



x^P{x,t)dx' 



(57) 



where 7 is a non-negative real number and P{x,t) is the 
so-called waiting room priority distribution. The result- 
ing waiting time distribution P{t) can be analytically 
found and it depends explicitly on the priority distribu- 
tion density function p{x). The scaling P{t) ~ 1/t is 
recovered when p[x) is exponentially distributed. 



C. The social web 
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Wide 
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Recently, a 
dynamics plays 
quickly gaining 

Web: Collaborative T aggin g (jCattuto et al\ . 
iGolder and Huberma 3, [200i). In web 

cations like del.icio.us (http://del.icio.us), 
Flickr (|http : //www . f li ckr . com^ Cite ULike 

(http://www.citeulike.org), users manage, share 
and browse collections of online resources by enriching 
them with semantically meaningful information in the 
form of freely chosen text labels {tags). The paradigm 
of collaborative tagging has been successfully deployed 
in web applications designed to organize and share 
diverse online resources such as bookmarks, digital 
photographs, academic papers, music and more. Web 
users interact with a collaborative tagging system 
by posting content {resources) into the system, and 
associating text strings {tags) with that content. Fig.l^ 
reports the frequency-rank distributions for the tags 
co-occurring with a few selected ones. The high-rank 
tail of the experimental curves displays a power-law be- 
havior, signature of an emergent hierarchic al structure , 
corresponding to a generalized Zipf's law (jZipj . Il949( ) 
with an exponent between 1 and 2. Since power laws 
are the standard signatur e of se l f-organization a nd of 
human a c tivity ( Barabasil 120051 : iMitzenmacheil . 120031 : 
iNewmanl . I20Q5D . the presence of a power-law tail is 
not surprising. The observed value of the exponent, 
however, deserves further investigation, because the 
mechanisms usually invoked to explain Zipf's la w and 
its generalizations (jZanette and Montemurrol . l2005i ) 
do not look very realistic for the case at hand, and a 
mechanism grounded on experimental data should be 
sought. Moreover, the low-rank part of the frequency- 
rank curves exhibits a flattening typically not observed 
in systems strictly obeying Zipf's law. Several aspects 
of the underlying complex dynamics may be responsible 
for this feature: on the one hand this behavior points 
to the existence of semantically equivalent and possibly 
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FIG. 21 Frequency-rank plots for tags co-occurring with a 
selected tag. Experimental data (black symbols) are shown 
for del.icio.us (circles for tags co-occurring with the popu- 
lar tag blog, squares for ajax and triangles for xml). For the 
sake of clarity, the curves for ajax and xml are shifted down 
by one and two decades, respectively. All curves exhibit a 
power-law decay for high ranks (a dashed line corresponding 
to the power law R~^^* is provided as an aid for eye) and 
a shallower behavior for low ranks. Gray (red online) sym- 
bols are theoretical data obtained by computer simulations 
of th e stochastic process described in (Cattuto et al., 2007, 
|2006| V (Inset) th e same graph for the much younger system 
Connotea. From l|Cattuto et aZ.1. 120071 '). 



competing high-frequency tags (e.g. hlog and blogs). 
More importantly, this flattening behavior may be as- 
cribed to an underlying hierarchical organization of tags 
co-occurring with the one we single out: more general 
tags (semantically speaking) will tend to co-occur with 
a larger number of other tags. 

At the global level the set of tags, though determined 
with no explicit coordination, evolves in time and leads 
towards patterns of terminology usage. Hence one ob- 
serves the emergence of a loose categorization system 
that can be effectively used to navigate through a large 
and heterogeneous body of resources. It is interesting to 
investigate the way in which users interact with those sys- 
tems. Also for this system a hyperbolic l aw for the user 
access to the system has been observed ( Cattuto et all 
[2007). In particular if one looks at the temporal auto- 
correlation function for the sequence of tags co-occurring 
with a given tag (e.g. blog), one observes a l/(i -I- r) be- 
havior, which suggests a heavy-tailed access to the past 
state of the system, i.e., a memory kernel for the user 
access to the s ystem. On this basis a stochastic model of 
user behavior ( Cattuto et 'all. I2OO7D has been proposed, 
embodying two main aspects of collaborative tagging: (i) 
a frequency-bias mechanism related to the idea that users 
are exposed to each other's tagging activity; (ii) a notion 
of memory (or aging of resources) in the form of a heavy- 
tailed access to the past state of the system. Remarkably, 
this simple scheme is able to account quantitatively for 
the observed experimental features, with a surprisingly 



high accuracy (see Fig. [21]) . This points to the direc- 
tion of a universal behavior of users, who, despite the 
complexity of their own cognitive processes and the un- 
coordinated and selfish nature of their tagging activity, 
appear to follow simple activity patterns. 

The dynamics of information access on the web rep- 
resents another source of data and several experiments 
have b een performed in the las t few y ears. In (JohanscJ, 
I2QOII : iJohansen and Sornettd . I2OOOI ) the dynamic re- 
sponse of the internauts to a point-like perturbation as 
the announcement of a web intervie w on stock market 
crash e s has been investigated. In ( Chessa and Mur"r3 . 
l2004l . |2006^ a cognitive model, based on the math- 
ematical theory of point processes, has been pro- 
posed, which extends the r esults of ('J ohanseiil . I2OOII : 



iJohansen and Sornette 



namics. In ( Dezso et 



2OOOI) t o download relaxation dy- 



2000) the visitation patterns of 



news documents on a web portal have been considered. 
The proble m of collective attent i on ha s been recently ad- 
dressed in ( Wu and Hubermanl . l2007i ) in the framework 
of the website Digg.com. Digg allows users to submit 
news stories they discover from the Internet. A new sub- 
mission immediately appears on a repository webpage 
called Upcoming Stories, where other members can find 
the story and, if they hke it, add a digg to it. A so-called 
digg number is shown next to each story's headline, which 
simply counts how many users have digged the story in 
the past. The distribution of diggs accumulated by the 
news stories has been found to be well approximated by a 
log-normal curve. This behaviour has been explained in 
terms of a simple stochastic model where the number of 
diggs acquired by a news story at time t can be expressed 
as: 



Nt = {l + rtXt)Nt^ 



(58) 



where Xi, X2, Xt are positive i.i.d. random vari- 
ables with a given mean and variance. The parameter rt 
is a time-dependent novelty factor consisting of a series 
of decreasing positive numbers with the property that 
ri = 1 and rt — > as i — s- cx). Experimentally it has 
been found that rt follows a stretched-exponential rela- 
tion with a characteristic time of about one hour. 



IX. SOCIAL SPREADING PHENOMENA 

Opinion dynamics deals with the competition between 
different possible responses to the same political ques- 
tion/issue. A key feature is that the alternatives have 
the same or at least comparable levels of plausibility, so 
that in the interaction between two agents each of them 
can in principle influence the other. In phenomena like 
the propagation of rumors or news, the interaction is in- 
trinsically asymmetric: possible states are very different 
in nature. The flow is only from those who know to 
those who do not. The propagation of rumors or news 



46 



is an instance of the vast class of social spreading phe- 
nomena, that includes the diffusion of fads, the adop- 
tion of technological innovations, the success of consumer 
products mediated by word-of-mouth. The investigation 
ot such types of dyna mics has a long tradition in so- 
ciology and econorn ics (|Bassl Il969t iBikhchandani et al\ . 
ll992i : lRoger^ . l2003D . Work along this line has been per- 
formed recentl y by statistical physicists on the d i ffusio n 
of innovations ( Guardiola et "oLI. I2OO2I: lllas et a/.l . [200l . 
the occurrence of information cascad e s in social and eco- 
nomic systems ( Centola et all 12007: IWa tts. 200. 2), dis - 
aster spreading in infrastructures ( Buzna et all I2OO6I) . 



the emergence o f "hit s" in markets (jHidalgo et aLl . 120061: 
ISinha and Paul l2006t ). Many of the models introduced 
for these phenomena assume that a local threshold in 
the fraction of active neighbors must be overcome for 
the spreading process to occur. Here we will, somewhat 
arbitrarily, review only the activity on the problem of 
rumor spreading. This phenomenon is modeled without 
a threshold: hearing a rumor from a single neighbor is 
generally enough for a single agent to start spreading 
it. It is then clear that rumor spreading bears a lot of 
resemblance with the evolution of an epidemics, with in- 
formed people playing the rol e of infected agents and un- 
informed of susc eptible ones (j Goffman and Newili 119641 : 
iRapoportl . Il953 ^. Obviously there are crucial qualitative 
differences: rumor/news spreading is intentional; it usu- 
ally involves an (at least perceived) advantage for the 
receiver, etc. However most of such differences lie in the 
interpretation of parameters; the analogy is strong and 
the field is usually seen as closer to epidemiology than to 
opinion dynamics. 

When considering rumor spreading some of the rele- 
vant questions to address are similar to those for epi- 
demiology: How many people will eventually be reached 
by the news? Is there an 'epidemic threshold' for the rate 
of spreading, separating a regime where a finite fraction 
of people will be informed from one with the info remain- 
ing confined to a small neighborhood? What is the de- 
tailed temporal evolution? Other issues, more connected 
to technological applications, deal with the cost of the 
spreading process and its efficiency. Rumor dynamics 
has also appealing connections with the search for robust 
scalable communication protocols in large distributed 



systems (jKermarrec et 

"viral" strategies in marketing (jLeskovec et al 



O rotocols m large c Listrip uted 
2003tlVogels 6^0/.!. I2003D a nd 



2QM)- 



Detailed applications of the common models for epi- 
demics to the investigation of empirical data o n the 
dissemination o f ideas exist ( Bettencourt et al\ . 120061 : 
iGoffmanl . I1966D . but the mo st popular model for ru - 
mor spreading, introduced in ( Daley and Kendall . Il964 ) 
(DK), has an important differen ce. As in the SIR model 
for epidemiology (Anderson a nd Mavl . Il991i) , agents are 
divided in three classes: ignorants, spreaders and stiflers, 
i.e., those who have lost interest in diffusing the informa- 
tion/rumor. Their role is exactly the same of the suscep- 
tible, infected, recovered agents of the SIR model. The 
only difference is that while for an epidemics infected (I) 



people become recovered or removed (R) spontaneously 
with a certain rate, typically people stop propagating a 
rumor when they realize that those they want to inform 
are already informed. Hence the transition to state R is 
proportional to the density of spreaders s{t) in the SIR 
model, while it is proportional to s{t)[s{t) + r{t)] in the 
DK model, where r(t) is the density of stiflers. 

The DK model has been studied analytically in the 
case of homogeneous mixing, revealing that there is no 
threshold: for any rate A of the spreadin g process a finit e 
fraction Too of people would be informed ( Sudburvl . ll985( ). 
given by the solution of 



1 



-(1-|-A/Q)r„ 



(59) 



where a is the proportionality constant of the transition 
rate to state R. The same r esult holds for the very simi - 
lar Maki-Thompson model ( Maki and Thompsonl . Tl973l ). 



Hence the nonlinear transition rate removes the thresh- 
old of the SIR model. Clearly when both mechanisms 
for the damping of the propagation are present (self- 
recovery and the nonlinear DK mechanism) a thresh- 
ol d is recovered, since the linear term prevails for small 
s (iNekovee et adl2007l ). 

In the context of statistical physics the focus has 
been o n the behavior o f the DK mode l on co mplex net- 
works ( Liu et all . 120031 : iMoreno et al. I. I2004al lbl). When 
scale- free networks are considered, the fraction Too of peo- 
ple reached decreases compared to homogeneous nets. 
This occurs because hubs tend to become stifiers soon 
and hence hamper the propagation. However, if one con- 
siders the efficiency E of the spreading process, defined 
as the ratio between Too and the total traffic L gener- 
ated, it is found that for any value of the parameters 
scale-free networks are more efficient than homogeneous 
ones, and in a broad range of parameters they are more 
efficient than the trivial broadcast spreading mechanism 
(i.e., each node transmits the message to all its neigh- 
bors). 

A remarkable phenomenon occurs when the DK dy- 
namics takes p lace on the small-w orld Watts-Strogatz 
(WS) network (IZanettd. I2OO1I . l2002t) . In this case there 
is an 'epidemic' transition depending on the rewiring pa- 
rameter p. For p > pc the rumor propagates to a finite 
fraction Too of the network sites. For p < pc instead the 
rumor remains localized around its origin, so r^c vanishes 
in the thermodynamic limit. Notice that pc is finite for 

^ 00, at odds with the geometric threshold charac- 
terizing the small-world properties of WS networks, that 
vanishes in the limit of infinite network. Hence the tran- 
sition is dynamic in nature and cannot be ascribed to a 
pure geometric effect. 

The diffusion of corruption has also been modeled as an 
epidemic-like process, with people accepting or practicing 
a corrupt behavior corresponding to infected individuals. 
The main difference with respect to usual epidemiologi- 
cal models is that the chance of an individual to become 
corrupt is a strongly nonlinear function of the number 
of corrupt neighbors. Other modifications include global 
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coupling terms, modeling the process of people getting 
corrupt because of a perceived high prevalence of cor- 
ruption in the society and the response of the society as 
a whole, which is proportional to the fraction of non- 
corr upt people. The resulting phenomenology is quite 
rich (|Blanchard et aLl . l2005l ). 

It is also worth mentioning the mod el of Do d ds an d 
Watts for social contagion (D odds and Wattsl . l2004t ). 
that introduces memory in the basic SIR model thus 
taking into account the effect of repeated exposure to 
a rumor. Also i nformation sp r eadin g in a population of 
diffusing agents ( Aghari et aLl . [2006| ) has been studied. 

Finally, some activity has been devoted to the related 
problem of gossip spreading. While rumors are about 
some topic of general interest so that they may poten- 
tially extend to all, gossip is the spreading of a rumor 
about some person and hence it is by definition a local 
phenomenon; it may concern only people close to the sub- 
ject in the social network. If only nearest neighbors of 
the subject can spread, the fraction of them reached by 
the gossip exhibits a minimum as a function of the degree 
k for some empirical and model social networks. Hence 
there is an ideal nu mber of connection s to minimize the 
gossip propagation ( Lind et a/] . l2007al lbl) . 



X. COEVOLUTION OF STATES AND TOPOLOGY 

All models considered in the previous sections are de- 
fined on static substrates: the interaction pattern is 
fixed and only opinions, not connections, are allowed to 
change. The opposite case is often considered in many 
studies of network formation: vertices are endowed with 
quenched attributes and links are formed or removed de- 
pending on such fixed node properties. 

In fact, real systems are mostly in between these two 
extreme cases: both intrinsic properties of nodes (like 
opinions) and connections among them vary in time over 
comparable temporal scales. The interplay of the two 
evolutions is then a natural issue to be investigated. More 
interestingly, in many cases the two evolutions are explic- 
itly coupled: if an agent finds that one of his contacts is 
too different, he tends to severe the connection and look 
for other interaction partners more akin to his own prop- 
erties. 

The investigation of the coevolution of networks 
and states has recently attracted interest in several 
contexts, including self-organ iza tion in Boolean 
networks (iBornholdt and Rohlj. |2000|) . synchro- 
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opments is ( Gross and Blasiui . |2008( ). that also contains 
a useful attempt to classify generic patterns of behavior 
in the field. For opinion and cultural dynamics the 
study of adaptive networks is still at the beginning, but 
it promises to be a very active field in the next years. 



In coevolving systems, three ingredients are to be spec- 
ified. The dynamics taking place on nodes usually be- 
longs to standard classes of opinion or cultural dynamics 
(voter-like, majority-rule, continuous opinion, etc.). Link 
dynamics can be of many types. Examples exist of pre- 
scriptions allowing only link deletion or independent link 
addition or removal. A more common, simple yet quite 
realistic, assumption is that links are rewired: an agent 
unhappy about one connection cuts it and forms a new 
link with another agent. In this way the average con- 
nectivity (fc) is conserved. The third ingredient is the 
relative rate of link or node updates, generally modeled 
by a probability 0: a node is chosen at random; if it 
is equal to one of its neighbors (or all, depending on the 
node dynamics) nothing happens. Otherwise with proba- 
bility (j) it rewires one of its connections; with probability 
1 — it updates its state leaving the topology unchanged. 

If node dynamics leads toward consensus, it is reason- 
able to expect that for small 4> there will be asymptoti- 
cally a single network with all nodes in the same state. 
In the opposite limit (j) ^ 1 instead, rewiring will quickly 
split the network in disconnected components that will 
separately reach a different consensus state. These two 
limit behaviors are hence separated by a topological 
phase transition for some critical value (jjc- Such type 
of tra nsition was first reported in ( Holme and Newmanl . 
I2OO6D . where Potts variables assuming G different values 
are defined on the nodes of a network. G is proportional 
to the number of vertices, so that 7 = G/N is constant. 
At each time step a node and a neighbor are selected and 
with probability 1 — the node picks the opinion of the 
neighbor. With probability instead, the node rewires 
the link to a new vertex chosen randomly among those 
having its same opinion. Dynamics continues up to com- 
plete separation in components, within which there is full 
consensus. For large only rewiring is allowed and the re- 
sulting topological clusters trivially coincide with the sets 
of initial holders of individual opinions. The distribution 
of their sizes is multinomial. For small (j) practically only 
opinion changes are allowed and the final clusters are the 
components of the initial graph. The phase transition oc- 
curring at (j)c is characterized by a power-law distribution 
of cluster sizes, with an exponent about —3.5 (Fig. E^ . 
which differs from the one at the threshold of the gi- 
ant component formation in random graphs (mean field 
percolation class) . Finite size scaling gives additional ex- 
ponents that are universal (do not depend on (fc) and 7) 
and different from percolation. Critical slowing down is 
also observed. 

A similar dy namics, but for a bin ary state variable, is 
investigated in ( Vazquez et aLl . [20081) . A master equation 
approach gives an absorbing phase-transition between an 
active phase for (j) < (f)^ — ((fc) — 2)/((fc) — 1), and a 
frozen phase with zero active links for cj) > <j)c- In finite 
systems, numerical simulations show that in the absorb- 
ing phase the network is composed of two disconnected 
components, one with all spins up and the other with 
all spins down, while the active phase corresponds to all 
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FIG. 22 Coevolution of opinions and networks. Histograms 
of cluster sizes in the consensus state for values of cf> above, 
at and below the critical point in panels (a), (b) and (c), 
respectively. In panel (b) the distribution appears to follow 
a power law for part of its range with exp onent —3.5 ± 0.3 , 
as in dicated by the solid line. From (|Holme and Newmanl 
120061 '). 



nodes belonging to the same giant component. 

Nardini et al. ( Nardini et al\ . [20071 ) investigate a very 
similar model, with only the slight difference that links 
are rewired to nodes selected completely at random, with 
no preference for similar vertices holding the same opin- 
ion. They focus on the approach to consensus in the 
limit of small rewiring rate, finding opposite behaviors if 
the selected node picks the opinion of the neighbor (di- 
rect voter dynamics) or the neighbor is convinced by the 
node (reverse voter) . In the first case the consensus time 
becomes logarithmic in the system size N; in the second 
it diverges exponentially with N . A mean field approach 
clarifies the origin of the difference and also allows to un- 
derstand why the consensus time remains proportional 
to N for the AB-model, independently of the direct or 
reverse update prescription. 

Differently from the models just discussed, in the 
first p apers on opinio n dyn a mics in adaptive net- 
works (|Gil and Zanettd . l2006t IZanette and Gill . I2006D 
links can only be deleted. In the initial configuration, 
binary opinions are distributed randomly on the nodes 
of a fully connected network. If two neighbors disagree, 
one of them is set equal to the other with probability pi 
(voter dynamics). With probability (1 — pi)p2 instead 
they get disconnected. Iterating this dynamics a station- 
ary state is reached depending only on the combination 
q — P\l\p\ -|- (1 — Pi)p2]- For small q the system breaks 
down in two communities of similar size and opposite 
opinions, with a large fraction of internal connections. 
For large q there are two possibilities: a single commu- 
nity with the same opinion or one well connected commu- 
nity with a set of poorly connected smaller communities. 
In correspondence to an intermediate value gc, the total 



density of links exhibits a minimum rc((7c); both and 
(/c vanish for large system size. 

S tarting from t h e usu al Axelrod model (see Sec. lIV.Ap 
in ( Centola et al\ . l2007t ) a further step is added: if the 
overlap w between two nodes is exactly zero, the link is 
removed and one of the agents connects to another ran- 
domly chosen vertex. In this way the transition gc be- 
tween a monocultural state and fragmentation is moved 
to much larger values: coevolution favors consensus. 
At qc both a cultural and a topological transition take 
place: the system becomes separated in cultural groups 
that also form topologically disconnected network sub- 
sets. For even higher values of the variability of the ini- 
tial state q ^ q* > qc another transition occurs, involv- 
ing only the network structure. For q > q* the system 
remains culturally disordered but a giant component is 
formed again. In this regime it is likely that each ver- 
tex is completely different from its neighbors, therefore 
it continuously breaks links and looks (unsuccessfully) 
for new more similar par tners. The transition occurring 
at qc can be explained ( Vazquez et all |2007[ ) in terms 
of the competition between the temporal scales of cul- 
tural and topological evolution. The topological transi- 
tion occurring at q* can be instead seen as the value of q 
such that the temporal scale for reaching a topologically 
stationary state is maximum. Another model of adap- 
tiv e network coupled to vectorial o pinions is introduced 
in ( Grabowski and Kosihskil [2006a[) . 

Network rewiring has also been considered for the dy - 
namics of the Deffuant model ( Kozma and Barra3 . l2008h . 
At each time step with probability 1 — w a. step of the 
usual opinion dynamics is performed with confidence 
bound rf, otherwise one agent breaks one link and recon- 
nects it to a randomly chosen other node. By changing 
w it is then possible to go from pure opinion dynamics 
in a static environment to fast topological evolution in a 
quenched opinion state. Coevolution has opposite effects 
on the two transitions exhibited by the model on static 
ER networks. The confidence bound threshold di, above 
which consensus is found, grows with w. The thresh- 
old c?2 of the transition between a polarized state (for 
d2 < d < di) and a fragmented one with no macroscopic 
domains (for c? < ^2) goes instead to zero: d2{'w > 0) = 0. 
The fragmented state disappears because, even for small 
d, a node can rewire its connections and find other agents 
with whom to reach an agreement. Another coevolv- 



ing ge neralization of Deffuant model is (jStauffer et all 
l2006aD . 

Finally, a Glauber zero-tem perature spin d y namic s has 
been recently considered in ( Benczik et all l2007f) . to- 
gether with a particular link update rule, different from 
all previous cases: at each time step two nodes carry- 
ing equal spins are connected with probability p, while 
nodes with opposite values of the spin variable are con- 
nected with probability 1—p. No topological fragmenta- 
tion occurs in this case. Consensus is always reached on 
finite systems but exceedingly long disordered metastable 
states exist. 
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XI. OUTLOOK 

With this review we made a first attempt to summa- 
rize the many activities in the field of the so-caUed social 
dynamics. Our point of view has been that of reviewing 
what has been done so far in this young but rapidly evolv- 
ing area, placing the main emphasis on the statistical 
physics approach, i.e., on the contributions the physics 
community has been giving to social oriented studies. 

Though it is generally very difficult to isolate the con- 
tribution of a given community to an intrinsically inter- 
disciplinary endeavor, it is nevertheless useful to identify 
the contribution the physics community has been giving 
and the role it could play in the future. In this perspec- 
tive it is clear how the statistical physics role in social 
dynamics has been mainly focused on modeling, either 
by introducing brand new models to capture the main 
features of a given phenomenology or performing detailed 
analysis of already existing models, e.g., focusing for in- 
stance on the existence of phase transitions, universal- 
ity, etc.. An inspiring principle has been provided by 
the quest for simplicity. This has several advantages. It 
allows for discovering underlying universalities, i.e., re- 
alizing that behind the details of the different models 
there could be a level where the mathematical structure 
is similar. This implies, in its turn, the possibility to per- 
form mapping with other known models and exploit the 
background of the already acquired knowledge for those 
models. In addition physicists placed a great emphasis 
on the role of scales (system sizes, timescales, etc.) as 
well as on the topology (i.e., the network of interactions) 
underlying the observed phenomenology. 

Closely related to modeling is the data analysis ac- 
tivity, both considering synthetic data coming from sim- 
ulations and empirical data gathered from observations 
of real systems or collected in the framework of newly 
devised experiments. Data analysis is very important 
not only for the identification of new phenomenologies 
or surprising features, but also for the validation of the 
models against empirical data. In this way a positive 
feedback mechanism could be triggered between the the- 
oretical and the experimental activities in order to make 
the results robust, well understood and concrete. 

Methodologically we can identify several important di- 
rections the research in this area should possibly follow. 
It would be crucial fostering the interactions across disci- 
plines by promoting scientific activities with concrete mu- 
tual exchanges among social scientists, physicists, mathe- 
maticians and computer scientists. This would help both 
in identifying the problems and sharpening the focus, as 
well as in devising the most suitable theoretical concepts 
and tools to approach the research. 

As for the modeling activity it would highly desirable 
to identify general classes of behavior, not based on mi- 
croscopic definitions, but rather on large-scale universal 
characteristics, in order to converge to a shared theoreti- 
cal framework based on few fundamental paradigms. The 
identification of which phenomena are actually described 



by the theoretical models must become a priority. For 
instance, the celebrated Axelrod model has not yet been 
shown to describe at least semi-quantitatively any con- 
crete situation. Without applications the intense activity 
on modeling risks to be only a conceptual exercise. 

In this perspective a crucial factor will be most likely 
represented by the availability of large sets of empiri- 
cal quantitative data. The research carried out so far 
only rarely relied on empirical datasets, often insufficient 
to discriminate among different modeling schemes. The 
joint interdisciplinary activity should then include sys- 
tematic campaigns of data gathering as well as the de- 
vising of new experimental setups for a continuous moni- 
toring of social activities. From this point of view the web 
may be of great help, both as a platform to perform con- 
trolled online social experiments, and as a repository of 
empirical data on large-scale phenomena, like elections 
and consumer behavior. It is only in this way that a 
virtuous circle involving data collection, data analysis, 
modeling and predictions could be triggered, giving rise 
to an ever more rigorous and focused approach to so- 
cially motivated problems. A successful example in this 
perspective is the study of traffic and pedestrian behav- 
iors, that in the last few years has attained a high level 
of maturity, leading to rel i able q uantitative predictions 
and control ([Helbing et all |2007^ (see also Sec.lVll). 

We conclude this review by highlighting a few interest- 
ing directions that could possibly have a boosting eflfect 
on the research in the area of social dynamics. 



1. Information dynamics and the Social Web 

Though only a few years old, the growth of the World 
Wide Web and its effect on the society have been as- 
tonishing, spreading from the research in high-energy 
physics into other scientific disciplines, academe in gen- 
eral, commerce, entertainment, politics and almost any- 
where where communication serves a purpose. Inno- 
vation has widened the possibilities for communication. 
Blogs, wikis and social bookmark tools allow the imme- 
diacy of conversation, while the potential of multimedia 
and interactivity is vast. The reason for this immediate 
success is the fact that no specific skills are needed for 
participating. In the so-called Web 2.0 (|0'Reillvl . I2005D 
users acquire a completely new role: not only information 
seekers and consumers, but information architects, coop- 
erate in shaping the way in which knowledge is struc- 
tured and organized, driven by the notion of meaning 
and semantics. In this perspective the web is acquir- 
ing the status of a platform for social computing, able 
to coordinate and exploit the cognitive abilities of the 
users for a given task. One striking example is give n 
by a series of web games ( von Ahn and Dabbishl . |2004| ). 
where pairs of players are required to coordinate the as- 
signment of shared labels to pictures. As a side effect 
these games provide a categorization of the images con- 
tent, an extraordinary difficult task for artificial vision 
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systems. More generally, the idea that the individual, 
selfish activity of users on the web can possess very use- 
ful side effects, is far more general than the example 
cited. The techniques to profit from such an unprece- 
dented opportunity are, however, far from trivial. Spe- 
cific technical and theoretical tools need to be developed 
in order to take advantage of such a huge quantity of 
data and to extr act from this noisy source solid and us- 
able information! Arrowl l2003t iHuberman and Adamid . 
I2OO4I ). Such tools should explicitly consider how users 
interact on the web, how they manage the continuous 
flow of data they receive (see Sec. IVIlI]) . and, ultimately, 
what are the basic mechanisms involved in their brain 
activity. In this sense, it is likely that the new social 
platforms appearing on the web, could rapidly become 
a very interesting laboratory for social sciences. In par- 
ticular we expect the web to have a strong impact on 
the studies of opinion formation, political and cultural 
trends, globalization patterns, consumers behavior, mar- 
keting strategies. 



2. Language and communication systems 

Language dynamics is a promising field which encom- 
passes a broader range of applications with respe ct to 
what described in Sec, rvl ijLoreto and SteeT^ . I2007D . In 
many biological, technological and social systems, a cru- 
cial problem is that of the communication among the dif- 
ferent components, i.e., the elementary units of the sys- 
tems. The agents interact among themselves and with 
the environment in a sensorial and non-symbolic way, 
their communication system not being predetermined nor 
fixed from a global entity. The communication system 
emerges spontaneously as a result of the interactions of 
the agents and it could change continuously due to the 
mutations occurring in the agents, in their objectives 
as well as in the environment. An important question 
concerns how conventions are established, how commu- 
nication arises, what kind of communication systems are 
possible and what are the prerequisites for such an emer- 
gence to occur. In this perspective the emergence of a 
common vocabulary only represents a first stage while it 
is interesting to investigate the emergence of higher forms 
of agreement, e.g., compositionality, categories, syntactic 
or grammatical structures. It is clear how important it 
would be to cast a theoretical framework where all these 
problems could be defined, formalized and solved. That 
would be a major input for the comprehension of many 
social phenomena as well as for devising new technologi- 
cal instruments. 



dividuals, their interests, behavior etc.. In this per- 
spective, the discovery of communities, i.e., the identi- 
fication of more homogeneous groups of individuals, is 
a major challenge. In this context, one has to distin- 
guish the communities as typically intended in social 
network analysis (SNA) (iFreemani . 12004 IScottl . l2000t 
IWasserman and Fausti Il994 ) from a broader definition 
of communities. In SNA one defines communities over 
a communication relationship between the users, e.g. if 
they regularly exchange e-mails or talk to each other. In a 
more general context, for e.g. providing recommendation 
strategies, one is more interested in finding communities 
of users with homogeneous interests and behavior. Such 
homogeneity is independent of contacts between the users 
although in most cases there will be at least a partial 
overlap between communities defined by the user con- 
tacts and those by common interests and behavior. Two 
important areas of research can be identified. On the one 
hand, there is the question of which observable features 
in the available data structures are best suited for infer- 
ring relationships between individuals or users. Selecting 
a fe ature affects the method u sed to detect communi- 
ties ( Girvan and Newmanl . [2003) . which may be different 
if one operates in the context of recommendation sys- 
tems or in the context of semantic networks. On the 
other hand important advances are foreseeable in the do- 
main of coevolution of dynamics and the underlying so- 
cial substrates. This topic is still in its infancy, despite 
the strong interdependence of dynamics and networks in 
virtually all real phenomena. Empirical data on these 
processes are becoming available: it is now possible to 
monit or in detail the ev olution of large scale social sys- 
tems (jPalla et a/.l . l2007l) . 
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3. Evolution of social networks 

As real and online social systems grow ever larger, 
their analysis becomes more complicated, due to their 
intrinsic dynamic nature, the heterogeneity of the in- 
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